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the  Corps  of  Engineers,  water  depths  are  frequently  shallow  (in  the  sense  of  wind-wave  be¬ 
havior)  and  conditions  are  often  unsteady  and  nonuniform.  From  a  survev  of  recent  theoret¬ 
ical  considerations  and  associated  field  experiments,  it  is  evident  that  drag  coefficients 
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•in  shallow  and  nearshore  waters  depend  on  a  broader  suite  of  parameters  than  in  the  open 
ocean  case.  Among  these  are  more  detailed  characterizations  of  sea  state,  water  surface 
currents,  wind  direction,  water  depth,  and  distance  from  shore  as  well  as  buoyancy 
parameters. 

To  help  clarify  this  dependency,  a  two-part  experiment  has  been  conducted  wherein 
all  of  the  hypothetically  relevant  parameters  have  been  measured  ,-T'  The  experiments  were 
conducted  at  the  Coastal  Engineering  Research  Center's  Field  Research  Facility  (FRF)  on  the 
coast  of  the  Atlantic  Ocean  near  Duck,  North  Carolina,  where  advantage  could  be  taken  of 
low  land  relief,  reasonably  uniform  alongshore  topography  and  bathymetry,  and  a  pier,  600  m 
in  length,  which  could  act  as  a  stable,  over-water  instrument  platform. 

The  first  experiment  was  conducted  in  fall  1985  during  a  larger  scale  experiment 
known  as  DUCK85.  Its  purpose  was  to  provide  an  exploratory  survey  of  the  horizontal  vari¬ 
ability  of  winds  at  the  site  and  to  investigate  adequacy  of  various  platforms  for  detailed 
wind-stress  measurements.  Data  from  that  experiment  (which  includes  passage  of  Hurricane 
Gloria)  indicate  that  wind  differences  from  the  landward  end  of  the  pier  to  the  seaward  end 
are  finite,  varying  by  up  to  40  percent  in  speed,  but  highly  systematic.  Anemometers  on 
the  roof  of  the  FRF  building  and  near  the  pier  deck  at  the  seaward  end  of  the  pier  indi¬ 
cated  substantial  flow  distortion  due  to  the  proximity  of  wind-blocking  structures.  From 
this,  it  is  concluded  that  extreme  caution  is  necessary  when  mounting  anemometers  to  avoid 
platform  distortion  effects  and  that  careful  interpretation  is  necessary  when  using  beach- 
based  anemometers  to  characterize  winds  over  water. 

The  second  experiment  was  conducted  in  fall  1986  during  a  very  large-scale  experi¬ 
ment  known  as  SUPERDUCK.  Its  purpose  was  to  provide  detailed  wind-stress  measurements. 
Coupled  with  measurements  made  by  other  investigators  during  SUPERDUCK  of  other  relevant 
parameters,  data  from  this  experiment  can  be  used  to  develop  and  test  various  models  of 
drag-coefficient  variability.  Experiment  rationale  and  design  are  described  here.  Experi¬ 
ment  results  and  analysis  will  be  given  in  a  subsequent  report. 
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NEARSHORE  WIND-STRESS  MEASUREMENTS:  BACKGROUND , 


PRELIMINARY  FIELD  WORK,  AND  EXPERIMENT  DESIGN 

PART  I:  INTRODUCTION 

1.  Wind  plays  a  very  important  role  in  the  dynamics  of  the  upper  ocean 
and  shallow-water  bodies.  Momentum  imparted  by  the  wind  to  the  water  surface 
gives  rise  to  waves  and  currents,  and  through  turbulence,  a  much  increased 
diffusive  capabi'ty  over  pure  molecular  processes.  In  the  deep  ocean,  the 
effects  of  wind  are  readily  detectable  to  depths  on  the  order  of  100  m.  In 
waters  shallower  than  this,  such  as  coastal  regions,  harbors,  entrance  chan¬ 
nels,  and  lakes,  wind  can  affect  the  dynamics  of  the  whole  water  column.  The 
interaction  of  the  water  with  solid  or  semisolid  boundaries  is  of  particular 
importance  to  those  who  dwell  near  water.  Wind-generated  waves  radiate  energy 
shoreward  to  beaches,  coastal  dwellings,  and  shore  protection  structures.  In 
high  winds,  the  wind-generated  water  levels  can  surge  over  land  and  cause 
flooding.  Wind-induced  turbulence  can  add  to  wave-  and  current-induced  tur¬ 
bulence  and  thereby  enhance  the  diffusive  character  of  the  flow  which,  in 
turn,  has  a  direct  impact  on  coastal  erosion,  alongshore  sediment  transport, 
and  redistribution  of  biological  nutrients  and  waterborne  contaminants. 

2,  Wind  forcing  is  not  the  only  source  of  energy  at  a  given  shallow- 
water  site.  Tides,  planetary  scale  waves,  lateral  processes  from  large-scale 
currents  driven  by  distant  winds  (such  as  the  Gulf  Stream) ,  waves  in  the  wind- 
wave  frequency  band  radiating  in  from  distant  storms  (swell),  and,  in  enclosed 
bodies,  seiching  can  all  contribute  to  nearshore  dynamics.  In  addition,  local 
areas  where  water  runoff  from  land  occurs,  as  in  estuaries  and  near  river  out¬ 
falls,  are  influenced  strongly  by  currents  and  turbulence  from  this  source. 

At  any  given  time,  fv.  cing  by  any  of  these  processes  may  be  present.  However, 
high-energy  water  motions,  i.e.,  those  which  can  do  work  (or  damage)  at  the 
highest  rate,  are  associated  with  high  wind  conditions.  In  a  carefully  de¬ 
signed  and  highly  instrumented  study  of  coastal  dynamics  in  Lake  Erie, 

Schwab  et  al.  (1984)  found  from  an  analysis  of  four  storms  that  wind  forcing 
was  a  dominant  source  of  energy  throughout  the  nearshore  zone.  This  is  impor¬ 
tant  because  it  means  that  wind  has  a  strong  influence  not  only  in  deep  and 
transitional  water,  where  active  wave  generation  is  classically  thought  to 
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occur,  but  also  in  shallow  water  where  waves  shoal  and  break  and  where  wave- 
induced  currents  are  often  considered  dominant.  Schwab  et  al.  (1984)  and 
Huberts  found  wind-driven  currents  to  be  comparable  with  wave-induced 

currents.  Hence,  a  complete  dynamic  description  of  nearshore  processes  must 
include  a  characterization  of  wind  input. 

3.  To  do  this  accurately  and  directly  is  not  a  straightforward  task. 

The  simplest  estimator  of  the  intensity  of  wind  forcing  is  a  measurement  of 
the  average  wind  speed  just  above  or  near  the  water  surface  being  forced.  It 
is  intuitive  that  as  the  mean  wind  speed  increases  the  force  imparted  to  the 
water  increases.  All  observations  indicate  that,  for  all  other  system  vari¬ 
ables  held  constant,  this  is  approximately  (but  only  approximately)  true. 

TIi is  type  of  reasoning  prompted  Taylor  (1916)  to  propose  that  the  mean  wind 

stress  i  with  units  of  force  per  area,  the  air  density  p  with  units  of 

mass  per  volume,  and  the  mean  wind  speed  U  with  units  of  length  per  time 

(and  aligned  with  the  stress,  both  being  taken  here  to  be  in  the  x-direction) 

2 

form  a  closed  system.*  Dimensional  analysis  then  requires  t/oU  =  constant 
■  CD  ,  where  CD  is  called  the  drag  coefficient.  This  result  is  called  the 
quadratic  drag  law  or  the  bulk  parameterization  of  wind  stress  and  is  usually 
written  in  the  form 

t  =  pCDU2  (1) 

4.  If  Equation  1  is  true,  then  is  constant  and  measurements  would 

provide  its  value.  Table  1  shows  ranges  of  given  by  several  authors  from 

measurements  made  over  water.  Some  of  the  references  cite  the  authors  own 
observations.  Others  are  compendia  of  several  experiments.  For  example, 
Garratt  (1977)  cites  23  references  and  Huang  et  al.  (1986)  cite  21  with  very 
little  redundancy  in  citation.  Amorocho  and  DeVries  (1980)  cite  virtually  the 
same  data  as  these  last  two,  but  their  paper  is  included  for  completeness. 
Clearly  there  is  doubt  about  the  constancy  of  .  The  ratio  of  maximum  to 

minimum  value  for  the  table  as  a  whole  is  about  32.  The  average  from  individ¬ 
ual  authors  is  about  6. 

5.  This  is  less  than  satisfactory  for  accurate  modeling.  If  wind 


*  For  convenience,  symbols  and  abbreviations  are  listed  in  the  Notation 
(Appendix  A) . 
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Table  1 


Panees  of 

Drag  Coefficient  Over  Water  Measured  or 

Cited  by  Various  Authors 

Source 

Minimum 

Maximum 

Hsu  (1972) 

1  . 1 

6 . 4 

SethuRaman  and  Raynor 
(1975) 

0.5 

3.0 

Garratt  (1977) 

0.6 

4.2 

SethuRaman  (1979) 

0.3 

2.1 

Amorocho  and  DeVries 
( 19801 

0.2 

4.5 

U’u  (1980) 

0.8 

2.3 

Large  and  Pond  (1981) 

0.5 

2.3 

Geernaert  and  P'atsaros 
(1986) 

1.3 

Geernaert,  Katsaros, 
and  Richter  (1986) 

0.7 

3.2 

Huang  et  al.  (1986) 

0.2 

4.2 

stress  is  as  important  as  other  effects  in  nearshore  processes  and  its  magni¬ 
tude  is  not  known  to  within  a  factor  of  6,  then  it  can  be  expected  that  model 
results  will  vary  similarly.  Relying  on  such  a  model  for  design  of  coastal 
protection  could  result  in  exhorbitantly  expensive  overkill  or  woefully  inade¬ 
quate  protection  depending  on  how  the  model  is  used  or  interpreted. 

6.  The  range  of  values  in  Table  1  suggests  and,  in  fact,  most  of  the 
authors  cited  there  propose,  at  if  a  model  represented  by  Equation  1  is  to 
be  used,  then  must  implicitly  contain  all  of  the  physics  not  explicit  in 

Equation  1.  The  necessary  physics  are  contained  in  the  fields  of  boundary- 
layer  turbulence  and  air-sea  interaction.  A  number  of  investigations  in 
recent  years  have  led  to  improved  characterization  of  C^  in  terms  of  other 
measurable  wind-field  and  sea-state  parameters.  Some  are  well  founded  in 
turbulence  theory;  others  are  more  conjectural  but  physically  reasonable. 

Part  II  of  this  report  contains  a  brief  summary  of  this  latest  work.  It  also 
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describes  the  shortcomings  of  our  present  understand i.ie  and  gives  justifica¬ 
tion  for  ongoing  field  experiments  to  refine  ar.d  test  arametric 
rel at icnships . 

It  is  not  proposed  that  bulk  parameterization  (i.e.,  drag  and  ex¬ 
change  coefficients)  be  abandoned.  They  are  the  most  reasonable  of  practical 
tools  available  since  they  rely  on  a  finite  number  of  specifically  defined 
mean  environmental  parameters  of  a  quality  which  can  be  measured  at  remote 
coastal  sites,  instrumented  buovs,  and  vessels  of  opportunity.  Many  of  the 
models  used  at  the  Coastal  Engineering  Research  Center  (CERC)  rely  on  bulk 
parameterization  for  wind  effects  in  coastal  flooding  models  (e.g.,  Sheng  and 
Butler  1982)  and  nearshore  wave  and  current  models  (e.g.,  Hubertz  1987).  Al¬ 
ternatives  are  either  to  deploy  wind-stress  measui ing  devices  on  ships,  buoys, 
and  coastal  sites  or  to  employ  a  large-scale  numerical  model  to  represent  the 
balance  of  forces  present  in  coastal  meteorology.  Neither  of  these  is  as 
practical  for  large-scale,  long-term  studies  as  the  bulk  method  although  both 
are  attractive  conceptually . 

8.  Direct  measurements  of  wind  stress  require  highly  accurate,  three- 
dimensional  anemometers  with  good  frequency  response  to  about  5  Hz.  Indirect 
methods  (other  than  the  bulk  method)  require  either  arrays  of  low-frequency- 
response  temperature,  humidity,  and  horizontal  wind-speed  sensors  or  single 
high-frequency  (20  Hz  to  1  kHz,  depending  on  method)  sensors  for  the  same 
variables.  Though  impractical  for  remote,  long-term,  unattended  deployment, 
such  instrumentation  and  methods  are  necessary  to  quantify  the  interrelation¬ 
ships  among  parameters  used  in  the  models  which  form  the  basis  of  the  bulk 
method.  Part  IV  of  this  report  gives  a  brief  description  of  currently  used 
methods,  instruments,  and  platform  requirements  and  limitations  tor  making 
wind-stress  measurements.  Part  V  of  this  report  describes  an  experiment  which 
tests  various  platforms  at  CERC's  Field  Research  Facility  for  adequacy  in 
making  wind-stress  measurements  by  various  methods. 

9.  Meteorological  models  such  as  the  one  described  by  Resio,  Vincent, 
and  Corson  (1982)  for  the  open  ocean  would  be  attractive,  if  applicable,  for 
use  in  coastal  regions.  However,  as  described  by  Hsu  (1970,  1979)  in  studies 
near  the  coast  of  Texas,  and  Resio  and  Vincent  (1977)  in  analysis  of  winds 
over  the  Great  Lakes,  conditions  are  frequently  highly  nonuniform  in  both 
cross-shore  and  vertical  directions.  Highly  variable  surface  roughness  from 
water  to  land  (and  possibly  over  the  water  itself,  see  Part  IT)  and  sharp 


contrast  -  in  surface  heat,  moisture,  and  momentum  fluxes  give  rise  to 
baruclinic  c i rcula t ions  famong  them  the  sea  breeze)  and  the  formation  of 
internal  boundary  layers  (Stunder  and  SethuRaman  1985),  the  effects  of  which 
are  noticeable  for  kilometres  to  tens  of  kilometres  landward  and  seaward  of 
the  coast.  Furthermore,  such  models  require  characterization  of  the  surface 
fluxes  of  mass  and  momentum.  Conventionally ,  this  is  done  in  the  same  way 
that  the  surface  hulk  coefficient  models  are  derived.  Any  limitation  in  our 
abilitv  to  characterize  the  bulk  coefficients  will  certainly  be  reflected  as 
inadequacies  in  a  large-scale  meteorological  model. 

10.  Hence,  it  seems  most  productive  to  confine  attention  to  the  lowest 
few  tens  of  metres  of  the  atmosphere,  where  mean  wind  speed,  temperature,  and 
humidity  measurements  are  most  often  made,  and  to  deduce  the  laws  that  relate 
these  quantities  to  their  respective  fluxes.  In  consideration  of  observations 
and  hypotheses  given  in  Parts  II  and  III  of  this  report,  and  measurement 
methods  described  in  Part  IV,  an  experiment  design  is  described  in  Part  V 
which  should  add  significantly  to  our  ability  to  formulate  these  lavs. 
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PART  II: 


RECENT  WORK  IN'  BULK  COEFFICIENT  FORMULATION 


11.  ” '  'ic  count  for  the  broad  range  of  drag  coefficients  shown  in 
TabLe  I,  it  is  assumed  either  that  there  is  substantial  instrument  ard  experi¬ 
ment  error  or  that  the  drag  coefficient  depends  on  other  parameters.  While 
early  methods  of  wind-stress  measurement  were  often  inaccurate,  modern  methods 
and  instrumentation  are  substantially  better.  The  broad  range  of  investiga¬ 
tions  represented  by  Table  l,  some  of  which  are  quite  recent,  suggests  that  a 
parametric  dependency  of  drag  coefficients  is  legitimate. 

12.  Several  investigators  have  compiled  data  from  other  sources  to  ex¬ 
amine  behavior  of  drag  coefficients.  One  of  the  most  frequently  cited  papers 
is  that  by  Garratt  (1977).  In  that  paper,  a  large  number  of  investigations, 
selected  for  reasonable  accuracy  and  lack  of  atmospheric  stabilitv  effects,  is 
analyzed  for  a  dependence  of  the  drag  coefficient  on  the  wind  speed  at  a  10-m 
elevation  above  the  sea  surface.  Though  the  scatter  of  results  is  high,  there 
is  a  clear  correlation.  Furthermore,  Garratt  (1977)  was  able  to  show  the 
major  trend  of  the  data  with  a  simple  atmospheric-surface-layer  model,  a  major 
step  in  accounting  for  some  of  the  variability  of  the  drag  coefficient. 

13.  The  model  described  by  Garratt  (1977)  and  modifications  to  it  will 
be  described  below.  In  the  paragraphs  that  follow  some  of  the  equations  and 
scaling  laws  of  turbulence  theory  will  be  used,  often  without  pretext  or  cita¬ 
tion.  Substantial  background  material  on  the  governing  equations  and  common 
scaling  arguments  can  be  found  in  such  classic  texts  as  those  by  Hinze  (1959), 
Phillips  (19691,  Kraus  (1972),  and  Tennekes  and  Lumlev  (1972).  Equations  will 
only  be  used  as  needed  for  clarity.  It  is  noted  that  because  of  its  complex¬ 
ity,  turbulence  is  a  highly  empirical  field  with  many  of  its  hypotheses  based 
on  scaling  arguments  (dimensional  analysis)  and  experiments.  Some  of  these 
are  well  founded,  others  less  so.  An  attempt  will  be  made  here  to  distinguish 
amone  them. 


Basic  Concepts 

14.  The  environment  for  which  the  drag  coefficient  formulae  are  de¬ 
signed  is  considered  to  be  horizontally  homogeneous,  turbulent  air  flow  over 
water  with  deviations  from  horizontal  homogeneity  assumed  negligible.  The 
wind  stress  to  be  defined  is  not  the  true  Newtonian  viscous  stress  but  a 
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vertical  advection  of  fluctuating  horizontal  momentum  known  as  a  Reynolds 
stress.  In  a  coordinate  system  with  the  x-axis  aligned  with  the  mean  wind  and 
the  vertical  axis  positive  upward  from  the  sea  surface,  the  total  velocitv  in 
the  x-direction  u  can  be  thought  of  as  having  mean  U  and  fluctuating  u' 
parts  where 

T 

1  a 

u  -  f  u  dt  (?) 

a  0 

T  is  an  appropriate  averaging  time  (see  Part  III)  and  integration  is  with 
respect  to  time.  The  total  vertical  velocity  w  has  no  mean  value  in  this 
idealized  flow,  but  only  a  fluctuating  part  w'  . 

15.  The  total  momentum  per  unit  volume  in  the  x-direction  is  ou  =  D 
(U  +  u')  where  o  is  the  air  density  (mass  per  unit  volume).  The  mean  ver¬ 
tical  advection  of  this  bv  the  fluctuating  (turbulent)  vertical  velocity  is 
the  correlation  between  ou  and  w'  .  This  is  usually  found  to  be  negative 
so  the  negative  correlation  is  considered  a  positive  stress  x  as 

T 

i  ^  _ 

x  =  —  —  C  Puw’  dt  =  -  D  u’w’  (3) 

3  0 

where  the  overbar  represents  the  time  average,  U  does  not  contribute  because 
it  is  constant  for  this  period  by  Equation  2,  and  fluctuating  parts  of  the 
density  field  are  assumed  negligible  through  the  Boussinesq  approximation. 
Equation  3  is  considered  a  stress  because  it  is  a  flux  of  momentum,  albeit 
advective  instead  of  molecular.  It  is  positive  because,  on  the  average,  when 
ou'  is  positive,  w'  is  negative.  This  represents  a  downward  momentum  flux 
which  tends  to  accelerate  fluid  layers  below,  including  the  water  surface. 

When  pu'  is  negative,  w'  tends  to  be  positive.  This  tends  to  decelerate 
fluid  layers  above  (including  the  mean  air  flow), 

16.  The  mean  flow  U  is  a  strong  function  of  z  in  the  vicinity  of 
the  water  surface.  At  the  water  surface,  the  air  molecules  are  coupled  to 
water  molecules  by  molecular  viscosity  such  that  the  mean  air  velocity  equals 
the  mean  water  velocity.  Well  above  the  water  surface  the  air  is  forced  by 


very  large-scale  pressure  ind  gravitational  forces  and  often  moves  much  faster 
than  the  wate-.  In  between,  the  mean  flow  is  sheared,  i.e.,  varies  with 
respect  to  the  vertical  coordinate  z  .  Fluctuations  in  vertical  velocity 
tend  to  carry  faster  air  down  and  slower  air  up,  which  suggests  a  coupling 
between  the  mean  shear  3U/?z  and  the  stress  given  by  Equation  3.  Since  u' 
and  w'  are  correlated,  organized  rotary  motions,  or  eddies,  are  suggested. 
The  eddies  cannot  penetrate  the  water  surface,  so  a  maximum  size  for  eddies  in 
the  vicinity  of  the  water  surface  is  likely  to  be  the  distance  z  to  the 
water  surface.  If  x  ,  0  ,  and  z  form  a  closed  set,  then  by  dimensional 
analysis 


z 

u 


•k 


3U 

?z 


constant 


_1 

< 


(4) 


where 

I/O 

u*  =  (t/ p)  “  friction  velocity 
<  =  von  Karman  constant 

17.  The  model  given  by  Equation  4  will  not  work  for  air  flow  in  the 
troughs  between  wave  crests  because  the  flow  is  not  horizontally  homogeneous 
and  becomes  extremely  complex.  This  is  shown  theoretically  by  Chalikov  (1986) 
and  Jacobs  (1987)  and  in  laboratory  work  by  Papadimitrakis ,  Hsu,  and  Street 
(1984)  and  Papadimitrakis,  Hsu,  and  Wu  (1986).  In  analogy  with  flow  over 
solid  roughness  elements  (Schlichting  1979),  it  is  usually  assumed  that  flow 
distortion  and  stress  redistribution  by  wave  surfaces  decays  by  about  three 
wave  heights  above  the  surface.  Equation  4  should  hold  for  some  distance 
above  that,  i.e.,  the  air  acts  as  if  it  is  flowing  over  a  horizontally  uniform 
surface.  This  provides  one  constraint  on  the  region  of  applicability  of  Equa¬ 
tion  4,  i.e.,  z  >  3H  ,  where  H  is  a  characteristic  wave  height. 

s  s 

18.  Equation  4  provides  a  very  powerful  tool  for  air-sea  interaction 

2 

studies  since  it  allows  a  relation  between  the  wind  stress  pu^  and  the  mean 
velocity  U  at  some  distance  above  the  water  surface.  The  point  of  this 
discussion  is  that  by  the  continuum  hypothesis  the  flux  is  continuous  through 
the  level  z  at  which  measurements  are  made.  If  the  flow  is  horizontally 
homogeneous,  then  the  only  way  momentum  can  get  to  the  water  surface  from  the 
air  is  vertically  through  elevation  z  .  Sufficient  knowledge  of  this  and  the 
water  response  is  all  that  is  needed  to  form  a  closed  model.  Detailed 
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knowledge  of  the  air  flow  in  and  near  the  wave  troughs  is  then  unnecessarv 
which  is  an  enormous  simplification  of  the  problem. 

19.  The  bulk  of  recent  work  focuses  on  two  main  topics.  One  is  the 
specification  of  the  constants  of  integration  in  the  solution  of  Equation  4. 
The  other  is  accounting  for  atmospheric  stability,  i.e.,  variations  of  air 
density  due  to  fluxes  of  heat  and  water  vapor,  which  is  known  to  modify  Equa¬ 
tion  4  from  studies  of  steady,  horizontally  uniform  flow  over  land.  These  are 
discussed  in  turn  and  are  followed  by  considerations  of  unsteadiness  and  hori¬ 
zontal  nonuniformity . 

20.  If  the  wind  stress  is  constant  or  a  weak  function  of  distance  from 

the  sea  surface,  then  Equation  4  can  be  integrated  simply  between  two  levels 

of  flow.  Long  (1981)  in  a  similarity  solution  of  an  idealized,  neutrally 

stable  planetary  boundary  layer  model  given  by  Businger  and  Arya  (1974)  showed 

that  the  stress  is  approximately  linear  in  the  lowest  part  of  the  boundary 

laver  with  a  slope  3t/3z  on  the  order  of  -6pu.f  ,  where  f  is  the 

*  c  c 

Coriolis  parameter  (»  2 S3  sin  A  ,  where  (2  is  the  radian  rotation  rate  of 

Earth  and  A  is  latitude! .  The  height  at  which  t  is  still  90  percent  of 

its  boundarv  value  is  then  on  the  order  of  u^/bOf  .  For  a  standard  anemom- 

*  c 

eter  height  of  10  m,  this  imposes  a  constraint  that  u.  >  60f  *  10  m 

*  c 

=  0.06  m/sec  or  U(10  m)  >  1  m/sec  based  on  the  drag  coefficient  curve  of 
Garratt  (1977).  The  corresponding  constraint  for  an  anemometer  height  of  20  m 
is,  at  most,  U(20  m)  _>  2  m/sec.  Since  wind  speeds  of  interest  are  usually 
much  greater  than  this,  i.e.,  5  m/sec  and  higher,  it  is  reasonable  to  treat  t 
as  approximately  constant  in  the  lower  few  tens  of  metres  of  the  atmosphere. 

As  a  result,  this  region  is  often  called  the  "constant  stress"  or  "constant 
flux"  layer  of  the  atmosphere  with  the  understanding  that  some  caution  is  to 
be  exercised  at  very  low  wind  speeds. 

21.  Equation  4  is  then  expected  to  apply  for  z  £  3Hg  ,  as  noted  previ¬ 
ously,  and  for  z  £  u^/60fc  .  An  observation  height  of  10  m  should  be  ade¬ 
quate  for  wave  heights  less  than  about  3  m.  For  seas  to  6  m,  an  elevation  of 
about  20  m  is  necessary.  Violations  of  these  constraints,  e.g.,  10-m  observa¬ 
tions  during  hurricanes  where  seas  exceed  3  m,  could  be  expected  to  yield 
suspicious  results. 

Logarithmic  Profiles  and  Limits  of  Applicability 

22.  With  t  ,  and  therefore  u^  ,  assumed  constant  in  the  atmospheric 


11 


surface  layer.  Equation  4  can  be  integrated  between  two  arbitrary  levels  of 
the  flow  such  that 

U* 

l'(z)  -  lT  =  —  In  —  (5) 

O  K  Z 

o 

where  lT  and  z^  are,  mathematically,  constants  of  integration.  For  gen¬ 
eral  application,  they  must  be  specified  as  constants,  functions  of  other  flow 
parameters,  or  specific  velocity  and  level  within  the  region  where  Equation  5 
applies.  From  Equations  1,  3,  and  5,  along  with  the  definition  of  u^  fol¬ 
lowing  Equation  4,  a  drag  coefficient  expression  can  be  written  as 


where,  as  discussed  below,  the  inclusion  of  modifies  the  form  given  in 

Equation  1. 

23.  From  Equation  6  it  is  seen  that  a  drag  coefficient  is  a  function  of 
z  in  general.  Setting  a  standard  anemometer  height  of,  say,  z  ■  z  -  10  m 
compensates  for  this  dependency.  Equation  6  is  often  used  to  correct  for  data 
collected  at  heights  other  than  z  *  10  n  to  standardize  comparisons. 

24.  In  standard  form,  C_  depends  on  ic  ,  U  ,  and  z  .  The  value 

D  oo 

of  <  is  usually  taken  to  be  0.41.  Frenzen  and  Hart  (1983)  give  a  list  of 
recent  estimates  based  on  atmospheric  studies.  Nowell  and  Long  (1983)  in  a 
detailed  and  comprehensive  analysis  of  pipe-flow  data,  find  pc  ■  0.408  t  0.004 
at  95-percent  confidence.  However,  in  a  very  well  executed  atmospheric  exper¬ 
iment  over  land  (Haugen,  Kaimal,  and  Bradley  1971),  the  data  analysis  by 
Businger  et  al.  (1971)  found  <  *  0.35  .  This  value  is  given  some  weight  be¬ 
cause  of  the  high  quality  of  the  experiment.  On  the  other  hand,  a  reanalvsis 
of  the  same  data  by  Wierenga  (1980),  who  took  flow  blockage  effects  into  ac¬ 
count  (see  Part  IV) ,  found  <  =  0.41  .  The  difference  is  important  since  < 
is  squared  in  Equation  6.  A  17-percent  difference  in  <  yields  about  a 
34-percent  difference  in  and  any  stress  computed  using  .  Note  that 

data  scatter  reported  by  an  investigator  who  uses  a  single  value  of  <  to 
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estimate  C  is  not  affected  bv  this,  but  the  application  of  a  model  of  C 

IJ  -  r  r  d 

to  estimate  stress  is  affected.  The  concensus  of  investigations  shown  bv 
Frenzen  and  Hart  ''1^83)  supports  the  higher  value  of  k  .  While  the  debate 
continues  fsee  Hogstrom  (1985)  and  a  reply  by  Telford  and  Businger  f  1986) ,  for 
example)  it  is  assumed  here  that  <  =  0.41  . 

25.  In  studies  of  fluid  flow  over  solid,  impermeable  boundaries.  Equa¬ 
tion  6  is  generallv  applied  with  U  =0  and  z  ,  called  the  surface  roueh- 

o  o 

ness  scale,  is  the  level  at  which  an  extrapolation  of  the  logarithmic  profile 
toward  the  boundary  gives  U(z  )  =  0  .  As  discussed  above  concerning  the  com¬ 
plex  flow  in  the  region  very  near  a  water-wave  surface,  extrapolation  of  a 
logarithmic  profile  to  a  region  very  near  a  solid  surface  does  not  provide 
legitimate  velocity  estimates  in  this  so-called  sublayer  region.  Hence,  zq 
is  an  artifice  both  of  integration  of  Equation  4  and  of  application  of  the 
no-slip  boundary  condition.  By  the  continuum  hypothesis,  however,  the  mean 
velocity  must  be  continuous  into  the  sublayer  such  that  z^  can  be  considered 
a  measure  of  the  elevation  at  which  the  mean  velocity  in  the  sublave’-  matches 
the  velocity  in  the  overlying  logarithmic  layer.  Since  the  matching  height  is 
the  lowest  region  of  applicability  of  the  logarithmic  region,  the  turbulent 
stress,  i.e.,  Reynolds  momentum  flux,  must  also  be  continuous  at  this  height. 

26.  Within  the  sublayer,  the  flow  can  be  quite  complicated  depending  on 

the  physical  geometry  of  the  solid  boundary.  If  the  boundary  is  smooth 

enough,  then  momentum  transfer  at  the  boundary  is  purely  by  Newtonian  viscous 

shear  stress.  In  this  case  x  =  u3U/3z  ,  where  u  is  the  dynamic  viscosity, 

2 

or  in  equivalent  form,  u^  ■  v3U/3z  ,  where  v  =  u/p  is  the  kinematic  viscos¬ 
ity  (with  dimensions  of  length  squared  per  time) .  The  velocity  scale  of  this 

flow  is  u^  and  the  length  scale  is  v/u^  .  The  mean  flow  itself  is 

2 

U(z)  =  u“z/v  .  This  profile  has  been  found  empirically  (Schlichting  1979)  to 
match  an  overlying  logarithmic  profile  at  z  =  10  v/u^  ,  such  that,  from  this, 
z  =  v/9u.  .  This  flow  is  called  "dynamically  smooth,"  and  occurs  if  the 
solid  boundary  has  no  deformations  or  protuberances  with  length  scales  larger 
than  about  5  v/u^  .  Note  that  the  total  vertical  flux  of  momentum  is 
considered  constant  through  the  sublayer.  The  turbulent  advection  of  momen¬ 
tum,  which  scales  with  z  ,  decreases  toward  the  boundary,  but  is  compensated 

2 

by  viscous  stress  which  increases  to  pu*  at  the  boundary. 

27.  If  the  physical  boundary  contains  roughness  elements  of  length 
scale  k^  >>  v/u*  ,  then  momentum  is  transferred  through  the  sublayer  to  the 
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boundarv  net  only  by  viscous  shearing  stress  (skin  friction)  but  also  by  the 

net  dynamic  pressure  (normal  or  form  drag)  forces  induced  by  the  fluid  as  it 

is  diverted  around  and  over  roughness  elements.  In  this  case,  known  as 

"dvnamically  rough,"  the  mean  flow  is  very  complex  and,  in  general,  has  no 

simple  solution.  The  length  scale  of  this  flow  is  dominated  by  the  size  of 

the  roughness  -lements  k  such  that  for  boundaries  with  roughness  of  a 

single  scale  it  can  be  expected  that  zq  and  the  sublayer  matching  height 

both  scale  with  k  .  This  was  verified  empirically  (Schlichting  1979)  for 

flow  over  close-packed  sand  grains  of  diameter  kR  .  It  was  found  that 

z  =  k.,/30  and  the  sublayer  matching  height  was  about  3  k  .  While  this 
o  R  K 

scale  for  z  is  often  assumed  for  dynamically  rough  flows  in  general,  it  is 
o 

important  to  note  that  over  surfaces  with  multiple  length  scales  the  roughness 
height  alone  is  not  sufficient.  For  instance,  Schlichting  (1979)  finds 
k  /200  5  z  ^  k  / 1 0  for  flow  over  elements  of  height  kB  but  of  different 
separation  and  pattern  scales.  That  is,  length  scales  of  surface  geometry 
appear  to  dominate  the  definition  of  z^  ,  but,  for  surfaces  of  multiple 
scales,  in  complex  combination.  The  total  momentum  flux  into  the  sublayer  at 
the  matching  height  is  t  =  pu^  as  noted  above.  Within  the  sublayer,  the 
division  of  this  flow  into  form  drag  and  skin  friction  is  a  complicated  func¬ 
tion  of  local  surface  geometry,  local  flow  scales,  and  viscosity.  For  total 
force  on  a  rigid  surface  (like  the  inside  of  a  pipe)  it  is  not  necessary  to 
know  the  detailed  flow  in  the  sublayer  to  relate  the  stress  to  the  velocity 
above  the  sublayer,  if  the  net  mean  roughness  scale  zq  is  known. 

28.  The  discussion  above  is  necessary  as  background  because  the  success 
of  the  reasoning  used  for  flow  over  solid  boundaries  suggests  that  similar 
reasoning  can  be  used  for  flow  over  a  water  surface.  Such  reasoning  is  bound 
to  be  incomplete,  however,  because  a  water  surface  differs  from  a  rigid  sur¬ 
face  in  that  it  deforms  under  the  action  of  wind  stress.  In  geophysical 
flows,  this  deformation  is  characterized  by  a  broad  range  of  length  scales. 

In  dynamic  conditions,  vertical  length  scales  vary  from  millimetres  (for  the 
amplitudes  of  capillary  waves)  to  10  m  in  storm  conditions  at  sea,  shoaling  of 
long  waves,  or  standing  tidal  motions.  Horizontal  length  scales  vary  from 
centimetres  for  capillary  waves  to  hundreds  of  kilometres  for  tidal  and 
planatary  scale  waves.  In  calm  conditions,  length  scales  are  very  small, 
displacements  being  on  the  order  of  molecular  mean  free  paths  as  governed  by 
surface  tension  and  thermodynamic  constraints  on  intermolecular  forces. 
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Surface  Roughness 


2°.  Relating  water  surface  displacement  and  motion  to  the  mean  and 
fluctuating  parts  of  the  driving  wind  is,  of  course,  the  very  heart  of  air-sea 
interaction.  A  primary  objective  in  this  field  is  to  determine  which  scale  or 
combination  of  scales  of  the  sea  surface  characterizes  the  surface  roughness 
scale  zq  in  Equation  5  so  that  from  Equation  6  can  be  determined. 

30.  Remember  that  Equation  5  is  assumed  legitimate  in  steady,  horizon¬ 
tally  uniform,  neutrally  stratified  flow.  As  for  flow  over  solid  surfaces, 
the  sea  surface  is  considered  uniform  if  statistics  of  surface  displacements 
and  motions  do  not  vary  substantially  in  any  horizontal  direction.  Further¬ 
more,  if  the  sea  is  considered  to  consist  of  wind-driven  waves  only,  and  these 
waves  are  assumed  linear  (i.e.,  all  motions  are  purely  cyclic),  then  there  is 
no  mean  motion  at  the  air-sea  interface.  In  this  case  Uq  =  0  in  Equations  5 

and  6  and  C  depends  entirelv  on  z  for  a  fixed  reference  elevation  z 
D  r  o  r 

and  reference  wind  speed  U  =  U(z  )  . 

r  r 

31.  Equation  6  can  be  solved  for  z  which  results  in 


z 

o 


z  e 
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z  e 
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(7) 


From  a  collection  of  observations  of  wind  stress  and  wind  speed  made  prior  to 

1970,  Kitaigorodskii  (1973)  noted  that,  apparently,  z  varies  by  over  seven 

-7  -1  ° 

orders  of  magnitude  from  about  10  m  to  about  10  m  for  virtually  the  same 
value  of  wind  stress.  He  suggested  that  much  of  the  scatter  was  due  to  un¬ 
certainty  in  investigative  method  and  cited  the  need  for  proper  parameteriza¬ 
tion  of  z 

o 

32.  A  simple  parameterization  of  surface  roughness  was  proposed  by 
Charnock  (1955).  His  argument  was,  in  essence,  that  gravity  waves  (i.e., 
motions  whose  restoring  force  is  gravity,  characterized  by  gravitational  ac¬ 
celeration  g  )  being  driven  by  wind  stress  'with  principal  varying  part  being 
uA  since  P  is  more  nearly  constant)  might  be  expected  to  have  a  surface 
roughness  dictated  by  the  two  parameters  u^  and  g  .  The  only  possible 
dimensionally  correct  form  is 
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where  '  -  c  universal  constant.  From  his  own  observations  in  a  16-m-deep 
reservoir  with  lateral  dimensions  on  the  order  of  1  km,  Charnock  f 1^55)  esti¬ 
mated  a  =  0.0067  . 

33.  Following  a  slightly  different  reasoning  but  arriving  at  the  form 
given  by  Equation  8,  Kitaigorodskii  (1973)  treated  his  collection  of  scattered 
(pre-1970)  data  as  samples  from  a  random  process.  By  averaging  zq  estimates 
within  bands  of  ,  he  deduced  a  mean  a  =  0.035  with  80-percent  confidence 
intervals  spanning  about  three  orders  of  magnitude  in  z^  .  This  a  differs 
substantially  from  that  given  by  Charnock  (1955). 

34.  A  third  estimate  of  a  is  given  by  Garratt  (1977),  mentioned  pre¬ 
viously.  Citing  experiments  from  the  late  1960's  to  the  mid-1970's  (and, 
hence,  largely  independent  of  data  cited  by  Kitaigorodskii  (1973)),  he  found 
a  =  0.0144  .  This  differs  substantially  from  both  of  the  estimates  given 
above.  His  data  are  grouped  by  reference  wind  speed  U  before  averaging. 

His  estimate  of  confidence  is  the  standard  deviation  of  resulting  esti¬ 

mates.  Use  of  the  two  a  estimates  given  above  yields  C^  estimates  within 
Garratt 's  error  bars  for  wind  speeds  below  about  20  m/sec.  Above  20  m/sec, 
they  lie  outside  his  error  bars.  However,  his  data  are  sparse  at  speeds  above 
20  m/sec  so  the  alternate  curves  might  not  exceed  a  90-percent  confidence 
interval  (which  is  not  given),  and  so  may  not  be  significantly  different. 

35.  At  =  40  m/sec,  C^  varies  by  a  factor  of  about  2  for  a 

varying  from  0.0067  to  0.035.  Hence,  the  difference  is  important  in  high 

winds.  Note  that  though  a  varies  by  a  factor  of  about  5  in  this  comparison, 

the  difference  in  C  is  less.  This  is  because  it  is  the  logarithm  of  z  , 

D  o 

and,  hence,  a  ,  that  enters  Equation  5  so  that  C^  is  less  sensitive  to 
large  variations  in  estimates  of  a  .  At  *  10  m/sec,  C^  varies  by  about 
50  percent  for  the  two  a's  cited.  Garratt 's  (1977)  results  lie  between 
these  cases. 

36.  A  fourth  work,  Wu  (1980),  discusses  values  of  a  which  cover  the 
middle  part  of  the  range  given  above.  Citing  his  own  research  entirely,  Wu 
(1980)  states  that  he  has  found  a  «  0.012  from  laboratory  work  alone, 

a  =  0.0156  as  a  mean  value  joining  laboratory  and  field  data,  a  =  0.0185 
for  data  excluding  reference  winds  less  than  about  5  m/sec  (for  which  z^  may 
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be  affected  by  viscosity  or  surface  tension),  and  a  =  0.027  based  on  scaling 

arguments  for  wind  flow  separation  over  the  crests  of  oceanic  waves  of  lew 

phase  speed.  He  concludes  that  a  general  expression  for  z  should  include 

o 

parameters  other  than  u^  and  g  ,  and  notes  the  need  for  further  research. 

He  suggests  i  =  0.0185  for  the  z^  regime  he  believes  is  most  nearlv  gov¬ 
erned  by  u^  and  g  .  This  differs  by  about  25  percent  from  the  a  =  0.0144 
found  by  Garratt  (1977). 

37.  Both  Kitaigorodskii  (1973)  and  Garratt  (1977)  averaged  results  from 
many  experiments  to  obtain  statistically  meaningful  estimates  of  a  from  data 
with  much  scatter.  To  test  the  hypothesis  of  a  constant  a  ,  Amorocho  and 
DeVries  (1980)  conducted  wind-stress  measurements  in  a  wind  tunnel  and  over  an 
aquaduct.  They  report  a  broad  range  of  a  from  about  0.001  to  1.0  in  their 
wind  tunnel  and  from  about  0.0001  to  1.0  for  the  aquaduct.  For  a  given  wind 
speed,  the  scatter  in  the  data  is  large,  being  about  one  order  of  magnitude  in 
their  lab  work  and  about  two  orders  of  magnitude  in  their  field  work. 

Amorocho  and  DeVries  (1980)  use  these  results  to  argue  in  favor  of  a  simple 
three-range  drag  law  wherein  is  constant  within  each  of  two  ranges  of 

,  and  a  transition  range  joins  them.  However,  the  scatter  of  measure¬ 

ments  remains  high.  One  interpretation  of  their  results  is  that  Charnock's 
model  does  not  include  all  of  the  parameters  necessary  to  define  zq  in 
general . 

38.  A  more  straightforward  approach  to  the  definition  of  z  was  de- 

o 

scribed  by  Kitaigorodskii  (1973).  He  tested  the  hypothesis  that  wave  dis¬ 
placements  act  solely  as  roughness  elements.  In  any  sea  of  interest,  a  typi¬ 
cal  wave  height  is  >>  v/u^  so  the  sea  surface  is  everywhere  dynamically 
rough  in  this  scheme.  The  roughness  parameter  zq  should  then  be  propor¬ 
tional  to  some  measure  of  surface  displacement.  Kitaigorodskii  (1973)  chose 
sea-surface  variance  as  such  a  measure.  However,  he  found  no  general  in¬ 

dications  of  correlation  between  z^  measured  from  wind  profile  data  and  a 

39.  This  is  a  reasonable  result  since,  as  Munk  (1955)  pointed  out, 

basing  zq  on  simple  displacement  would  suggest  that  a  1-m  wind-driven  sea 

would  have  the  same  roughness  as  a  1-m  tide.  Kitaigorodskii  (1973)  next  made 

the  hypothesis  that  surface  slope  is  important  since  it  is  a  measure  of  the 

extent  to  which  a  wavy  surface  presents  a  bluff  obstruction  to  wind  flow.  He 

estimated  a  mean  slope  as  a  /X  ,  where  X  is  the  mean  wavelength  deduced 

no  o  ? 

from  mean  wave  period  Tq  by  the  deepwater  dispersion  relation  Xq  =  gTV2n  . 


With  z  ,  "  ,  and  c  / A  as  variables  in  a  -et  (assumed  closed),  he  sought 

Q  r  Q 

a  correlation  between  z  /c  and  o  /X  but  found  none. 

on  no 

40.  Results  of  these  two  tests  do  not  necessarily  mean  that  surface 
displacement,  with  or  without  surface  slope,  is  not  important.  They  indicate 
that  these  variables  are  incomplete  descriptors  of  z^  ,  possibly  because  of 
the  smoothing  effect  of  using  mean  wave  properties  in  their  definitions  or 
because  such  definitions  treat  the  sea  surface  as  static,  whereas,  in  nature, 
it  is  dynamic. 

41.  Following  the  logic  of  this  statement,  Kitaigorodskii  (1973)  rea¬ 
soned  that  if  roughness  elements  are  mobile,  then  it  is  wind  speed  relative  to 
speed  of  roughness  elements  that  is  important  in  their  interaction.  For 
roughness  elements  as  waves  with  phase  speed  C  (and  still  assuming  no  mean 
velocity  at  the  water  surface)  the  logarithmic  profile  of  Equation  5  can  be 
expressed  relative  to  C  as 


Lr(z) 


(9) 


where  z  is  the  elevation  at  which  U(z  )  =  C  ,  Manipulation  of  Equation  9 

3.  3 

enables  transformation  from  a  coordinate  system  moving  at  speed  C  to  one 
that  is  fixed.  The  result  is 


uCz)  ,  m  -  Jc7u^  do) 

z  e 
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Comparison  with  Equation  5  shows  m  z ^  exp  (-icC/u^)  .  On  the  asymptote  of 
C/u^  -*  0  (i.e.,  static  roughness)  the  expression  z^  =  z^  is  recovered. 

42.  For  finite  C/u^  ,  Kitaigorodskii  (1973)  made  the  somewhat  tenuous 

assumption  that  z  for  a  monochromatic  wave  is  proportional  to  the  wave  am- 

a  2 
plitude.  If  a  spectrum  of  waves  is  present,  then  (I/2)za  is  the  variance  of 

the  sea  surface  in  a  narrow  band  of  frequencies  Au  centered  at  frequency  w 

(in  radians  per  second),  i.e., 
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where  S(-)  =  spectral  density.  Since  phase  speed  is  uniquely  related  to 
frequency  for  a  given  water  depth  d  (at  least  in  linear  wave  theory  where 
C  =  (g/w)  tanh  („d/\)),  then  the  expression  for  ,  above,  is  frequency 

dependent,  Fitaigorodskii  (1973)  makes  the  further  assumption  that  the  rough¬ 
ness  associated  with  each  spectral  band  is  a  linear  contributor  to  the  total 
roughness  such  that 


-2k  C(oj)/u^ 
2S(u)  e  Aw 


(12) 


For  a  continuous  spectrum  in  deep  water  (where  C  =  g/w) ,  integration  over  all 

frequencies  and  solution  for  k  yields 
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(13) 


Total  roughness  is  thus  a  weighted  sum  of  all  spectral  components.  It  is 
emphasized  that  Equation  13  is  derived  for  deep  water  and  tacitly  assumes  all 
waves  are  propagating  in  the  direction  of  the  wind.  For  kR  >>  v/u*  it  is 
assumed,  in  this  model,  that  z  is  proportional  to  k  . 

43.  The  estimate  of  roughness  given  by  Equation  13  depends  strongly  on 
the  definition  of  S(w)  used  in  evaluation  and  on  the  size  of  the  exponent  in 
the  weighting  function.  For  very  low  frequencies,  w  +■  0  ,  the  weighting 
function  prevents  much  contribution  from  the  spectrum.  The  most  emphasis  is 
given  to  high-frequency  waves.  To  evaluate  Equation  13  in  a  simple  way, 
Kitaigorodskii  utilized  Phillips'  (1969)  equilibrium  spectral  model  for  high- 
frequency  wind  waves.  This  takes  the  form 


S(gj) 
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where 


^  =  spectral  peak  frequency 

5  =  constant  (according  to  Phillips  (1969)) 

This  spectrum  represents  a  state  of  saturation  of  a  wave  field  wherein  energy 
added  from  wind  is  transferred  to  turbulence  through  wave  breaking  or  to  low- 
frequency  waves  through  nonlinear  interaction.  This  last  mechanism  tends  to 
lower  w  in  Equation  1 4  during  a  condition  where  a  wind-driven  sea  is 
building . 

44.  Kitaigorodskii  (1973)  identified  three  general  regimes  dependent  on 
the  variables  involved  when  Equation  14  is  substituted  in  Equation  13.  First 
is  the  case  wherein  a  high  wind  acts  on  a  relatively  undeveloped  sea.  In  this 
case,  uA  is  much  greater  than  the  phase  speeds  of  any  spectral  components 
(g/w) ,  such  that  the  exponential  weighting  function  in  Equation  13  is  essen¬ 
tially  unity  for  all  frequencies.  Roughness  elements  thus  act  as  if  they  are 

2  2  4 

rigid  and  k"  =  (l/2)8g“/uJ  ,  which  is  simply  the  variance  of  a  sea  surface 

K  O 

with  a  spectrum  given  by  Equation  14. 

45.  If  the  wind  continues  to  act,  then  the  spectrum  broadens,  w  is 

o 

reduced  and  the  phase  speed  at  the  spectral  peak  (Cq  =  g/u^)  becomes  more  com¬ 
parable  to  u^  .  In  this  second  regime,  Kitaigorodskii  treats  the  exponential 
weighting  function  in  Equation  13  as  a  slowly  varying  function  (relative  to 
w  ■*)  characterized  by  oj  =  .  In  this  case,  the  surface  roughness  length 

squared  is  the  spectral  variance  attenuated  by  the  exponential  function  eval- 

2  2  A 

uated  at  u  ■  u  .  Thus,  k_  =  (1/2)  8g  w  exp  (-2icg/u.w  ) 

^  O  K  O  O 

=  ( 1 /2 ) Sg^oj  exp  (~2kC  /u.)  .  This  is  an  important  theoretical  result  since 

it  is  one  of  the  first  to  show  that  sea-surface  dynamics  (i.e.,  moving  rough¬ 
ness  elements)  have  a  direct  influence  on  the  structure  of  flow  in  the  over- 
lying  atmosphere.  Furthermore,  it  shows  a  reduction  in  surface  roughness  for 
the  dynamic  case  as  compared  with  the  static  roughness  case.  This  cor¬ 
responds,  at  least  qualitatively,  with  observations. 

46.  The  third  regime  represents  the  somewhat  idealized  case  where  the 
sea  has  developed  to  the  point  where  w  ->■  0  .  Equation  13,  using  the  spec¬ 
trum  of  Equation  14,  can  then  be  evaluated  directly.  The  result  is 

, —  2  2 

kR  =  ■'38/2<  u^/g  .  Comparison  of  this  result  with  Equation  8  shows  that  this 
has  the  same  parametric  dependency  as  Charnock's  (1955)  z  .  If,  in  analogy 
with  flow  over  uniform  sand  grains,  zq  =  k^/30  »  then  zq  *  /36/60<  u^/g  > 
and  Charnock's  coefficient  is  a  =  /Tb/60k^  .  SethuRaman  (1979)  evaluated 
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this  expression,  using  <  =  0.35  (from  Businger  et  al.  1971)  and  8  =  0 . 0 1  I ' 
(as  suggested  by  Phillips  (1969)).  He  found  a  =  0.026  .  This  is  between  the 
values  i  =  . 0 1 14  (given  by  Garratt  (1977))  and  a  =  0.035  (deduced  by 

Kitaigorodskii  '!9~3)).  Remember  that  the  two  latter  values  were  found  by 
averaging  large  amounts  of  data  without  discriminating  regimes.  Variation  in 
j  is  quite  possible  due  to  variations  in  Phillips'  coefficient  6  . 

Geernaert,  Katsaros,  and  Richter  (1986)  note  that  typical  measured  values  of 
;  range  from  0.008  to  0.020,  and  found  values  spanning  this  range  from  wave 
measurements  in.  16-m  water  depth  in  the  North  Sea.  With  the  more  frequentlv 
cited  <  =  0.41  and  still  assuming  z  -  k  / 30  ,  this  yields  Charnock's 
coefficient  in  she  range  0.015  <  a  <  0.024  .  This  is  consistent  with  the 
range  cited  above  and  suggests  that  parameters  which  affect  8  such  as  fetch, 
water  depth,  or  degree  of  saturation  (Geernaert  et  al.  1986)  will  have  a 
direct  effect  on  a  .  The  same  effects  can  also  violate  the  assumptions 
leading  to  this  regime  as  is  discussed  below. 

47.  In  the  first  two  regimes,  Kitaigorodskii  (1973)  found  a  reasonable 

correlation  between  his  model  and  observations  in  cases  where  waves  were  only 

of  local  wind-driven  type.  Estimates  of  z  scattered  about  the  theoretical 

o 

curves  by  about  one  order  of  magnitude.  Compared  with  the  seven  orders  of 
magnitude  discussed  above,  this  is  a  major  improvement.  He  classified  his 
regimes  by  ratio  of  u^  to  spectral  peak  phase  speed  =  g/u^  such  that 
for  u./C  >>  1  (regime  of  initial  wave  growth)  and  u,/C  -  1  (intermediate 
wave  growth  regime)  his  hypotheses  were  justified. 

48.  However,  for  the  last  stages  of  wave  growth  (u./C  <<  1),  i.e.,  the 

regime  for  which  Charnock's  relation  was  derived,  Kitaigorodskii  did  not  find 
a  clear  correlation  between  his  model  and  his  observations  in  the  Central 
Pacific  Ocean.  Specifically,  for  u./C  <  0.05  ,  he  found  z  to  be  highlv 
variable,  generally  larger  than  the  exponentially  weighted  sea-surface  stan¬ 
dard  deviation,  and  inconsistent  with  a  constant  coefficient  in  Charnock's 
expression.  Part  of  the  reason  for  this  il-r  in  the  way  Kitaigorodskii ' s 
model  is  derived.  Equation  9  relates  wind  speed  and  wave-phase  speed.  Noting 

from  Table  l  that  drag  coefficients  at  standard  elevation  (z  =  10  m)  are 

22  -3  -3  r 

typically  in  the  range  C^  =  u^/U^  =  0.5  *  10  to  3  *  10  ,  it  is  seen  that 

u./U  =  0.02  to  0.05  .  Comparing  this  with  u./C  -  0.05  indicates  that  the 
*  r  *o 

point  of  inconsistency  in  Kitaigorodskii ' s  theory  occurs  when  C^  =  U  . 
Arguments  leading  to  Equation  9  essentially  treat  the  sea  surface  as  rigid  and 


moving  at  phase  speed  C  .  If  C  =  U  ,  then  there  would  be  not  mean  shear 
and,  hence,  no  stress.  Finite  stresses  were  reported  bv  Kitaigorodskii  0973'i 
under  these  conditions  so  that  additional  physical  mechanisms  apparently  be¬ 
come  dominant  when  seas  become  more  fully  developed.  His  results  suggest  that 
his  model  is  valid  for  initial  and  intermediate  stages  of  growth  as  identified 
bv  the  parameter  u./C  .  The  inverse  of  this,  C  /u.  ,  called  the  saturation 
parameter  by  Geernaert  and  Katsaros  (1986),  is  considered  significant  for 
rejection  of  Kitaigorodskii ’ s  model  when  it  is  on  the  order  of  20  or  larger. 

49.  In  the  discussion  of  the  work  of  Charnock  and  Kitaigorodskii,  the 

surface  roughness  parameter  z  has  been  discussed  in  terms  of  six  integral 

o 

parameters  associated  with  waves.  Thev  are:  u.  ,  g  ,  v  ,  C  ,  a  ,  and 

*  o  n 

.  Along  with  z^  this  forms  a  set  of  seven  variables  in  which  length  and 
time  are  the  only  dimensions.  From  dimensional  analysis  a  set  of  five  dimen¬ 
sionless  groups  can  be  formed  from  these  variables.  Combinations  formed  from 

2 

discussions  above  are  z  g/u.  (Charnock's  equation,  if  constant),  z  u./v 

o  *  o  * 

(dvnamicallv  smooth  flow,  if  constant),  C  /u.  (saturation  parameter),  z  / z 
(Kitaigorodskii ' s  initial  wave  growth  formula,  if  constant),  and  "r/’o  (wave 
steepness) .  This  is  a  complicated  set  with  different  parameters  or  combina¬ 
tions  of  parameters  important  in  different  regimes,  only  a  few  of  which  have 
been  discussed  above.  Other  combinations  of  parameters  have  been  proposed  by 
other  investigators. 

50.  For  example,  Hsu  (1972)  made  wind-stress  measurements  over  water  of 
depth  less  than  5  m  on  the  Florida  coast  of  the  Gulf  of  Mexico  and  found 
rather  large  drag  coefficients  (see  Table  1).  In  consideration  of  water  shal¬ 
lowness  and  its  effect  on  wave  steepness,  he  proposed  (Hsu  1974,  1976)  that 
Charnock's  relationship  be  modified  by  the  steepness  of  the  dominant  waves  in 
the  form 

,  2 
4o  u. 

z  <*  t -  (15) 
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where  4o^  is  significant  wave  height  assuming  a  Rayleigh  distribution  of 
wave  heights.  He  did  not  test  this  expression  in  shallow  water,  but  noted 

that  for  linear  deepwater  waves,  X  g  =  2  such  that  Equation  15  becomes 

o  o 
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in  deep  water.  Hsu  (  1974')  evaluated  this  expression  citing  19  sources,  but 
only  34  data  points,  all  from  deepwater  observations.  Data  were  grouped  in 
five  ranges  of  values  of  the  right  side  of  Equation  16  and  averaged  geometri¬ 
cally.  He  found  the  mean  values  to  correlate  very  well  with  Equation  16  using 
a  coefficient  of  1.0.  Maximum  deviation  by  mean  values  was  about  20  percent. 
Standard  deviations  of  the  geometric  averages  of  z^  were  about  half  an  order 
of  magnitude,  similar  tc  what  Kitaigorodskii  (1973)  reported  for  his  inter¬ 
mediate  stage  roughness. 

51.  In  a  second  test,  Hsu  (1976)  found  similar  results  from  four  addi¬ 
tional  sources.  It  is  of  interest  to  note  that  for  most  of  the  data  cited  in 

this  test  the  saturation  parameter  was  in  the  range  20  <  C  /u.  <  90  ,  i.e.,  a 

o  * 

region  where  Kitaigorodskii ' s  (1973)  model  was  found  to  be  inadequate.  A 

major  difference  between  the  two  models  is  that  Kitaigorodskii ' s  model  becomes 

independent  of  C  for  C  / u.  >>  1  ,  whereas,  Hsu's  model  does  not.  This 
o  O  * 

suggests  that  Hsu's  model  (or  at  least  a  continued  dependence  on  C  /u.)  mav 

o  * 

be  more  appropriate  in  the  latter  stages  of  sea  development. 

52.  SethuRaman  (1978)  reports  a  dependence  of  drag  coefficient  on 

C^/u^  from  measurements  in  water  depths  of  33  m  at  a  location  5  km  from  the 

south  shore  of  Long  Island.  For  C  (m  varying  from  about  15  to  about  50, 

o  * 

the  drag  coefficient  (reference  elevation  z  =  8  m)  dropped  from  about 

-3  -3  r 

2  *  10  to  about  0.5  x  10  .  In  accordance  with  Equation  7,  a  drop  in 

corresponds  to  a  reduction  in  z  .  SethuRaman  did  not  include  c  in  his 

o  n 

analysis,  so  Hsu's  model  was  not  tested.  There  is  sufficient  scatter  in  the 
data  to  suoport  the  hvpothesis  of  C,  =  constant  for  C  /u.  >  35  ,  i.e.,  C_ 
independent  of  C^/u^  •  Hence,  these  data  support  the  ideas  of  general  depen¬ 
dence  of  0  on  C  /u  ,  but  cannot  be  used  to  distinguish  validity  for 
either  Hsu's  or  Kitaigorodskii ' s  model. 

53.  A  theoretical  and  experimental  treatment  of  these  models  was  re¬ 
ported  recently  by  Huang  et  al.  (1986).  In  this  work,  they  used 
Kitaigorodskii ' s  definition  of  physical  roughness,  given  here  as  Equation  13, 
coupled  with  a  wave  spectral  deftnition  for  the  band  of  frequencies  for  which 
waves  are  fully  developed — known  as  the  simplified  Wallops  spectrum 
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(Huang  et  al.  198!).  This  spectrum  has  the  form 


where  m  is  a  function  in  the  simplified  Wallops  spectrum 

U^'2  -  vol2 

m  =  - t - = - — 

log  2 

8  =  (m  -  1)(2tto  /\  )2 
n  o 

and  is  the  spectral  peak  frequency  (in  radians  per  second) .  The  spectral 

shape  is  fixed  when  o)q  and  o^/X^  are  specified.  Again,  by  using  the  deep¬ 
water  dispersion  relation  in  the  exponent  of  Equation  13,  substituting  Equa¬ 
tion  17  in  Equation  13  and  integrating  yields  the  general  result 


where  G  ,  a  function  in  the  simplified  Wallops  spectrum,  is  defined  as 
follows : 


where  y  is  the  incomplete  gamma  function  (Abramowitz  and  Stegun  1970) . 

54.  Huang  et  al.  (1986)  note  that  although  Equation  18  is  a  somewhat 

complicated  function,  the  result  relates  two  models  through  a  third  under  the 

assumption  z  «  k  .  To  within  a  coefficient  (assumed  constant)  the  function 
o  R 

Go  /X  corresponds  to  Charnock's  (1955)  coefficient  a  in  Equation  8.  Also 
no 
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to  within  a  coefficient,  the  function  G  corresponds  to  a  generalized  coeffi¬ 
cient  for  Hsu's  ( 1972)  model,  given  here  as  Equation  15.  Both  are  found  from 
Kitaigorodskii ' s  (1973)  integral. 

55.  The  importance  of  the  work  of  Huang  et  al.  (1986)  is  that  their 

spectral  model,  derived  from  other  considerations,  suggests  one  possible  set 

o 

of  relationships  among  the  dimensionless  parameters  z  g/u7  and  a  !\  when 

0*^0 

incorporated  in  Kitaigorodski i ' s  integral.  It  shows  that  generalized  forms  of 
both  Charnock's  and  Hsu's  models  can  evolve  from  this.  On  the  other  hand,  it 
illustrates  the  sensitivity  of  Kitaigorodskii ' s  roughness  definition  to  the 
spectral  form  employed  (compare  Equation  18  with  the  result  obtained  above 
when  Phillips'  spectrum.  Equation  14,  is  used).  It  is  also  restricted,  as 
discussed  above  for  Kitaigorodskii ' s  results,  to  developing  seas,  i.e., 

Cq/u^  <  20  ,  and  to  deepwater  conditions. 

56.  The  experiment  program  of  Huang  et  al.  (1986)  had  some  curious  fea¬ 
tures.  They  required  a  coefficient  of  proportionality  between  z  and  k  , 

o  R 

the  latter  represented  here  as  Equation  18.  To  find  this,  they  used  data  from 

a  wind-wave  tank.  In  one  set  of  experiments,  the  relation  k  /z  was  found 

for  waves  generated  by  wind  alone.  In  six  other  sets,  paddle-generated  waves 

of  varying  amplitude  and  spectral  content  were  established  in  the  water  during 

wind-stress  observations.  Seven  mean  values,  one  for  each  set  of  runs,  of 

k^/z^  were  found  to  lie  in  the  range  11  <_  kj,/zo  <,  30  with  a  net  mean 

k„/z  =  20  .  This  is  close  to  k_/z  =  30  for  rough  flow  over  uniform  sand 
R  o  R  o 

grains  and  well  within  the  range  10  <  k_/z  <  200  described  by  Schlichting 
(1979)  for  less  regular  roughness.  It  is  also  close  to  the  value 
k  /z  =  17  ,  found  by  Kitaigorodskii  (1973). 

57.  However,  for  each  of  the  seven  sets  of  runs  the  standard  deviation 

was  large,  varying  from  41  to  123  percent  of  the  mean,  with  a  net  standard 

deviation  of  87  percent.  Furthermore,  there  was  a  rough  indication  of  a 

decrease  in  k  /z  for  increasing  background  wave  amplitude,  though  the  scat- 

ter  was  so  large  that  no  statistical  significance  could  be  assigned  to  this. 

Huang  et  al.  (1986)  were  careful  to  examine  data  having  C  / u.  <  10  which 

would  correspond  to  initial  wave  growth  regions  in  both  their  model  and  that 

of  Kitaigorodskii  (1973).  In  this  regime  z  should  vary  with  a  directlv. 

o  n 

The  large  scatter  suggests  either  that  the  model  is  overly  simple  or  that  sec¬ 
ondary  effects  in  the  experiment  regime  (horizontal  nonuniformity,  secondary 
flow  in  either  air  or  water,  unnatural  side  or  end  boundary  conditions) 
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modified  Cheir  results.  The  indicated  variation  with  background  wave  condi¬ 
tions  suggests  that  when  wind  and  waves  are  not  in  equilibrium  the  effective 
surface  roughness  is  modified. 

58.  The  data  of  Huang  et  al.  (1986),  in  combination  with  laboratory  and 

field  observations  from  eight  other  campaigns  and  when  plotted  in  the  form 

z  /c  versus  (C  /u.)  1  ,  show  a  distinct  decrease  with  increasing  C  /u.  as 
o  r  o  *  o  * 

was  found  bv  Kitaigorodskii  (1973)  and  Hsu  (1976).  The  scatter  in  z  /o  is 

o  n 

about  one  order  of  magnitude  as  has  been  noted  also  in  results  of  these  latter 

two  investigations.  The  model  results  of  Huang  et  al.  (1986)  pass  through 

these  data  and  so  are  viable.  It  is  noted,  however,  that  the  data  source  for 

the  largest  C  /u.  used  by  Huang  et  al.  (1986)  is  the  same  as  was  used  bv  Hsu 
o 

(1976)  so  that  within  the  variability  of  the  data  the  two  models  are 
indistinguishable . 


Model  Comparisons,  Observations,  and  Additional  Parameters 


59.  In  an  effort  to  compare  several  models  for  drag  coefficients, 
Geernaert,  Katsaros,  and  Richter  (1986)  conducted  wind  stress  and  wave  mea¬ 
surements  from  a  fixed  mast  in  water  of  16-m  depth  located  28  km  offshore  of 
West  Germany  in  the  North  Sea.  In  this  study  they  noted  that  tides,  waves, 
and  wind  all  induce  mean  currents  such  that  the  water  surface  mean  flow  bound¬ 
ary  condition  is  not  Uq  =  0  in  Equations  5  and  6  as  has  been  assumed  in  all 

investigations  discussed  above.  Thev  estimated  lT  from  tidal  theory,  local 

o 

tide  records,  and  estimates  of  sea-surface  slope.  The  drag  coefficient  was 
then  estimated  from  mean  wind  speed,  the  component  of  in  the  direction  of 

the  wind  and  measured  wind  stress  using  Equation  6,  first  equality.  Data  they 
used  were  constrained  to  speeds  above  6  ra/sec  from  a  direction  of  infinite 
fetch  and  were  also  corrected  for  atmospheric  stability  as  discussed  below. 

Six  models  of  were  then  tested  by  incorporation  in  Equation  6,  second 

equality,  and  comparison  with  measured  values  of  .  Comparison  was  done  in 

several  ways,  and  in  each  case  the  models  of  Charnock  (1955)  (Equation  8), 
Kitaigorodskii  (1973)  (Equation  13),  and  Hsu  (1974)  (Equation  15)  performed 
best  based  on  mean  square  deviation  from  measurements.  Measured  rather  than 
idealized  wave  spectra  were  incorporated  in  Kitaigorodskii ' s  model. 

60.  In  the  first  comparison,  the  single  free  coefficient  required  for 
each  model  was  determined  from  the  author's  own  work  or  from  published  best 


estimates  from  other  investigators.  For  example,  Charnock's  a  of  0.0144 
from  Garratt  (1977)  and  0.0185  from  Wu  (1980)  were  both  compared.  In  this 
test,  Charnock's  model  with  ol  =  0.0185  and  Hsu's  model  had  the  lowest  root 
mean  square  (RMS)  differences  of  about  10  percent  of  the  mean  measured  drag 

coefficient  given  by  Geernaert,  Katsaros,  and  Richter  (1986)  (this  was 

-3 

CQ  -  1.6  *  10  ).  Charnock's  model  with  a  =  0.0144  was  off  by  about  12  per¬ 
cent  and  Kitaigorodskii's  model  was  off  by  about  27  percent.  In  this  test  the 
measured  adjusted  for  the  estimated  drift  current  (but  not  the  tidal  cur¬ 

rent,  for  an  unexplained  reason)  to  be  nearer  the  zero-drift  surface  velocity 
assumed  by  the  modelers. 

61.  In  the  second  test,  the  coefficients  were  allowed  to  vary  to  mini¬ 
mize  differences  between  models  and  observations.  Again,  zero-drift  surface 
velocity  was  assumed.  Charnock's  coefficient  was  found  to  be  0.0192,  Hsu's 

coefficient  went  from  1.00  to  0.83,  and  Kitaigorodskii's  z  =  k  / 30  became 

o  R 

z^  =  k^/106  .  With  these  adjustments  all  three  models  had  about  10  percent 
RMS  differences  from  observations. 

62.  In  the  third  test,  surface  drift  velocity  was  included,  and  best 
fit  coefficients  were  computed  again.  In  this  test,  the  models  still  had 
about  10-percent  RMS  differences  from  observations,  but  the  coefficients 
changed.  Charnock's  a  became  0.0288,  Hsu's  coefficient  became  1,76,  and 
Kitaigorodskii's  zq  became  k^/71  • 

63.  Evidently,  the  presence  of  surface  currerts  has  an  effect  on  model 
results.  This  may  account  for  some  of  the  scatter  in  results  discussed  above 
where  surface  currents  were  not  considered.  This  also  suggests  that,  in 
regions  where  surface  currents  are  strong,  as  occurs  frequently  nearshore, 
such  effects  need  to  be  incorporated  in  both  observations  and  models. 

64.  Geernaert,  Katsaros,  and  Richter  (1986)  conclude  that 

Kitaigorodskii's  formulation  performed  best,  to  within  a  coefficient,  but  that 

results  are  too  close  to  make  a  strong  distinction  among  the  various  models. 

As  they  note,  their  data  are  from  one  site,  are  limited  in  number  (53  points 

were  used),  and  are  limited  in  range  of  wind  speed  (6  m/sec  <  U  <20  m/sec). 

r  *— 

Additional,  careful  field  observations  are  required  to  discern  the  best  of 
existing  models. 

65.  However,  Geernaert,  Katsaros,  and  Richter  (1986)  point  out  two 
additional  characteristics  of  their  observations  that  may  be  of  great  impor¬ 
tance  in  coastal  and  other  shallow-water  environments.  The  first  of  these  is 
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temporal  variability  of  the  drag  coefficient  and  the  second  is  an  apparent 
tendency  fcr  drag  coefficients  to  increase  in  shallow  water. 

66.  As  regards  the  first  effect,  Geernaert,  Katsaros,  and  Richter 
(19861  made  additional  measurements  after  their  primary  experiment  at  a  time 
when  their  wave-measuring  devices  had  been  removed.  Detailed  observations 
were  made  at  ti-  surrounding  the  passage  of  a  storm  front  wherein  wind 
changed  direction  by  roughly  60  deg  and  speed  dropped  from  about  20  m/s ec  to 
about  12  m/sec  in  the  course  of  about  2  hr.  Since  they  lacked  sea  state  in¬ 
formation,  they  predicted  a  drag  coefficient  based  on  a  regression  of  measured 
drag  coefficients  on  wind  speed  difference  (10-m  wind  speed  minus  water  sur¬ 
face  speed)  from  their  primary  experiment.  These  data  had  about  10-percent 
scatter  about  the  regression  line.  However,  when  they  compared  predicted  Cp 
during  the  frontal  passage,  they  found  variations  of  +20  percent,  which  is 
significant. 

67.  Other  observers  have  reported  similar  findings.  SethuRaman  (1978) 
reported  a  100-percent  change  in  drag  coefficient  in  the  course  of  2  hr  when 
wind  changed  direction  by  about  80  deg.  His  observations  included  some  wave 
information.  He  noted  that  when  the  sea  was  in  the  process  of  adjusting  to 

the  new  wind  direction,  the  saturation  parameter  (C  /u.)  rose  from  about 

o  * 

13  immediately  following  the  wind  shift  (when  doubled)  to  about  25  an 

hour  later  when  CQ  dropped  again.  This  behavior  is  qualitatively  consistent 
with  models  discussed  above  (all  of  them,  in  fact,  excepting  Charnock's  model) 
wherein  C  depends  strongly  on  C  /u.  . 

68.  A  third  observation  is  described  by  Boyle,  Davidson,  and  Spiel 
(1987)  based  on  an  experiment  in  deep  water  in  the  Gulf  of  Alaska.  They  also 
estimated  CD  from  a  purely  wind-speed-dependent  model  (Large  and  Pond  1982) , 
corrected  for  atmospheric  stability,  and  used  this  to  estimate  wind  stress. 
Measured  wind  stresses  were  then  compared  with  estimates.  They  found  very 
good  correlation  between  the  two  during  steady  wind  conditions.  However, 
immediately  preceding  passages  of  large-scale  weather  systems  (fronts,  pres¬ 
sure  troughs,  and  pressure  ridges)  they  observed  that  measured  stresses  ex¬ 
ceeded  predicted  stresses  by  a  factor  of  two  to  three.  They  conclude  that 
adjustments  of  sea  state  (of  which  they  had  no  quantitative  observations)  are 
extremely  important  in  unsteady  conditions. 

69.  The  importance  of  these  observations  to  modelers  of  waves  and  cur¬ 
rents  is  twofold.  It  is  common  practice  in  modeling  to  employ  a  wind  drag 
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coefficient  based  on  wind  speed  alone  (see,  e.g.,  Sheng  and  Butler  (1982)  cr 
Hubertz  (1987)).  Many  investigations  have  proposed  formulae  for  such  depen¬ 
dence.  A  few  of  these  have  been  referenced  in  this  report:  Hsu  (1972), 

Garratt  (1977),  Ku  (1980),  Amorocho  and  DeVries  (1980),  and  Large  and  Pond 
(1982).  However,  if  wind  stress  depends  on  other  parameters,  such  as  those 
associated  with  sea  state,  then  the  more  elementary  models  will  be  incorrect 
at  least  part  of  the  time.  The  times  at  which  such  models  are  likely  to  be 
remiss  are  those  of  most  importance  in  coastal  dynamics,  i.e.,  during  storms, 
where  observations  noted  above  indicate  large  variations  in  .  Clearly, 

further  research  is  needed  to  relate  to  sea  state  and  wind  speed,  either 

to  verify  existing  models  or  to  identify  parametric  dependencies  not  yet 
proposed . 

70.  The  other  important  consequence  of  temporal  variation  is  a  little 

more  heuristic  but,  nonetheless,  plausible.  The  changes  in  observed  dur¬ 
ing  rapid  changes  in  wind  conditions  suggest  that  will  vary  whenever  wind 

and  waves  are  not  in  equilibrium.  The  models  described  above  were  almost 
entirely  in  consideration  of  unobstructed,  steady  winds  acting  in  the  initial 
and  intermediate  stages  of  wave  growth  in  deep  water  and  with  wind  primarily 
in  the  direction  of  wave  propagation.  If  the  wind  changes  direction  alone 
then  there  will  be  a  tendency  for  a  modification  of  the  existing  sea  and  ini¬ 
tiation  of  a  sea  in  the  new  wind  direction. 

71.  A  simplified  application  of  this  idea  can  be  used  in  conjunction 
with  Kitaigorodskii ' s  (1973)  model  as  represented  by  Equation  9.  If  the  wind 
speed  is,  say,  twice  the  phase  speed  such  that  U  =  2C  and  U  -  C  =  C  ,  then 
a  simple  180-deg  reversal  of  wind  causes  U  -  C  *  -2C  -  C  =  -3C  since  the 
wave  phase  speed  vector  is  now  opposite  the  wind  vector.  The  ratio  of  magni¬ 
tudes  of  relative  winds  is  a  factor  of  3.  The  square  of  this  is  roughly  the 
ratio  of  stresses,  which  is  9,  with  the  larger  stress  acting  until  the  sea 
adjusts  to  the  new  wind  direction.  While  this  is  an  oversimplification  of  the 
problem,  the  qualitative  concept  is  clear. 

72.  A  similar  effect  would  be  expected  for  wind  waves  entering  or  being 
generated  in  shoaling  waters  such  as  coastal  regions.  Water  of  finite  and 
spatially  varying  depth  is  known  to  modify  waves  independently  of  the  wind  and 
so  can  change  the  equilibrium  conditions  described  above  for  wind  acting  in 
deep  water.  First,  wavelengths  decrease  so  that  even  if  wave  heights  are  con¬ 
stant,  the  wave  steepness  increases.  This  means  that  if  Hsu's  (1974)  model, 
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Equation  15,  is  correct,  then  surface  roughness  and  therefore  the  drag  coeffi¬ 
cient  increases.  Second,  phase  speeds  (at  least  of  some  waves')  are  modified 

by  finite  depth.  In  the  limit  of  shallow  water,  all  phase  speeds  are 

1/2 

C  =  (gdt  where  d  is  water  depth.  In  Kitaigorodskii ' s  (1973)  model. 
Equation  12,  C  is  then  no  longer  a  function  of  cu  so  the  weighting  function 
can  be  removed  from  Equation  13.  The  solution  of  Equation  13  is  then 

1  /° 

-<(gd)  7“/u* 

=  2an  e  (19) 

Since  shallow-water  phase  speeds  are  always  less  than  deepwater  phase  speeds 

at  the  same  frequency,  then  k  (and  C  )  will  tend  to  increase  in  shallow 

K  U 

water  (for  the  same  a  )  bv  this  model. 

n 

73.  A  third  effect  occurs  if  the  wind  is  not  parallel  to  the  gradient 

of  the  bottom.  In  this  case,  at  least  some  waves  will  propagate  by  refraction 
in  a  direction  different  from  the  wind.  It  is  difficult  to  make  a  hypothesis 
about  how  this  affects  wind  stress.  In  the  case  of  nearly  shore-parallel 
winds,  the  dominant  wave  crests  can  be  nearly  parallel  to  the  wind  such  that 
the  component  of  dominant  wave-phase  speed  vector  in  the  direction  of  wind  is 
negligible.  This  suggests  that  may  drop  out  of  the  problem  such  that 

both  directional  and  spectral  characteristics  of  high-frequency  waves  are 
important.  There  are  no  direct,  quantitative  observations  of  these  effects. 
This  lack  of  knowledge  emphasizes  the  need  for  a  field  program  specifically 
designed  to  examine  wind  stress  in  shallow  water  under  a  variety  of  measured 
wind  and  sea  conditions. 

74.  A  fourth  effect  is  that  wind-  and  wave-induced  currents,  which  can 
reach  to  a  depth  on  the  order  of  30  m  in  deep  water  (Long  1981) ,  as  well  as 
tidal  and  other  currents  all  extend  through  the  whole  water  column  in  shallow 
water  and  so  are  affected  by  the  bottom.  Bottom  friction  tends  to  retard  the 
flow.  Wave  radiation  stresses  tend  to  drive  alongshore  flows  of  cross-shore 
varying  intensity.  All  of  these  effects  are  different  from  deepwater  behavior 
and  so  will  affect  Uq  in  Equations  5  and  6  with  a  consonant  effect  on  , 
Note  that  the  retarding  effect  of  bottom  friction  on  purely  wind-driven  mean 
flow  will  reduce  Uq  in  shallow  water  and,  by  Equation  6,  increase  . 

75.  A  fifth  effect  is  a  constraint  on  o  in  shallow  water.  There  is 

n 

a  tendency  for  waves  to  break  when  wave  height  and  depth  are  of  the  same 
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order.  Surface  roughness  estimates  based  on  sea  surface  variance.  Equation  19 
for  example,  would  then  be  reduced.  Hence,  would  be  expected  to  be 

smaller  in  the  breaker  zone  than  in  deep  water  due  to  is  effect. 

76.  There  is  little  direct  evidence  of  a  quantitative  relationship  be¬ 

tween  Cp  and  water  depth  though  there  is  some  indirect  evidence  that  water 
depth  is  important.  The  large  found  by  Hsu  (1972)  (see  Table  1)  in  verv 

shallow  water  may  be  due  to  shoaling  effects.  Storm  surge  models  sometimes 
require  larger  wind  stress  coefficients  to  correspond  with  observations.* 
Geernaert,  Katsaros,  and  Richter  (1986)  noted  that  their  regression  of  on 

wind  speed  had  larger  values  at  high  wind  speeds  than  curves  derived  from  ob¬ 
servations  by  Garratt  (1977)  and  Large  and  Pond  (1981).  These  latter  two 
investigations  were  in  the  open  ocean  whereas  Geernaert,  Katsaros,  and  Richter 
(1986)  made  their  observations  in  16-m  water  depth.  Their  observations  com¬ 
pared  favorably  with  several  other  investigations  in  shallow  water  of  similar 
depths.  They  conclude  that  C^/u^  is  an  important  parameter  in  coastal  re¬ 
gions  because  depends  implicitly  on  depth  and  also  on  fetch  which  is 

often  considered  an  important  parameter  of  wind  stress  (Wu  1980) . 

77.  While  there  is  some  merit  in  the  above  conclusion,  the  complete 
coastal  problem  may  be  more  complex  in  light  of  possible  misalignment  of  winds 
and  dominant  waves  due  to  shoaling,  as  mentioned  above.  In  this  case,  wind 
stress  may  depend  not  only  on  the  variance  in  the  wave  frequency  band  most 
strongly  coupled  to  the  wind  but  also  on  the  way  in  which  these  waves  are 
modified  by  (if  they  are  not  identical  to)  other  waves,  particularly  low- 
frequency  waves,  in  the  spectrum.  If  low-frequency  waves  are  important  to 
wind  stresses,  then  their  effect  will  be  more  dramatic  in  shoaling  water  be¬ 
cause  they  are  affected  by  the  sea  bottom  in  deeper  water  than  are  higher 
frequency  waves. 

78.  There  are  two  scaling  arguments  which  suggest  that  dominant  wind 
input  to  waves  occurs  at  frequencies  equal  to  or  higher  than  the  spectral  peak 
frequency.  The  first  of  these  is  based  on  Kitaigorodskii ' s  integral,  Equa¬ 
tion  13,  using  Phillips'  spectrum,  Equation  14.  As  discussed  by  Geernaert, 
Katsaros,  and  Richter  (1986),  the  frequency  of  maximum  contribution  to  surface 
roughness  can  be  found  by  differentiating  the  integral  of  Equation  13  with 

*  Personal  Communication.  H.  Lee  Butler,  1987,  Chief,  Research  Division, 
Coastal  Engineering  Research  Center,  US  Array  Engineer  Waterways  Experiment 
Station,  Vicksburg,  MS. 
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respect  to  .  setting  the  result  equal  to  ;ro,  and  solving  for  u.  , 

result  here  being  called  id  .  Using  Equation  14,  this  vields  w 

max  -  max 

=  0.4icg/u^  .  In  deep  water,  the  spectral  peak  phase  speed  C  is  related  to 

o 

peak  frequency  w  by  C  =  g/w  .  From  this,  <d  =  g/C  .  The  ratio  of 
0  0  0  0  0 

these  frequencies  is  w  /id  =  0.4<C  /u.  .  Since  Equation  14  allows  no  vari- 

max  o  o  * 

ance  for  w  <  d  ,  then  <d  =  w  for  0.4<C  /u.  <  1  or  C  /ux  <  6  (for 
o  max  o  o  *  o  *  ~ 

<  -  0.4).  This  corresponds  roughly  to  the  initial  stages  of  wave  growth  as 

identified  by  Kitaigorodskii  (1973)  (wherein  z  a  o  )  and  the  region  identi- 

o  n 

fied  by  Huang  et  al.  (1986)  using  the  simplified  Wallops  spectrum  (1  <  C  /u* 

<  10),  where,  likewise,  z 0  "  ari  •  In  the  latter  stages  of  wave  growth,  where 

C  / u.  -  20  ,  the  ratio  of  frequencies  becomes  id  /id  =  3  .  This  result  sug- 
o  *  max  o 

gests  that  for  much  of  wind  wave  growth  in  deep  water  <d  >  u>  and  the  pri- 

max  o 

mary  contribution  f  om  wind  occurs  at  frequencies  higher  than  the  spectral 
peak. 

79.  The  second  argument  arises  from  the  model  of  Charnock  (1955)  and  is 
discussed  by  Wu  (1980).  If  zq  =  ^u^/g  ,  with  a  in  the  range  0.01  £  n 

-  0.03  ,  and  if  z  11  k  (i.e.,  dynamically  rough  flow)  with  the  ratio  z  /k_ 
broadly  in  the  range  1/10  £  z  /k  ^  1/200  ,  as  discussed  previously,  then  k 

2  0  “  2  ° 

will  be  in  the  range  0.1  u^/g  <  <_  6  u^/g  *  Furthermore,  if  kR  is 

representative  of  the  heights  of  waves  at  saturation,  i.e.,  at  a  limiting 

steepness  of  k  /X  -  1/7  (Kinsman  1984),  then  the  wavelength  X  is  in  the 
2  K  2 

range  0.7  u^/g  <,  X  42  u*/g  .  In  deep  water,  X  can  be  expressed  in  terms 

2 

of  phase  speed  C  as  X  =  2^C  /g  ,  If  this  is  inserted  in  the  previous 

expression,  the  ratio  C  /u.  ,  where  C  is  the  C  of  maximum  interac- 

max  *  max 

tion,  is  found  to  be  in  the  range  0.3  <  C  /u.  <2.6  .  That  is,  under  the 

max  * 

assumptions  given,  the  waves  which  contribute  to  Charnock's  roughness  scale 

all  have  phase  speeds  C  of  order  u.  .  Since  C  /u.  based  on  spectral- 

peak  phase  speed  can  be  equal  to  or  larger  than  this  (values  to  at  least 

C  /u.  -  20),  then  C  /C  >  1  .  From  this,  (g/w  )/(g/u>  )  =  id  /<d  >,  1 

o  *  o  max  o  max  max  o 

such  that  contributing  frequencies  can  be  equal  to  or  up  to  20  times  the  spec¬ 
tral  peak  frequency. 

80.  This  last  scaling  argument,  with  the  result  that  C  ~  u*  »  has  a 
consequence  which  could  possibly  introduce  another  parameter  into  the  problem. 
According  to  linear  wave  theory,  if  surface  tension  is  included  with  gravity 
as  a  restoring  force,  then  there  occurs  a  minimum  in  wave  phase  speed  which  is 
on  the  order  of  0.2  m/sec  (Kinsman  1984).  This  corresponds  to  a  wavelength  of 
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about  1.7  cm  for  tvpirqi  air-water  intarfacial  surface  tension.  Phase  speeds 
increase  rapidly  for  shorter  wavelengths  (pure  capillary  waves)  and  increase 
more  gradually  for  longer  wavelengths  (pure  gravity  waves) .  For  C  /u. 

-  1  ,  then  uA  must  be  about  0.2  m/sec  before  momentum  exchange  occurs  with 
roughness  characterized  by  Charnock's  model  (subject,  of  course,  to  limita¬ 
tions  in  the  above  scaling  arguments).  At  lower  u*  surface  roughness  may  be 
characterized  bv  the  dynamically  smooth  regime  (z  “  v/u.)  or  possibly  by  a 
combination  of  this  with  scales  associated  with  pure  capillary  waves. 

81.  A  value  of  u*  =  0.2  m/sec  corresponds  to  10-m  wind  speeds  of 

5.6  m/sec  <  <  6.2  m/sec  using  Charnock's  zq  with  a  in  the  range 

0.01  <  a  <  0.03  .  For  wind  speeds  lower  than  this,  one  might  expect  a  varia¬ 
tion  in  characterization  in  surface  roughness  which  would  have  an  effect  on 
drag  coefficient  estimates.  At  higher  wind  speeds,  which  are  generally  the 
wind  speeds  of  interest,  u*  increases.  If  C  =  u^  ,  then  wave  phase  speeds 
of  interest  are  more  nearly  those  of  pure  gravity  waves. 

82.  Munk  (1955),  Phillips  (1969),  Wu  (1980),  Kinsman  (1984),  and  others 

have  discussed  the  importance  of  capillary  waves  (in  which  surface  tension  is 
important)  in  wave  growth.  Munk  (1955)  and  Geernaert,  Katsaros,  and  Richter 
(1986)  cite  experiments  in  which  a  reduction  in  surface  tension  (by  addition 
of  detergent  or  by  biological  activity)  serves  to  reduce  in  general.  By 

the  above  scaling  a  reduction  in  surface  tension  reduces  the  minimum  phase 
speed  which  corresponds  to  a  smaller  u^  (for  u^  11  C) .  This  leads  to  a 
smaller  z^  by  Charnock's  model  and,  by  Equation  6,  to  a  smaller  .  This 

should  affect  primarily  C  for  winds  with  u.  near  this  minimum  phase 
speed.  If  it  affects  for  a  wider  range  of  wind  speeds,  then  this  sug¬ 

gests  that  the  presence  of  capillary-gravity  waves  is  more  generally  important 
in  air-sea  interaction.  As  Wu  (1980)  points  out,  this  has  not  been  fully  re¬ 
solved  yet.  It  is  less  important  if  surface  tension  (and  therefore  the  mini¬ 
mum  wave  phase  speed)  is  constant.  However,  in  coastal  waters  subject  to  sea¬ 
sonal  planktonic  blooms  and  riverborne  surfactants  (biological  or  otherwise) 
surface  tension  may  vary  significantly.  It  should  then  be  considered  a  poten¬ 
tially  important  parameter  in  nearshore  drag  coefficient  estimation. 

83.  Whether  of  capillary  nature  or  not,  the  waves  interacting  most 
strongly  with  wind  are  generally  at  higher  frequencies  than  the  spectral  peak 
frequency  by  the  above  scaling  arguments.  The  models  of  Kitaigorodskii  (1973) 
and  Hsu  (1974),  both  justified  with  experimental  evidence,  predict  a  decrease 
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Since  a 


is  the  total  variance 


in  the  ratio  z  /?  for  increasing  C  /u.  . 

of  sea-surface  displacement,  this  result  simply  means  that  low-f requencv 
waves,  especially  those  moving  at  near  the  mean  wind  speed,  contribute  less 
effectively  to  total  surface  roughness.  This  argument  applies  in  deep  water 
where  wind  and  waves  of  all  frequencies  are  moving  in  roughly  the  same  direc¬ 
tion.  In  shallow  water,  where  low-frequency  wave  phase  speeds  are  retarded 
and  their  directions  are  altered  by  refraction,  it  is  possible  that  more  or 
less  of  the  total  sea-surface  variance  will  contribute  to  the  effective  rough¬ 
ness.  If  long  waves  are  slowed  and  are  in  the  wind  direction,  it  is  expected 
that  and  therefore  will  increase  relative  to  the  deepwater  case.  If 

long  waves  have  refracted  so  their  direction  of  propagation  is  normal  to  the 
wind  direction,  it  is  possible  that  they  make  no  contribution  to  surface 
roughness  (if  long  crested,  so  no  bluff  surface  is  presented  to  the  wind)  or 
that  they  make  maximum  contribution  to  surface  roughness  (if  short  crested,  so 
that  larger  surface  displacements  are  effectively  unmoving  in  the  wind 
direction) . 

84.  Clearly,  the  net  effect  depends  critically  on  the  wind  direction 

and  the  detailed  frequency  distribution  of  wave  propagation  directions  and 
variances,  i.e.,  the  directional  wave  spectrum.  The  arguments  leading  to  this 
dependence  are  based  on  extensions  of  arguments  applied  to  deep  water  and  are 
therefore,  hypothetical.  Verification  of  these  hypotheses  requires  observa¬ 
tions.  To  date,  no  experiment  testing  the  dependence  of  C.^  on  the  direc¬ 
tional  spectrum  in  shallow  or  coastal  waters  h?s  bee"  to  be  reported  in 

the  literature.  This  is  one  of  the  primary  reasons  for  conducting  the  re¬ 
search  program  described  in  this  report. 

85.  In  summary,  this  section  describes  some  of  the  previous  work  lead¬ 

ing  to  formulations  of  the  dependence  of  the  drag  coefficient  on  various  ob¬ 
servable  parameters  related  to  dynamics  of  both  air  and  water  for  steady,  uni¬ 
form  conditions  in  the  deep  ocean.  All  the  models  discussed  showed  to 

depend  solely  on  a  characterization  of  the  surface  roughness  parameter  z^  . 

The  models  suggest  z  depends  variously  on  the  parameters  v  ,  g  ,  o  , 

O  q 

C  ,  u  ,  and,  possibly,  surface  tension.  Observed  variations  of  these 
o  * 

models  for  shallow  water,  changes  in  sea  state,  and  variation  of  surface  cur¬ 
rents  suggest  that  additional  parameters  are  necessary.  Extensions  of  argu¬ 
ments  leading  to  and  derived  from  the  deepwater  models  indicate  the  nature  of 
the  added  parameters.  They  are  direction  of  reference  wind  6^  ,  water 
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surface  veloritv  LT  and  its  direction  6  ,  water  depth  d  ,  and  replacement 

o  o 

of  and  C  with  variances  and  phase  speeds  in  frequency  and  directional 

bands,  i.e.,  ^ ,  9 )  and  C(cj,6)  ,  respectively.  Dependence  of  CQ  on 

these  added  parameters  is  to  be  determined  empirically  through  detailed  obser¬ 
vations.  Simplified  formulae  for  purposes  of  application  can  then  be  derived 
from  the  empirical  correlations. 
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PART  III:  PROCESSES  WHICH  MODIFY  THE  LOGARITHMIC  PROFILE 


86.  The  models  described  in  Part  II  all  depend  critically  on  the  exis¬ 
tence  of  the  logarithmic  wind  profile,  Equation  5,  as  derived  from  the  dimen¬ 
sional  arguments  leading  to  Equation  4.  For  Equation  5  to  be  valid,  the  flow 
and  Koundary  conditions  must  be  horizontally  uniform  and  steady.  Further,  the 
air  density  must  be  constant.  Violations  of  these  conditions  impart  added 
variables  to  the  flow  dynamics  and  can  modify  the  wind  profile  shape.  In  this 
case  different  wind  stresses  can  occur  for  the  same  reference  wind  speed  so 
that  errors  occur  when  a  drag  coefficient  is  used  with  wind  speed  to  estimate 
stress.  It  is  important,  therefore,  to  account  for  processes  which  modify  the 
wind  profile.  Improved  accuracy  is  obtained  where  this  can  be  done  quantita¬ 
tively.  Where  this  cannot  be  done,  qualitative  uncertainty  or  limits  of  ap¬ 
plicability  of  the  model  can  be  assigned. 

87.  The  three  general  conditions  listed  above,  density  variations,  un¬ 
steadiness,  and  nonuniformity,  are  discussed  in  this  section.  Reasonabl;-  pood 
estimates  can  be  made  of  the  effects  of  density  variation,  and  these  are  de¬ 
scribed  first.  Qualitative  arguments  on  the  effects  of  the  remaining  two  con¬ 
ditions  are  then  presented. 


Stratification 


88.  Air  density  can  vary  in  both  space  and  time.  For  steady,  uniform 
flow  the  primary  effect  of  density  variations  in  the  surface  layer  of  the 
atmospheric  boundary  layer  is  due  to  changes  of  density  in  the  vertical  direc¬ 
tion.  If  density  changes  in  the  vertical  direction,  two  effects  occur. 

First,  turbulence  processes  are  altered  due  to  buoyancy  effects  and,  second,  a 
flux  of  mass  tends  to  occur.  Buoyancy  inhibits  vertical  turbulent  motions 
under  stable  conditions  (density  decreasing  upward)  since  fluid  elements  dis¬ 
placed  upward  tend  to  be  more  dense  than  the  surrounding  fluid  such  that  grav¬ 
ity  acts  to  retard  their  upward  motion.  Fluid  elements  displaced  downward 
tend  to  be  less  dense  than  the  surrounding  fluid  such  that  buoyancy  retards 
their  downward  movement.  Under  unstable  conditions  (density  increasing  up¬ 
ward)  buoyancy  forces  tend  to  enhance  vertical  motions.  Mass  flux  occurs  be¬ 
cause  there  is  a  gradient  in  density  which  tends  to  cause  diffusion.  The  mass 
fluxes  of  interest  here  are  turbulent  fluxes.  Purely  molecular  fluxes  are 


36 


generally  much  smaller.  Relationships  are  thus  sought  which  show  the  effect 
of  stratification  on  the  mean  velocity  profile  and  the  fluxes  in  terms  of  mea¬ 
surable  mean  <=tate  variables. 

8Q.  One  of  the  more  successful  theories  leading  to  such  relationships 
is  the  Yonin-Obukhov  similarity  theory  (Tennek.es  and  I.umley  1972).  Guidance 
for  modification  of  the  mean  velocity  profile  by  stratification  is  provided  bv 
the  turbulent  kinetic  energy  equation,  derived  from  the  total  fluid  kinetic 
energy  equation  by  Reynolds  averaging.  For  steady,  homogeneous  flow  in  the 
x-direction,  this  equation  is 


u'w'  -  -  o'w'  -D 

ez  — 


where 

£  =  viscous  dissipation  rate  of  turbulent  kinetic  energy 
; '  =  fluctuating  density 

D  =  turbulent  kinetic  energy  gradient  diffusion  term 
The  left  side  of  Equation  20  represents  production  of  turbulent  kinetic 
energy.  The  first  term  on  the  right  of  Equation  20  is  a  symbol  that  stands 
for  viscous  dissipation  of  turbulent  kinetic  energy.  Its  formal  expression  is 


3  3  --  ,  TT 

V  V 

u  3x .  !x. 

i-l  3=1  J  J 


where 

=  u'  ,  v'  ,  or  w'  ,  respectively,  for  i  =  1  ,  2  ,  or  3 
x^  =  x  ,  v  ,  or  z  ,  respectively,  for  j  =  1  ,  2  ,  or  3 
It  represents  the  conversion  of  turbulence  to  heat  by  the  small-scale,  high- 
frequency  mean-square  strain  rate  at  which  molecular  viscosity  is  effective. 
The  second  term  on  the  right  of  Equation  20  embodies  the  direct  effect  of 
buoyancy  in  the  turbulent  energy  balance.  It  acts  as  a  source  or  a  sink  term 
depending  on  the  sign  of  p ' w'  .  The  last  term  on  the  right  of  Equation  20  is 
called  the  diffusion  term  and  is  given  by 
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where 

p'  =  fluctuating  pressure 
2  2  ?  2 

e'“  =  l/2(u'  +  v  +  w  )  ,  the  instantaneous  turbulent  kinetic  energv 

"ci  unit  mass 

It  is  called  the  diffusion  term  since  it  contains  gradients  of  turbulent  prop¬ 
erties.  It  is  very  difficult  to  quantify  because  fluctuating  pressure  and 
triple  velocity  correlations  are  difficult  to  measure.  It  is  generally  con¬ 
sidered  negligible  in  unstratified  flow  and  more  dominant  in  stratified  flow 
(Long  1981,  Frenzen  1983). 

90.  Equation  20  can  be  converted  to  usable  form  by  rewriting  its  com¬ 
ponents  in  terms  of  measurable  quantities.  For  the  buoyancy  term,  the  den¬ 
sities  can  be  converted  to  measurable  quantities  through  an  equation  of  state 

for  air.  This  is 


P 


pR  T 
a  v 


(23) 


where 

p  =  total  atmospheric  pressure 

R  =  modified  ideal  gas  constant 
a 

T  =  virtual  temperature 

v 

Virtual  temperature  is  an  atmospheric  state  variable  which  includes  the  ef¬ 
fects  of  both  sensible  (thermometric)  temperature  and  humidity.  In  terms  of 
measured  temperature  and  humidity,  the  virtual  temperature  is 

T  -  T(1  +  0.61q)  (24) 

v 


where 

T  =  air  temperature,  °K 

q  =  specific  humidity,  kg  of  water/kg  of  air 
A  summary  of  temperature,  humidity,  and  pressure  relationships  is  given  bv 
Blanc  (1985).  Expanding  Equations  23  and  24  in  terms  of  mean  and  fluctuating 
parts,  averaging,  subtracting  the  averaged  equations  from  the  expanded 
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equations,  dividing  by  the  averaged  equations,  and  retaining  the  dominant 
terms  yields,  to  within  about  1  percent, 

o’ 


w’here 

T’,  -  T '  +  0.61Tq  ' 

=  fluctuating  virtual  temperature  to  within  about  1  percent 
T'  =  fluctuating  (sensible)  air  temperature 
T  =  mean  (sensible)  air  temperature 
q'  =  fluctuating  specific  humidity 
T  =  T( 1  +  0.61q) 

=  mean  virtual  temperature  to  within  about  1  percent 
q  =  mean  specific  humidity 

91.  The  buoyancy  term  in  Equation  20  can  be  written  using  Equation  25 
to  yield 


-  £  o V  -  w'T'  (26a) 

-  *  v 

P  T 

v 

=  (w'T'  +  0.61  T  w ' q '")  (26b) 

T 

v 

Equation  26a  shows  the  buoyancy  term  as  a  function  of  the  flux  of  virtual  tem¬ 
perature,  Equation  26b  shows  it  in  terms  of  the  actual  fluxes  of  temperature 
and  water  vapor. 

92,  In  Monin-Obukhov  theory  (Tennekes  and  Lumley  1972),  turbulent 

fluxes  of  mass  and  momentum  are  all  considered  to  vary  weakly  with  z  from 

the  top  of  the  sublayer  through  the  surface  layer  in  the  same  way  that  the 

model  for  neutral  flow  was  derived.  Thus,  u^  ■  -  u * w'  as  before.  When  this 

is  used  with  Equation  26a  in  Equation  20  and  the  result  is  normalized  with  the 
3 

expression  kz/u^  ,  then  a  result  similar  to  Equation  4  is  obtained.  This  has 
the  form 
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(27) 


where  the  terms  are  in  the  same  order  as  in  Equation  20.  The  term  on  the  left 
of  Equation  27  is  a  dimensionless  shear.  In  unstratified  flow  there  is  no 
mass  flux  and  the  diffusion  term  is  considered  negligible.  In  this  case  the 
term  on  the  left  is  equal  to  unity  by  Equation  4.  The  only  remaining  term  in 
Equation  27  is  the  dissipation  term  (first  term  on  the  right)  which  must  also 
be  unity  by  identity.  Thus, 


3 

-  — 

KZ 


(28) 


is  a  reasonable  estimate  of  dissipation  in  neutral  flow. 

93.  When  the  flow  is  stratified,  the  second  term  on  the  right  of  Equa¬ 

tion  27  becomes  important.  It  can  be  written  as  -z/L  where 
L  =  -u^  T^/xg  w'T^  is  caHed  the  Monin-Obukhov  length  scale.  Stratification 
is  important  when  L  is  of  the  same  order  as  z  in  the  surface  layer.  Note 

that  L  *  00  (z/L  -*•  0)  when  there  is  no  mass  flux,  i.e.,  in  neutral  condi¬ 

tions.  The  effect  of  stratification  on  dissipation  and  diffusion  is  not  gen¬ 
erally  understood  so  Equation  27  cannot  be  used  directly  to  resolve  the  effect 
on  3l'/3z  .  However,  it  suggests  that  the  dimensionless  shear  is  a  function 
of  z/L  in  the  form 


kz  3D  /z\ 
u*  3z  "  %\jJ 


(29) 


with  q  (0)  =  1  to  correspond  to  Equation  4.  The  form  of  $  must  be  deter- 
m  m 

mined  empirically.  Note  that  the  assumption  of  dependence  of  z/L  above  is 
the  simplest  assumption  that  can  be  made  based  on  Equation  27.  If  verified 
experimentally,  the  assumption  can  be  said  to  be  justified. 

94.  In  the  same  sense  that  the  turbulent  kinetic  energy  equation  is  an 
equation  for  the  rate  of  change  of  variance  of  the  total  velocity  field,  equa¬ 
tions  for  the  variances  of  temperature  and  humidity  can  be  derived  by  Reynolds 
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averaging  the  p.ass  conservation  equations.  Relationships  between  fluxes  and 
gradients  of  temperature  and  humidity  similar  to  that  for  momentum  (Equa¬ 
tion  27^  are  then  found.  The  assumption  of  z/L  dependence  of  the  fluxes 
then  requires  formulae  for  temperature  and  humidity  similar  to  Equation  29. 
They  take  the  forms 


w'T' 


5T 

3Z 


(30) 


and 


u..  z 


w  q 


is.  = 

3z 


,(!) 


(315 


wh  ere 

$,j,  =  dimensionless  mean  temperature  gradient 

6^  =  dimensionless  mean  humidity  gradient 

95.  Empirical  expressions  for  <J>  ,  <t>-  ,  and  4>  have  been  proposed 

ml  q 

by  numerous  investigators  from  observations  over  land.  There  is  no  reason  to 
believe  that  the  fundamental  physics  is  different  over  water,  so  reliable 
land-based  models  should  be  applicable  at  sea.  One  set  of  expressions  which 
has  found  broad  support  in  independent  experiments  and  which  is  used  widely 
has  the  general  form  for  stable  conditions  (z/L  >  0) 


=  1  + 


z 

31  L 


-  *2  +  *3  l 


(32) 


where  a^  ,  a?  ,  and  a^  are  dimensionless  similarity  model  constants  and 
for  unstable  conditions  (z/L  <  0) 
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(33) 


where  b^  ,  b,,  ,  and  b^  are  dimensionless  similarity  model  constants. 

There  is  agreement  among  several  investigators  on  the  functional  forms  given 
by  Equations  32  and  33  but  some  disagreement  on  the  values  of  the  constants 
3^  >  >  a3  >  bj  «  b?  ,  and  b^  .  Table  2  lists  constants  recommended  by 

individual  investigators.  Note  that  von  Karman's  number  <  is  important  as 
well  by  virtue  of  Equation  29. 

Table  2 

Constants  Proposed  bv  Various  Investigators  for  the 
Similarity  Functions  Given  by  Equations  32  and  33 


Stable 

Unstable 

Source 

< 

al 

a2 

a3 

bl  b2 

b3 

Dyer  and  Hicks  (1970)* 

0.41 

_ 

_ 

_  _ 

16  1.00 

16 

Businger  et  al.  (1971)* 

0.35 

4.7 

0.74 

4.7 

15  0.74 

9 

Dyer  (1974) 

— 

5 

1.00 

5 

16  1.00 

16 

Lo  and  McBean  (1978)*,** 

0.40 

— 

— 

— 

-  - 

— 

Wierenga  (1980)* 

0.41 

6.9 

1.00 

9.2 

22  1.00 

13 

Dyer  and  Bradley  (1982) 

0.40 

28  1.00 

14 

*  Experiment  considered 

only  heat 

flux  in 

stratification 

effects. 

**  With  k  =  0.40  found 

results  o 

f  Dyer  et 

al.  (1970) 

and 

Businger  et  al. 

(1971)  to  be  statistically  equivalent. 


9b.  As  noted  by  Dyer  (1974),  a  basic  difference  between  the  results  of 
Dyer  and  Hicks  (1970)  and  Businger  et  al.  (1971)  is  in  the  value  of  <  .  A 
comparison  of  these  two  models  for  unstable  stratification  by  Lo  and  McBean 
(1978)  using  the  same  data  base  as  used  by  Businger  et  al.  (1971)  indicated 
consistency  of  results  if  <  »  0.40  .  A  complete  reanalysis  of  these  data  by 
Wierenga  (1980)  indicated  a  different  set  of  constants  if  allowance  is  made 
for  wind-flow  distortion  by  the  instrument  platform.  More  recently,  results 
of  an  independent  experiment  were  reported  by  Dyer  and  Bradley  (1982).  They 
found  constants  which  were  very  similar  to  those  found  by  Wierenga  (1980). 

97.  Uncertainty  in  the  values  of  the  constants  represents  either  an 
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inadequacy  of  z/L  as  a  similarity  parameter  or  difficulties  in  making  accu¬ 
rate  measurements.  Lo  and  McBean  (1978)  note  that  flux  differences  of  20  to 
40  percent  occur  in  comparing  the  models  of  Dyer  and  Hicks  (1970)  and  Businger 
et  al.  (!°-!'  vith  much  of  the  difference  due  to  differences  in  k  .  They 
also  note  that  similar  differences  can  be  incurred  with  a  systematic  5-percent 
error  in  mean  wind  speed  alone  with  no  adjustment  of  temperature  or  flux  mea¬ 
surements.  This  does  not  resolve  the  problem.  It  is  usually  assumed  that  the 
theorv  is  adequate  and  that  resolution  of  the  constants  awaits  more  refined 
experiments.  Recent  applied  works  use  either  Dyer's  (1974)  recommendations 
with  <  =  0.40  (Large  and  Pond  1981,  Geernaert  and  Katsaros  1986)  or  the  rec¬ 
ommended  <  =  0.40  of  Lo  and  McBean  (1978)  with  the  constants  of  Businger 
et  al.  (1971)  (Blanc  1985,  Wu  1986).  Both  approaches  give  similar  results. 

98.  The  advantage  of  these  relationships  is  that  they  can  be  integrated 

to  render  relationships  between  fluxes  and  mean  air  and  wind  properties. 
Modification  of  the  drag  coefficient  can  then  be  estimated  or,  conversely, 
possible  errors  in  results  using  an  uncorrected  drag  coefficient  can  be  esti¬ 
mated.  The  parameters  added  to  the  problem  are  minimal  in  number.  They  are 
mean  air  temperature  T  at  reference  elevation  z  =  ,  mean  air  temperature 

T  at  the  air-sea  interface  (z  =  z  )  ,  mean  humidity  q  at  the  reference 
o  _  o 

elevation,  and  mean  humidity  q^  at  the  air-sea  interface  (z  =  z) .  The  lat¬ 
ter  quantity  is  often  assumed  to  be  the  specific  humidity  corresponding  to 
saturated  air,  i.e.,  100-percent  relative  humidity. 

99.  Noting  in  Table  2  that  all  investigations  yield  a?  =  b0  and  as- 


sutning  z^ 

is  the  same  for  all  variables. 

,  Equations  29, 

30, 

and  31  can  be 

integrated 

using  Equations  32  and 

33  to 

yield 

T  -  T  =  a. 

w'T' 

In  — 

- 

(34) 

o  2 

u* 

z 

0 

t\l/ 

1 

q  -  q  *  a. 

wr  q f 

_ 

In  — 
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and 


U  -  U  =  — 
o  < 


In  —  -  $ 

z  m  »  L  / 


where,  in  stable  stratification  (z/L  >  0) , 


*m(t)  =  "  al  I 


(36) 


(37) 


)  ’  *,(f)  -  5  f 


and,  in  unstable  stratification  (z/L  <  0) , 


♦Jr)  =2  In 

TIT 


1  +  ( 

v 1/4' 

1  -  bi  ?) 

+  In 

1  +  ( 

z\1/2" 
1  -bi  l) 
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(38) 


2  tan 
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1/4 


+  i  (39) 


r(!)  ‘  %(!) 


2  In 


1  +  U 


(l  -  b3  I) 
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(40) 


where 


=  dimensionless  integral  of  $ 

qQ  =  mean  specific  humidity  at  the  air-sea  interface 

=  dimensionless  integral  of 

•p  =  dimensionless  integral  of  $ 
m  m 

The  forms  of  Equations  39  and  40  were  given  by  Paulson  (1970).  Equations  34 
to  40  are  all  coupled  through  the  variable  z/L  .  In  terms  of  measured  prop¬ 
erties  (reference  height,  wind  speed,  water  speed,  temperatures,  and  humidi¬ 
ties)  and  modeled  surface  roughness  zq  ,  this  variable  is 
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o 

where  i_  is  used  in  the  denominator  to  represent  Equations  38  and  40.  Equa¬ 
tion  41  is  transcendental  in  z/L  and  so  must  be  solved  by  iterative  approxi¬ 
mation  or  through  graphs  or  tables.  Algorithms  can  be  v/ritten  readily  to  per¬ 
form  these  calculations  on  a  computer. 

100.  From  Equation  36  a  drag  coefficient  can  be  found  as 


Comparison  with  Equation  6  shows  the  effect  of  stratification,  for  this  model, 
to  be  due  to  the  function  ip  (z/L)  in  the  denominator  of  Equation  42.  Note 
that  transfer  coefficients  for  temperature  and  humidity  can  be  found  from 
Equations  34,  35,  and  42.  They  are  not  considered  here  because  the  primary 
concern  is  with  momentum  transfer. 

101.  To  illustrate  the  importance  of  stratification  in  wind-stress 
estimates  some  computations  were  made  using  the  equations  given  above.  In 
this  problem,  wind  stress  (i.e.,  u^)  was  held  constant  and  air-sea  tempera¬ 
ture  difference  was  varied  from  -10°  C  to  10°  C  .  No  variation  in  humidity 
was  considered.  Estimates  of  10-m  elevation  wind-speed  difference  and  drag 
coefficient  were  then  made.  The  constants  used  in  the  problem  were  those  pro¬ 
posed  by  Wierenga  (1980)  as  given  in  Table  2.  Surface  roughness  was  estimated 
by  Charnock's  model,  Equation  8,  with  a  -  0.02  ,  after  Wu  (1980). 

102.  Results  are  shown  in  Table  3  for  two  values  of  u*  •  A  low  u^ 
(0.20  m/sec,  corresponding  to  the  minimum  surface  wave  phase  speed  as  dis¬ 
cussed  in  Part  II)  and  a  somewhat  higher  u^  (0.90  m/sec)  are  indicated  in 
Table  3.  Results  are  somewhat  dramatic  in  the  low  u^  case.  For  the  temper¬ 
ature  range  shown  the  wind  speed  varies  from  4.9  to  8.2  m/sec  for  the  same 

2 

stress  (pu^) .  The  corresponding  drag  coefficient  varies  by  almost  a  factor  of 
-3  -3 

3  from  0.6  *  10  to  1.7  *  10  .  Conversely,  this  means  that  if 
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Table  3 


Some  Effects  of  Stratification  on  the  Relation  Between  Wind  Stress, 
Wind  Speed,  and  Drag  Coefficient 


u*  = 

0.20  m/sec 

u*  = 

0.90  m/sec 

T  .  -  T 

air  water 

°C 

u  -  u 

r  o 

m/sec 

CD  *  103 

U  -  U 

r  o 

m/ sec 

CD  .  103 

10 

8.2 

0.6 

20.2 

2.0 

5 

7.2 

0.8 

19.6 

2.1 

0 

5.6 

1.3 

19.1 

2.2 

-5 

5.1 

1.5 

18.7 

2.3 

-10 

4.9 

1.7 

18.3 

2.4 

stratification  effects  are  ignored,  then  wind  stress  estimates  can  vary  bv  as 
much  as  100  percent  from  true  wind  stresses  at  these  wind  speeds. 

103.  Less  dramatic  effects  are  seen  for  the  higher  u^  .  In  this  ca=e, 
the  wind-speed  difference  varies  by  about  i0  percent  and  the  drag  coefficient 
by  about  20  percent  over  the  temperature  range  shown.  The  reason  for  this  can 
be  seen  in  the  basic  definition  of  L  following  Equation  28.  For  about  the 
same  heat  flux  an  increase  in  u^  gives  a  larger  magnitude  of  L  and  so  a 
smaller  z/L  at  a  fixed  reference  height.  By  Equations  32  and  33  a  smaller 
z/L  results  in  less  effect  by  stratification.  Wind  speeds  shown  for  the 
higher  u*  are  on  the  order  of  20  m/sec  and  so  are  intermediate  between  low 
wind  speeds  and  hurricane  winds.  At  much  higher  wind  speeds,  effects  of  mod¬ 
erate  temperature  differences  between  sea  and  air  are  negligible  according  to 
this  model. 

104.  These  results  suggest  that  stratification  is  important  only  part 
of  the  time  over  the  full  suite  of  conditions  of  interest  in  coastal  dynamics. 
They  should  not  be  ignored,  however,  for  at  least  two  reasons.  First,  wind 
speeds  less  than  20  ra/sec  (for  which  stratification  is  important)  ire  present 
for  a  far  greater  proportion  of  time  than  are  higher  wind  speeds.  Proper 
modeling  of  day-to-day  wind-driven  waves  and  currents  along  with  processes 
which  depend  on  these  phenomena,  such  as  mean  sediment  transport  rates  and 
long-term  morphologic  change,  cannot  be  considered  accurate  if  the  primary 
forcing  (i.e.,  wind  stress)  is  unknown  to  within  a  factor  of  3.  Second,  even 
in  the  evolution  of  large  storm  systems  the  wind  during  the  initial  stages  of 
wave  growth  and  wind  fields  at  some  distance  from  the  region  of  maximum  storm 
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winds  are  less  than  20  m/sec.  If  the  areal  extent  of  these  lesser  winds  is 
large  relntve  to  the  area  of  intense  winds,  then  the  far  field  winds  tray  do 
significant  amounts  of  work  on  storm-driven  waves  and  currents.  The  rate  of 
doing  work  depends  or  the  intensity  of  the  wind  stress.  An  uncertainty,  bv  a 
factor  of  3,  will  then  have  a  significant  effect  on  total  storm  response. 

105.  These  arguments  indicate  that  vertical  variations  in  air  density 

can  have  a  significant  effect  on  relationships  between  wind  stress  and  wind 
profiles  upon  which  drag  coefficient  formulae  depend.  The  assumptions  made 
for  the  flux-profile  relationships  presented  here  are  among  the  simplest  that 
can  be  made  so  the  corrections  to  must  be  considered  of  first  order  onlv. 

Verification  of  these  effects  by  observations  over  land  indicates  that  it  is 
possible  to  reduce  the  uncertainity  in  Cp  from  about  300  percent  to  about 

30  percent  at  low  to  medium  wind  speeds.  This  is  a  significant  reduction. 

The  remaining  uncertainity  must  be  considered  to  be  the  limit  (due  to  physical 
arguments  or  observation  techniques)  of  current  understanding  of  this  process. 
It  is  suggested  that  a  complete  model  for  wind-driven  waves  and  currents  must 
include  consideration  of  stratification  effects  of  the  type  presented  here. 

106.  There  are  two  other  processes  of  air-sea  mass  exchange  which  ore 
important  dynamically  at  least  part  of  the  time.  Their  effects  sre  beyond  the 
scope  of  this  report,  but  their  existence  is  noted  here  for  completeness.  The 
first  is  precipitation.  Particulate  water  condensed  at  some  elevation  above 
the  sea  surface  and  falling  through  the  atmosphere  accumulates  vertical  momen¬ 
tum  from  gravity  and  horizontal  momentum  from  the  air  through  which  it  falls. 

On  striking  the  sea  surface,  some  of  this  momentum  is  imparted  to  the  sea  sur¬ 
face  either  as  a  direct  input  to  mean  sea  surface  velocity  or  indirectly 
through  the  formation  of  splash  waves,  which  would  increase  surface  roughness. 
The  addition  of  fresh  water  with  a  temperature  different  from  sea  temperature 
would  also  affect  heat  and  humidity  fluxes  in  the  surface  layer  of  the  atmo¬ 
sphere.  The  second  process  relates  to  sea  spray,  of  which  significant  amounts 
are  exchanged  with  the  atmosphere  in  very  high  winds.  Its  concentration  should 
decrease  with  elevation  above  the  sea  surface  so  it  would  have  a  stabilizing 
effect  (reduced  C^)  based  on  gradients  of  total  atmospheric  density.  A  coun¬ 
tering  effect  (increased  C^)  occurs  by  accumulation  of  momentum  in  the  spray 
from  the  wind  and  subsequent  impact  by  the  spray  on  the  sea  surface.  The  net 
effect  of  particulate  water  (rain  and  spray)  in  air-sea  momentum  exchange  is 
not  known.  It  may  be  important  and  so  should  be  considered  in  future  research. 
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107.  The  relationship  between  the  wind  stress  and  the  velocitv  profile, 
given  in  Equation  4,  can  be  modified  if  the  forcing  of  the  system  is  strongly 
time  dependent.  This  is  important  because  momentum  transfer  to  the  sea  sur¬ 
face  is  governed  by  flow  dynamics  (i.e.,  turbulence).  If  the  turbulence  can¬ 
not  keep  pace  with  changes  in  wind  speed,  then  the  velocity  profile  will  be 
modified.  Integration  of  the  true  velocity  profile  will  then  give  velocities 
at  a  given  reference  elevation  which  are  different  from  the  integration  of 
Equation  4.  Use  of  observed  wind  speeds  with  drag  coefficients  based  on  the 
integration  of  Equation  4  will  then  be  incorrect. 

108.  In  steady,  uniform  conditions,  the  primary  force  balance  in  the 
atmospheric  boundary  layer  is  among  Coriolis,  pressure  gradient,  and  fraction 
(stress  divergence)  forces.  In  natural  systems,  the  pressure  gradient  is  sel¬ 
dom  steady.  Fronts,  centers  of  high  and  low  pressure,  and  storm  systems  all 
migrate  in  time  such  that  the  boundary  layer  is  almost  constantly  adjusting  to 
new  conditions.  Such  systems  also  vary  in  space,  which  means  natural  systems 
are  not  completely  horizontally  homogeneous.  A  complete  treatment  of  natural 
systems  requires  modeling  (or  measuring)  all  terms  in  the  governing  equations, 
which  is  beyond  the  scope  of  this  work.  However,  a  simple  scaling  argument 
can  provide  at  least  an  estimate  of  conditions  where  wind  stress  estimates 
from  simple  drag  coefficient  models  can  be  considered  reliable. 

109.  In  typical  numerical  modeling  of  sea  state  evolution  it  is  conven¬ 
ient  to  treat  wind  input  at  each  time-step  as  if  it  were  in  an  equilibrium 
steady-state  force  balance,  even  though  the  mean  wind  changes  from  one  time- 
step  to  the  next.  That  is,  there  can  be  a  local  flow  acceleration  3l'/3t  due 
to  an  imbalance  of  the  other  forcing  terms.  Of  importance  is  the  magnitude  of 
the  acceleration  necessary  to  cast  doubt  on  the  reliability  of  the  velocity 
profile  given  by  Equation  4. 

110.  One  way  to  do  this  is  to  estimate  the  time  scale  of  adjustment  of 
the  atmospheric  boundary  layer  to  a  rapid  change  of  forcing  between  two  steady 
conditions.  While  the  velocity  scale  uA  in  steady  flow  is  identified  with 
the  stress,  it  can  also  serve  as  a  rough  estimation  of  turbulent  propagation 
rate.  That  is,  a  change  of  conditions  at  a  point  in  the  boundary  layer  can  be 
considered  to  diffuse  vertically  away  from  that  point  at  a  velocity  of  about 
u^  .  If  the  boundary  layer  has  a  thickness  6  ,  then  the  time  scale  T  for 


a  property  to  nierate  from  the  boundary  through  the  whole  boundary  layer 
thickness  is.  T  -  :/u.  .  A  typical  boundarv  laver  thickness  is  on  the  order 

of  1,000  r.  For  a  reasonably  small  u^  =  0.2  m/sec,  the  time  scale  is 

T  -  5,000  sec,  which  is  slightly  more  than  1  hr. 

II!.  For  the  shear  (3lT/?z)  in  Equation  4  to  be  unchanged  to,  say,  10  per¬ 
cent,  u^  must  vary  by  less  than  10  percent.  If  the  drag  coefficient  stavs 
about  the  same,  then  the  wind  speed  must  also  change  by  less  than  10  percent. 

If  the  acceleration  is  approximated  as  3U/3t  -  £U/At  with  wind  speed  change 
— U  -  0.1  U  and  time  scale  At  -  T  to  represent  the  tolerable  wind  speed 
change  in  the  time  it  takes  for  the  boundary  layer  to  readjust,  then  3U/3t 
-  0.1  U^/T^  .  That  is,  if  mean  reference  wind  changes  are  on  the  order  of 
10  percent  per  hour  or  less,  then  time  dependence  can  be  neglected  if  allow¬ 
ance  is  made  for  a  20-percent  uncertainty  level  in  drag  coefficient  estimates. 

112.  Note  that  over  the  sea  surface,  a  change  in  wind  stress  results  in 

a  change  in  sea  state.  This  changes  the  surface  roughness  and,  by  Equation  5, 
the  velocity.  This  effect  is  much  less,  however,  than  the  change  in  u* 

itself.  This  can  be  seen  bv  using  Charnock’s  model  for  z  (=  au^/g)  .  A 

o  * 

10-percent  change  in  u  changes  z  by  about  20  percent.  Using  Equation  5, 

*  o 

the  ratio  of  (u  /<)  ln(z/z  )  to  (u^/f)  ln(z/1.2  z  )  is,  with  some  manip- 
*  o  *  o 

ulation  of  these  expressions,  about  1  +  In ( 1 . 2) /ln(z fz  )  .  This  differs  from 

o 

1  by  approximately  1  percent  for  z  =  10  m  and  z^  =  0.01  m.  The  difference 

is  less  for  smaller  z  .  Thus,  the  primarv  effect  of  a  change  in  U  of 

o 

10  percent  is  a  change  in  u^  by  about  10  percent  with  an  additional  change 
of  about  1  percent  or  less  due  to  change  in  roughness.  This  argument  is  sub¬ 
ject  to  limitations  on  the  applicability  of  Charnock's  Zq  and  the  assumption 
that  sea  state  keeps  pace  with  changes  in  wind  stress.  If  sea  state  lags  a 
change  in  wind  speed  (as  it  does  because  wind  stress  must  work  on  the  sea  sur¬ 
face  over  a  finite  time  to  build  or  modify  a  sea),  then  change  in  roughness 
will  be  more  gradual  for  moderate  changes  in  wind  speed.  Dependence  of  z^ 
on  C  / u.  ,  as  suggested  in  the  models  of  Hsu  (1974)  and  Kitaigorodskii 
(1973),  is  also  affected  slightly.  A  10-percent  change  in  u*  gives  about 
10-percent  change  in  C  /u.  .  According  to  the  work  of  Huang  et  al.  (1986), 
which  incorporates  effects  of  both  these  models,  a  10-percent  change  in  0^/u^ 
gives  at  most  about  an  equivalent  change  of  10  percent  in  z0^an  •  By 
logarithmic  scaling  above,  the  difference  is  still  on  the  order  of  1  percent 
or  less  in  wind  speed. 
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113.  The  boundary  layer  must  also  adjust  if  the  wind  changes  direction. 
The  magnitude  of  acceptable  wind  direction  change  can  be  estimated  by  assurair.c 
the  wind  sreed  constant  and  a  change  in  direction  such  that  the  wind  component 
in  the  original  direction  varies  by  less  than  10  percent.  This  corresponds  to 
an  angle  with  a  cosine  of  0.9  or  about  25  deg.  If  wind-stress  direction  fol¬ 
lows  wind  direction,  then  surface  roughness  effects  are  of  the  same  magnitude 
as  given  above. 

114.  Through  these  rough  scaling  arguments  it  appears  that  wind-stress 
estimates  of  about  20-percent  accuracy  can  be  attained  for  hour-to-hour  wind 
speed  changes  of  less  than  10  percent  and  wind  direction  changes  of  less  than 
25  deg.  If  the  mean  wind  changes  by  more  than  this,  then  larger  errors  will 
occur.  The  gross  trend  of  the  error  can  be  anticipated,  but  its  magnitude 
must  rely  on  observations  in  the  absence  of  a  complete  turbulence  theorv. 

115.  The  argument  concerning  the  error  trend  is  based  on  the  relative 
rates  of  growth  and  decay  of  the  turbulence  which  couples  wind  shear  and 
stress.  In  conventional  turbulence  theory  (Tennekes  and  Lumley  1972)  fluctu¬ 
ating  flow  properties  exist  over  a  broad  range  of  frequencies  and  wave  num¬ 
bers.  Turbulence  is  produced  at  large  scales  (i.e.,  sizes  proportional  to 
distance  from  the  boundary)  and  dissipated  at  small  scales  where  the  fluctu¬ 
ating  strain  rate  is  high  enough  for  molecular  viscosity  to  convert  turbulent 
energy  to  heat.  The  sizes  of  these  small-scale  motions  can  be  estimated  from 
Kolmogorov's  scaling  argument  (Tennekes  and  Lumley  1972),  which  asserts  that 

the  length  scale  is  based  on  the  dissipation  rate  t  and  kinematic  viscosity 

3  1/4 

v  .  It  takes  the  form  (v  / z)  .  In  steady,  neutral  flow  it  is  often  as¬ 

sumed  that  the  rate  of  production  of  turbulence  equals  the  rate  of  dissipa¬ 
tion.  This  argument  was  noted  above  and  resulted  in  Equation  28  wherein 

3 

€  =  u*/<z  .  Using  this  with  typical  values  of  atmospheric  variables  yields  a 
dissipation  length  scale  on  the  order  of  1  mm  at  an  elevation  of  z  *  10  m. 
This  is  much  less  than  the  production  length  scale  which  is  on  the  same  order 
as  2  ,  i.e.,  10  u. 

116.  Turbulence  length  scales  intermediate  between  production  and  dis¬ 
sipation  scales  are  assumed  to  relate  to  fluid  motions  interacting,  through 
inertial  processes,  with  motions  at  nearby  scales.  In  this  case,  large-scale 
motions  (or  eddies)  decay  to  small-scale  motions  at  a  rate  governed  by  the 
magnitude  of  the  dissipation  term,  but  not  by  direct  viscous  processes.  The 
variances  at  wave  numbers  associated  with  these  length  scales  are  said  to  be 
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in  the  inertial  subrange  of  the  turbulence  spectrum.  Because  of  the  separa¬ 
tion  in  scales  of  production  and  dissipation,  .nere  is  a  finite  time  required 
for  large  eddies  to  decay  to  heat.  Production,  on  the  other  hand,  is  governed 
by  the  mean  shear  and  the  stress  in  accordance  with  the  left  side  of  Equa¬ 
tion  20.  If  the  mean  shear  changes,  then  the  rate  of  production  changes  imme¬ 
diately  even  if  the  stress  does  not  change.  Once  produced,  the  turbulent 
eddies  reside  at  intermediate  scales  as  they  decay  to  scales  small  enough  to 
be  dissipated  by  viscous  processes.  Since  this  takes  a  finite  time,  there 
will  tend  to  be  residual  turbulence  even  if  the  production  term  goes  to  zero. 

117.  The  effect  in  the  surface  layer  of  the  time-dependent  boundary 
layer  is  to  adjust  rapidly  to  flow  acceleration,  but  more  slowly  to  flow  decel 
eration.  If  the  pressure  gradient  forcing  the  flow  increases,  the  flow  tends 
to  accelerate.  Because  the  flow  at  the  boundary  must  move  at  the  boundary 
velocity  it  does  not  accelerate  as  rapidly.  This  increases  the  shear  at  the 
boundary,  which  increases  turbulent  production  of  eddies.  These  act  to  ex¬ 
change  momentum  with  the  boundary.  Hence,  the  stress  increases  at  the  bound¬ 
ary.  This  retards  the  flow  just  above  the  boundary,  increases  the  shear  there 
and  thereby  enhances  momentum  exchange.  This  readjustment  moves  outward  from 
the  boundary  at  a  rate  which  can  be  estimated  roughly  by  u^  and  reaches  ref¬ 
erence  elevation  z  in  time  z  / u.  .  For  z  =  10  m  and  u.  -  0.2  m/sec, 

r  r  *  r  * 

this  is  about  50  sec  which  is  rapid  relative  to  the  readjustment  time  T^_  of 
the  whole  boundary  layer  as  discussed  above.  An  anemometer  at  10-m  elevation 
would  sense  a  high  velocity  relative  to  the  new  boundary  stress  for  about  a 
minute,  and  would  then  come  into  adjustment.  If  averaged  for  a  time  much 
longer  than  a  minute,  the  anemometer  output  would  be  only  slightly  high  -"'r 
the  new  equilibrium  boundary  stress,  based  on  a  steady  drag  coefficient  model. 

118.  On  the  other  hand,  if  the  pressure  gradient  forcing  the  flow  de¬ 
creases,  then  the  flow  decelerates  and  the  shear  at  the  boundary  decreases. 
Though  local  shear  production  drops,  the  residual  turbulence  from  the  previous 
higher  velocity  flow  continues  to  act  to  transfer  momentum  to  the  boundary. 
Until  the  residual  eddies  equilibrate  to  the  new  flow  condition,  the  boundary 
stress  wil1  be  higher  than  that  deduced  from  the  instantaneous  velocity. 

Since  eddies  of  the  size  of  reference  elevation  from  the  boundary  are  dis¬ 
tributed  through  the  upper  boundary  layer  flow,  the  time  scale  of  readjustment 
Is  likely  to  be  as  long  as  T^_  deduced  above  for  readjustment  of  the  whole 
boundary  layer.  Evidence  of  this  behavior  was  provided  in  the  extreme  case  of 
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oscillating  boundary  layer  flow  by  Long  C1981)  in  an  analysis  of  water  tunnel 
data  reported  by  Jonsson  and  Carlsen  (1976).  In  this  case,  the  boundary 
stress  behaved  as  if  it  were  the  product  of  one  u*  based  on  the  ensemble 
averaged  instantaneous  mean  shear  (i.e.,  from  Equation  4)  and  another  u^ 
based  on  the  maximum  mean  shear.  At  times  when  the  instantaneous  mean  shear 
was  20  percent  of  the  maximum  shear,  the  stress  estimated  from,  the  velocity 
profile  (i.e.,  Equation  5)  was  low  by  a  factor  of  5.  This  means  that  the 
scaling  which  led  to  Equation  4  is  incorrect  for  this  case  and  results  deduced 
from  it  yield  significant  errors. 

119.  While  natural  atmospheric  conditions  are  seldom  this  extreme,  the 
results  of  the  arguments  presented  here  suggest  that  some  caution  should  be 
exercised  in  application  of  drag  coefficients  to  rapidly  changing  flow  condi¬ 
tions.  Data  from  experiments  designed  to  obtain  sequences  of  wind  stress  and 
mean  velocity  should  be  examined  for  biases  in  drag  coefficient  estimates  dur¬ 
ing  accelerating  and  decelerating  flow.  This  can  be  done  with  data  from  an 
experiment  like  that  described  in  Part  V  of  this  report. 

Horizontal  Nonuniformity 

120.  A  flow  is  horizontally  nonuniform  if  the  forcing  or  boundary  con¬ 
ditions  vary  in  any  horizontal  direction.  In  such  a  flow,  properties  such  as 
momentum  or  mass  can  accumulate  or  be  depleted  in  a  control  volume  by  horizon¬ 
tal  divergences  of  these  properties.  In  this  case  the  flux  of  momentum  is  not 
purely  in  the  vertical  direction,  and  so  violates  the  assumptions  leading  to 
Equation  4. 

121.  In  the  case  of  air  flow  over  coastal  regions,  conditions  are  never 
completely  uniform.  The  primary  cause  of  nonuniformity  in  such  a  system  is 
the  differences  between  wind  flowing  over  land  and  wind  flowing  over  water. 

Land  has  surface  roughness  elements  which  are  generally  fixed,  of  a  broad 
variety  of  scales,  and  usually  different  from  those  of  water.  Land  can  have 
isolated,  large-scale  features  which  distort  air  flow  in  their  vicinity  and 
thereby  cause  local  inhomogeneities.  Sources  of  moisture  and  heat  are  differ¬ 
ent  over  land  from  those  over  water.  This  means  that  stratification  effects 
which  distort  velocity  profiles  can  cause  horizontal  variations  in  mean  veloc¬ 
ity.  It  also  means  that  differential  acceleration  of  air  due  to  varying  den¬ 
sity  (thermal  wind)  may  be  important.  These  features  are  in  addition  to  normal 
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horizontal  variability  of  large-scale  weather  systems  which  drive  the  wind. 

122.  Another  effect  which  may  affect  horizontal  uniformity  over  water 
is  the  change  in  roughness  of  the  sea  surface  due  to  shoaling  effects.  That 
is,  if  shoalir.g-induced  wave  direction  and  steepness  change  induce  strong 
horizontal  gradients  in  zq  ,  then  there  will  be  a  tendency  for  the  mean  wind 
to  vary  significantly. 

L23.  It  is  important  to  obtain  an  idea  of  the  areal  extent  of  varia¬ 
tions  in  flow  properties  which  can  lead  to  significant  errors  in  drag  coeffi¬ 
cient  estimates  of  wind  stress.  It  is  also  useful  to  anticipate  the  nature  of 
deviations  so  that  such  effects  can  be  accounted  for  in  observations. 

124.  Terms  in  the  governing  equations  which  become  important  in  non- 
uniform  flow  are  the  advective  terms  of  momentum,  mass,  and  turbulence  proper¬ 
ties  along  with  any  added  pressures  necessary  to  divert  flow  around  obstacles. 
A  prominant  advective  term  in  the  mean  wind  direction  is  U  3U/3x  which  rep¬ 
resents  the  apparent  acceleration  due  to  advection  of  a  wind  field  that  varies 
in  the  x-direction.  To  estimate  the  size  this  term  acquires  before  signifi¬ 
cant  deviations  of  Equation  5  are  incurred,  a  rough  estimate  is  made  of  the 
horizontal  distance  required  for  the  boundary  layer  to  adapt  to  a  change  in 
conditions  at  the  surface. 

125.  In  the  section  on  time  dependence,  a  time  scale  T  =  <$/u.  was 

introduced  as  an  estimate  for  temporal  readjustment  to  a  uniform  change  of 
surface  condition.  If  a  change  occurs  in  space,  then  the  new  conditions  will 
migrate  outward  from  the  boundary,  as  before,  but  will  also  be  advected  down¬ 
wind.  If  is  a  scale  for  wind  speed,  then  the  boundary  layer  adjusts  in  a 

distance  of  Ax  -  U  T  =  U  5/u.  »  6 / /cl~  .  For  6  on  the  order  of  1  km  and 
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Cp  in  the  range  1  *  10  to  4  *  10  ,  then  Ax  is  on  the  order  of  15  to 

30  km.  If  it  is  required  that  wind  speed  varies  by  less  than  10  percent  over 

this  distance  such  that  AU  <  0.1  U  ,  then  U  3U/3x  -  U  AU/Ax  <  0.1  U  /T 
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Note  that  this  is  the  same  relationship  as  was  found  for  time  dependence.  For 
a  reference  wind  of  10  m/sec,  this  means  that  the  wind  must  vary  by  less  than 
1  m/sec  in  15  to  30  km. 

126.  While  variations  less  than  1  m/sec  may  occur  over  the  open  ocean, 
this  is  a  rather  extreme  requirement  for  coastal  regions.  A  distance  of  15  to 
30  km  encompasses  the  width  of  the  coastal  region  in  many  locations.  Resio 
and  Vincent  (1977),  in  a  study  of  winds  over  the  Great  Lakes,  noted  that  wind 
speeds  over  water  were  often  20  percent  greater  than  wind  speeds  over  land. 
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From  this,  it  would  appear  that  coastal  regions  are  frequently  subject  to  ex¬ 
cessive  nonuniformity.  While  this  is  true  locally,  the  effect  may  not  be  as 
extensive  as  this  gross  scaling  would  suggest.  For  instance,  if  an  onshore 
wind  is  uniform  up  to  the  beach  such  that  boundary  layer  adjustments  occur 
only  inland,  then  the  observations  of  Resio  and  Vincent  (1977)  would  be  sup¬ 
ported  and  the  air  flow  over  the  sea  would  be  uniform.  This  is,  in  fact,  sug¬ 
gested  by  the  Shore  Protection  Manual  (SPM)  (1984)  which  recommends  that  wind 
measurements  immediately  adjacent  to  a  water  body  are  the  same  as  winds  over 
water  for  onshore  flow.  However,  if  sea  state  varies  as  the  beach  is  ap¬ 
proached  or  if  the  boundary  layer  begins  to  adjust  upwind  of  the  beach,  then 
this  hypothesis  may  be  overly  simple. 

127.  Several  investigators  have  made  valuable  theoretical  and  empirical 
investigations  into  the  relations  between  winds  measured  over  land  and  winds 
that  exist  over  water.  Among  these  are  Resio  and  Vincent  (1977)  and  Hsu 
(1981,  1986).  However,  knowledge  of  the  wind  speed  and  direction  over  water 
is  not  sufficient  to  estimate  accurately  the  wind  stress  over  water,  if  the 
drag  coefficient  is  not  known  or  is  applied  in  a  region  where  it  is  invalid. 
The  intent  here  is  to  identify  regimes  wherein  the  drag  coefficient  deduced 
from  Equation  4  (or,  as  corrected  for  stratification,  from  Equation  29)  is 
approximately  valid.  In  these  cases,  reasonable  wind-stress  estimates  can  be 
expected  from  knowledge  of  wind  speed,  wind  direction,  sea  state,  and  mean 
vertical  density  change,  as  discussed  above.  In  cases  where  the  simple  model 
does  not  apply,  wind  stress  will  be  related  to  a  larger  set  of  parameters. 

128.  In  analogy  with  arguments  given  above  for  time-dependent  flow,  it 
is  possible  to  assume  that  the  atmospheric  surface  layer  adjusts  more  rapidly 
to  varying  conditions  than  the  whole  boundary  layer  because  of  the  close 
proximity  of  the  boundary.  In  this  case,  the  time  scale  for  altered  condi¬ 
tions  at  the  surface  to  reach  an  elevation  z  »  z  is  z  /u^  .  The  advective 
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length  scale  for  reference  mean  velocity  U  is  U  z  /u^  *  z  //C~  •  For  a 

r  r  r-3  r  D  _3 

10-m  reference  elevation  and  in  the  range  1  *  10  to  4  x  10  ,  this 

distance  is  on  the  order  of  150  to  300  m  (or  300  to  600  m  for  a  reference  ele¬ 
vation  of  20  m) .  The  constraint  of  less  than  10-percent  variation  in  velocity 
over  distance  of  this  order  requires  variation  of  less  than  10  percent  in  u^ 
or  z^/L  (for  the  stable  case,  i.e.,  the  most  sensitive)  and  roughly  less  than 

a  factor  of  2  in  z  for  cases  wherein  only  one  of  these  variables  is 
o 

changing. 
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129.  These  constraints  are  much  less  stringent  than  those  associated 
with  readjustment  of  the  whole  boundary  layer.  It  seems  true  intuitively  that 
the  constraints  are  met  at  least  part  of  the  time  in  open  coastal  waters.  Ex¬ 
ceptions  would  occur  for  very  rapid  spatial  variation  in  forcing  cr  rapid 
change  in  bottom  topography  which  leads  to  strong  shoaling  effects  (wave 
steepening,  refraction)  over  distances  short  compared  with  the  advective 
length  scale.  As  with  time-dependent  flow,  differences  would  be  expected  in 
transitions  over  decreasingly  rough  surfaces  compared  with  increasingly  rough 
surfaces.  If  roughness  increases  in  space,  the  turbulence  is  expected  to  keep 
pace  (i.e.,  be  locally  in  balance)  since  production  rates  are  increasing.  In 
cases  of  decreasing  roughness,  relict  turbulence  from  upwind  is  expected  to 
couple  the  wind  more  strongly  to  the  sea  surface  so  that  the  actual  stress  is 
higher  than  that  deduced  from  a  wind-speed-dependent  drag  coefficient. 

130.  Nearer  to  the  coast,  within  several  advective  length  scales  of  the 
beach,  the  flow  may  become  more  complicated.  The  surface  roughness  parameter 
z^  can  increase  by  orders  of  magnitude  over  land  as  compared  with  over  water. 
If  roughness  elements  on  land  are  not  themselves  on  the  order  of  the  reference 
elevation  z^  ,  then  the  problem  is  simplified  somewhat.  A  primary  variable 
in  this  case  is  the  direction  of  the  wind  relative  to  the  beach.  Extreme 
cases  are  winds  onshore  and  offshore.  For  onshore  winds,  the  air  flow  is 
retarded  as  it  progresses  inland.  While  some  effect  is  expected  upwind  of  the 
beach,  it  is  not  expected  to  extend  very  far  seaward.  The  reason  for  this  is 
that  the  roughness  over  land  (in  this  case)  is  not  drastically  different  from 
that  over  water.  If  shoaling  effects  are  not  too  abrupt,  it  can  be  expected 
that  onshore  winds  can  be  considered  uniform  up  to  the  beach. 

131.  If  the  wind  is  offshore  (i.e.,  air  moving  seaward),  then  some 
variation  in  wind  stress  over  water  will  occur.  This  is  because  turbulence 
from  upwind  is  advected  over  water.  This  will  be  complicated  by  sea  state. 

If  there  are  no  waves,  then  the  air  will  tend  to  accelerate  over  the  water. 

If  a  swell  is  running  and  zq  depends  on  wave  amplitude  and  phase  speed,  then 
it  is  possible  for  the  sea  to  be  more  rough  than  land.  In  this  case,  an  off¬ 
shore  wind  will  decelerate  over  water.  Such  flows  can  still  be  considered 
uniform  if  the  change  in  roughness  does  not  incur  a  change  in  u^  or  the 
characteristic  fluctuating  air  density  by  more  than  10  percent. 

132.  If  these  criteria  are  not  met,  then  conditions  can  be  considered 
nonuniform.  Wind  stress  and  velocity  profiles  will  not  scale  together. 
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Accurate  estimates  of  wind  stress  must  then  rely  on  more  elaborate  turbulence 
theory  along  with  empirical  correlations  among  measured  field  variables. 

133.  If  the  primary  interest  is  in  generation  by  wind  of  nearshore 
waves  and  currents,  then  nonuniformity  of  wind  blowing  offshore  may  not  be 
important.  Fetch  and  water-depth  limitations  result  in  much  smaller  wave  and 
current  generation  close  to  a  beacn  by  winds  blowing  offshore  than  by  winds 
blowing  onshore.  If  offshore  blowing  winds  are  neglected  in  the  interest  of 
high-energy,  wind-induced  nearshore  dynamics,  then  the  problem  is  simplified 
somewhat.  This  is  done  here. 

134.  For  onshore  blowing  winds,  there  are  two  additional  conditions 
which  can  affect  the  assumption  of  uniformity  near  a  beach.  The  first  is  the 
condition  wherein  physical  roughness  elements  on  land  have  a  height  scale  much 
larger  than  nearshore  waves.  Sand  dunes,  buildings,  trees,  and  mountainous 
cliffs  are  elements  of  this  type.  A  wind  encountering  such  elements  from  a 
smoother  environment  is  diverted  over  or  around  them  by  pressure  (form  drag) 
forces  on  their  upwind  faces.  The  modification  of  mean  streamlines  is  af¬ 
fected  some  distance  upwind  due  to  the  incompressible  nature  of  air  at  low 
Mach  numbers.  In  analogy  with  flow  around  an  isolated  solid  (Batchelor  1970), 
the  streamline  diversion  can  be  assumed  to  exist  at  least  three  roughness 
heights  upwind.  Blanc  (1983)  cites  work  showing  roughly  25  percent  deviation 
of  mean  horizontal  wind  speed  at  least  one  roughness  elevation  upwind  of  the 
toe  of  an  idealized  beach  with  a  1:4  slope.  Streamlines  diverted  upward  will 
advect  momentum  out  of  the  surface  layer  and  thus  reduce  shear-induced  momen¬ 
tum  transfer  to  the  sea  surface. 

135.  The  seaward  effect  of  land  roughness  on  horizontal  uniformity  over 
water  depends  on  land  roughness  height  (by  this  argument) .  Land  of  5-m  relief 
will  not  block  the  wind  very  far  seaward.  A  stand  of  30-m  trees  will  have  a 
direct  influence  roughly  100  m  offshore.  Low  mountains  of  500-m  relief  will 
affect  the  flow  more  than  1  km  offshore.  From  this  it  can  be  seen  that  con¬ 
sideration  of  nearshore  land  relief  can  be  of  importance  in  estimating  near¬ 
shore  wind  stress  during  times  of  onshore  winds. 

136.  The  second  condition  wherein  nonuniformity  may  be  important  is 
related  to  the  density  field.  If  the  vertical  air  density  distribution  over 
water  is  different  from  that  over  land,  then  a  pressure  gradient  is  induced 
across  the  intervening  space  due  to  the  unequal  weights  of  the  two  air  col¬ 
umns.  This  baroclinic  condition  occurs  frequently  near  the  shore  because  the 
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boundary  conditions  for  temperature  and  humidity  can  be  very  different  for 
land  as  compared  with  the  sea.  A  common  cause  of  this  is  a  differential  re¬ 
sponse  of  land  and  sea  surfaces  to  solar  radiation.  In  this  case,  wind  aris¬ 
ing  due  to  differential  pressure  is  called  sea  breeze  (or  land  breeze,  depend¬ 
ing  on  direction).  It  is  usually  transient,  acting  over  less  time  than  is 
required  for  geostrophic  equilibrium,  so  winds  tend  to  be  directed  down  the 
pressure  gradient,  i.e.,  cross  shore.  If  the  phenomenon  lasts  for  more  than  a 
few  hours,  then  Coriolis  effects  will  become  important,  and  wind  direction 
will  vary. 

137.  If  density  differences  result  in  forcing  that  is  uniform  with 
depth,  then  the  effect  is  no  different  than  a  barotropic  pressure  field  in  the 
surface  layer.  As  long  as  horizontal  variations  in  barotropic  pressure  gra¬ 
dient  are  less  than  20  percent  in  about  300  m  (as  discussed  above) ,  then  the 
flow  can  be  considered  uniform.  However,  horizontal  variations  in  density 
generally  cause  vertical  variations  in  pressure  gradient.  The  differential 
forcing  in  the  vertical  can  cause  shear  (3U/3z)  in  the  wind  independent  of  any 
boundary  layer  effect.  If  this  happens  in  the  atmospheric  surface  layer,  then 
reference  wind  speed  can  be  altered  without  changing  surface  stress.  If  the 
change  is  significant,  then  reference  wind  coupled  with  a  drag  coefficient 
given  by  Equation  42  will  give  an  incorrect  estimate  of  stress. 

138.  The  magnitude  of  this  effect  can  be  estimated  from  a  simple  force 
balance  wherein  flow  acceleration  is  caused  by  a  pressure  gradient  in  the 
absence  of  friction.  In  this  case, 


3U  _  1  3£ 

3t  *  ~  p  3x 


(43) 


If  the  pressure  p  is  hydrostatic  then  P  *  PQ  -  Pgz  where  p^  is  surface 
pressure  (i.e.,  at  z  *  0) .  The  horizontal  gradient  of  this  has  two  parts, 
one  due  to  variation  in  surface  pressure  and  the  second  due  to  horizontal 
variation  in  density.  The  magnitude  of  the  second  part  is  of  interest  here. 
This  part  causes  Equation  43  to  take  the  form 


HI  »  &£  l£.  gz  3T 

3t  p  3x  T  3x 


(44) 
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where  density  variations  are  assumed  to  arise  principally  from  temperature 
variations . 

139.  Equation  44  is  z-dependent,  indicating  that  in  time  increment  A t 

the  velocity  change  AU  at  elevation  z  depends  on  gravity,  a  reference 

(absolute)  temperature  and  the  horizontal  temperature  gradient.  Normal  mean 

horizontal  temperature  gradients  over  land  and  over  the  open  ocean  are  small 

enough  to  be  neglected.  Near  a  coast,  however,  they  may  become  large  enough 

to  be  significant.  For  example,  Hsu  (1970),  in  studies  of  the  Texas  coast  of 

the  Gulf  of  Mexico,  found  horizontal  temperature  gradients  on  the  order  of 

1°  K/km  to  persist  for  up  to  6  hr.  For  T  =  300°  K  ,  z  =  z  =  10  m  and 

r 

At  =  3  hr.  Equation  44  indicates  AU  =  3  m/sec.  That  is,  at  z  =  10  m,  the 
wind  speed  is  3  m/sec  different  from  the  wind  speed  at  the  surface  in  the 
absence  of  any  other  forcing. 

140.  The  change  is  significant.  For  a  normal  wind  speed  of  10  m/sec 
and  a  stress  unaffected  by  baroclinity,  the  wind  speed  could  vary  from  7  to 
13  m/sec  (depending  on  direction  of  the  temperature  gradient).  Estimated 
stress  based  on  this  range  of  wind  speeds  would  vary  by  the  ratio  of  the 
squares  of  these  speeds,  i.e.,  a  factor  of  more  than  3.  Relative  to  the 
actual  stress,  the  range  of  error  is  from  about  50  percent  low  to  about 

60  percent  high. 

141.  This  effect  is  expected  to  occur  only  during  times  of  low  wind 
speeds  and  high  land-sea  temperature  differences.  At  higher  wind  speeds,  the 
percentage  of  difference  is  less.  Furthermore,  increased  turbulence  would 
tend  to  equilibrate  land  temperature  more  rapidly. 

142.  Several  mechanisms  have  been  proposed  here  which  can  violate 
strongly  the  assumption  of  horizontal  uniformity  necessary  for  the  stability- 
corrected  logarithmic  velocity  profile  to  relate  directly  to  surface  stress. 
From  rough  scaling  arguments  it  appears  that,  if  fluxes  and  surface  roughness 
do  not  change  significantly  over  distances  of  a  few  hundred  metres,  then  the 
assumption  of  horizontal  uniformity  is  justified.  Variations  from  this  are 
possible  for  offshore  winds  over  distances  of  several  tens  of  kilometres, 
onshore  winds  in  regions  of  high  relief,  rapidly  changing  bottom  bathymetry, 
or  irregular  coastlines  at  scales  on  the  order  of  1  km.  Variations  can  also 
occur  for  winds  of  low  to  moderate  speed  in  regions  where  large  horizontal  air 
density  changes  exist.  These  effects  are  in  addition  to  errors  due  to  large 
horizontal  gradients  in  forcing  which  can  occur  in  intense  storms. 
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143.  The  arguments  given  here  are  physically  meaningful  but  largely 
speculative.  If,  as  indicated  in  the  introduction,  wind  stress  is  an  impor¬ 
tant  forcing  term  in  nearshore  wave  and  current  generation,  then  it  seems 
imperative  that  over-water  wind-stress  measurements  be  made  within  a  few 
hundred  metres  of  a  beach.  Drag  coefficients  computed  from  such  observations 
can  then  be  correlated  with  wind  direction,  land-sea-air  density  difference, 
land  relief,  and  nearshore  bathymetry  to  deduce  the  viability  of  various  drag 
coefficient  formulations.  An  experiment  of  the  type  described  in  Part  V  of 
this  report  can,  at  least  partially,  address  this  matter. 
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PART  IV:  METHODS  OF  ESTIMATING  WIND  STRESS  FROM  MEASUREMENTS 


144.  In  the  course  of  making  calculations  of  wind-generated  waves  and 
currents  it  is  generally  inefficient,  if  not  impossible,  to  resolve  the  wind 
field  at  all  space  and  time  scales  necessary  to  model  directly  the  fluctuating 
velocities  which  constitute  turbulent  stresses.  This  step  can  be  bypassed  if 
wind  stress  is  found  to  be  a  unique  function  of  other  integral  parameters 
which  can  be  measured  or  modeled  more  directly.  For  instance,  if  wind  stress 
were  directly  proportional  to  square  of  mean  wind  speed  (at  some  reference 
elevation),  then  Taylor's  (1916)  model  would  be  valid.  A  single  experiment  of 
high  quality  wherein  wind  speed  and  wind  stress  are  measured  simultaneously 
would  then  determine  the  universal  constant  coefficient  of  proportionality 
(0^)  ,  and  t..e  problem  would  be  solved. 

145.  However,  the  references  cited  heretofore  in  this  report  suggest 
that  Cp  is  not  constant  for  all  wind  speeds  and  is  not  independent  of  other 
conditions  such  as  surface  roughness  and  atmospheric  stability.  This  means 
that  further  experimental  and  theoretical  work  is  necessary  to  find  the  cor¬ 
relation  between  CQ  and  the  variables  upon  which  it  depends.  Much  work  has 
been  done  on  this  subject.  Some  of  it  has  been  discussed  in  this  report  in 
the  framework  of  simple  dimensional  arguments  and  supporting  observations. 

The  particular  model  presented  here  for  CQ  over  water  is  not  unique,  but  is 
as  consistent  with  laboratory  and  land-based  observations  as  any  other  set  of 
model  equations.  It  is  reasonable  to  assume  that  it  applies  in  the  marine 
environment . 

146.  Verification  of  hypotheses  embedded  in  the  model  relies  on  accu¬ 

rate  and  complete  observations  of  prototypical  conditions.  At  a  minimum,  this 
means  that  observations  must  be  made  of  wind  stress  and  wind  speed  so  that 
can  be  determined  from  the  most  primitive  model.  In  addition,  certain  ancil¬ 
lary  parameters  must  be  measured  to  support  or  reject  hypotheses  of  dependence 
of  on  these  parameters.  These  hypotheses  include  such  things  as  air-sea 

temperature  and  humidity  differences,  water  surface  mean  velocity,  and  sea 
state  parameters  as  outlined  above.  Verification  of  these  hypotheses  also 
requires  checking  to  confirm  that  assumptions  concerning  the  model,  for  exam¬ 
ple,  steadiness  and  horizontal  uniformity,  are  realized  in  the  prototype.  If 
not,  the  model  must  be  allowed  to  tolerate  some  error  to  compensate  for  these 
effects  or  must  be  augmented  to  account  for  these  effects. 
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147.  Of  the  variables  mentioned  here,  some  of  the  most  difficult  to 
measure  are  the  turbulent  stresses,  or,  more  generally,  the  turbulent  fluxes. 
It  is  because  of  this  difficulty  that  drag  coefficient  formulae  exist;  if 
stresses  can  be  expressed  in  terms  of  variables  which  are  simple  to  measure, 
then  direct  stress  measurements  are  not  necessary.  On  the  other  hand,  model 
verification  depends  critically  on  accurate  stress  measurements.  Histori¬ 
cally,  part  of  the  problem  in  fixing  model  constants  or  interdependency  rela¬ 
tionships  has  been  related  to  inaccuracies  in  stress  measurements.  In  this 
section,  several  common  methods  of  direct  or  indirect  stress  measurement  are 
described  along  with  some  constraints  and  hazards  associated  with  their  use. 

148.  Attention  is  focused  here  on  measurements  made  within  the  atmo¬ 
spheric  surface  layer  (i.e.,  the  lowest  50  m  or  so  of  the  atmosphere).  This 
is  because  of  the  relative  ease  of  use,  generally  lower  cost,  and  higher 
availability  of  such  data.  Platforms  from  which  these  measurements  can  be 
made  include  surface  ships,  buoys,  and  towers.  Other  platforms,  such  as  bal¬ 
loons,  aircraft,  and  satellites,  are  further  removed  from  the  surface  layer 
prrblem,  ua</e  higher  costs,  and  generally  provide  less  data.  These  are  not 
considered  here. 

149.  Major  problems  associated  with  instrumentation  are  twofold.  First 
is  the  adequacy  of  an  instrument  used  to  resolve  the  necessary  space  and  time 
scales  of  the  phenomenon  being  observed.  Second  is  proper  mounting  or  place¬ 
ment  of  an  instrument  such  that  the  sensing  volume  primarily  detects  the  de¬ 
sired  variable  and  not  one  that  has  been  distorted  by  the  mounting  platform  or 
instrument  frame.  Both  problems  are  important  in  measurements  associated  with 
wind-stress  estimation  and  are  considered  here. 

150.  There  are  four  common  methods  of  inferring  wind  stress  from  mea¬ 
sured  variables.  They  are  called  the  bulk  method,  the  profile  method,  the 
direct  method,  and  the  turbulence  inertial  subrange  (or,  sometimes,  dissipa¬ 
tion)  method.  All  are  based  on  arguments  given  within  this  report. 

Bulk  Method 


151.  The  bulk  method  is  the  primary  topic  of  this  report.  It  assumes  a 
relation  between  wind  stress  and  the  square  of  the  air-sea  surface  speed  dif¬ 
ference  through  a  drag  coefficient  which  is  a  function  of  other  mean  variables 
affecting  the  physics  of  air-sea  interaction.  In  the  model  presented  here  the 
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drag  coefficient  is  given  by  Equation  42.  For  a  given  reference  elevation,  it 

deoends  on  the  surface  roughness  scale  z  and  the  stratification-related 

o 

length  scale  L  under  the  assumptions  that  the  flow  is  steady  and  uniform. 

The  roughness  length  can  be  estimated  from  the  models  of  Charnock  (1955), 
Kitaigorodskii  (1973),  or  Hsu  (1974)  but  may  depend  on  a  more  detailed  de¬ 
scription  of  sea  state,  as  discussed  above.  The  stratification  length  can  be 
determined  through  Equations  37  to  41,  given  the  air-sea  temperature  and  hu¬ 
midity  differences. 

L52.  This  set  of  model  equations,  which  is  considered  here  to  be  sub¬ 
ject  to  refinement,  is  but  one  of  a  suite  of  bulk  coefficient  equation  sets 
proposed  by  various  investigators  who  based  their  formulae  either  on  physical 
considerations  similar  to  those  used  here  or  on  pure  correlations  of  stress 
observations  with  wind  speed  and  density  measurements.  Stress  estimates  from 
various  schemes  vary  widely.  For  example,  Blanc  (1985)  conducted  a  simple 
intercomparison  of  10  stress-estimating  schemes  proposed  within  the  last 
15  years,  using  as  input  over  2,000  mean  wind,  temperature,  and  humidity  ob¬ 
servations  from  the  open  North  Atlantic  Ocean.  Typical  mean  scheme-to-scheme 
variation  was  in  the  range  of  15  to  25  percent,  which  is  reasonable.  However, 
maximum  scheme-to-scheme  variation  was  between  40  and  80  percent  which  is 
somewhat  excessive.  Since  no  direct  stress  measurements  were  made,  it  could 
not  be  discerned  which,  if  any,  of  the  models  was  correct.  Since  the  same 
data  were  used  in  all  compaiisons,  the  differences  had  to  be  a  result  of  the 
model  equations.  This  suggests  that  reliable  stress  estimation  by  this  method 
requires  refinements  of  some  or  all  of  the  models,  combined  with  rigorous 
testing  to  distinguish  the  best.  This  is  especially  true  in  coastal  regions 
where  effects  of  finite  depth  and  cross-shore  nonuniformity  are  not  fully 
understood. 

153.  Bulk  formulations  are  somewhat  sensitive  to  the  accuracy  of  data 
used  in  them.  In  the  simplest  formulae,  this  relates  to  measurements  of  wind 
speed,  air  temperature,  water  temperature,  and  air  humidity  at  a  reference 
level  z^_  .  Accuracy  is  degraded  by  inadequate  frequency  response,  insuffi¬ 
cient  averaging  time,  poor  instrument  calibration,  platform  motion,  and  dis¬ 
turbance  of  the  signal  by  the  platform  itself. 

154.  Rough  estimates  of  length  and  time  scales  for  the  resolution  of 
various  attributes  of  a  turbulent  flow  can  be  deduced  from  reported  measure¬ 
ments.  Errors  due  to  a  known  range  of  calibration  uncertainty  can  be  deduced 
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from  the  model  equations.  Platform  motion  can  bias  results  due  to  spurious 
apparent  signals  and  due  to  motion  of  sensors  through  a  field  with  curvature 
in  the  measured  entity.  These  problems  must  be  resolved  by  analysis  of  indi¬ 
vidual  platforms.  Of  course,  ideal  platforms  for  Eulerian  observations  do  not 
move  relative  to  a  fixed  coordinate  system.  Signal  distortion  due  to  the 
presence  of  a  platform  in  a  flow  field  must  also  be  considered  by  analvsis  of 
individual  platforms,  if  corrections  are  to  be  made.  However,  rough  scaling 
arguments  can  be  employed  to  locate  instruments  so  as  to  reduce  substantially 
this  effect. 

155.  In  regard  to  sampling  length  and  time  scales,  the  hypothesis  gen¬ 

erally  employed  is  that  dominant  turbulent  eddies  in  the  surface  layer  scale 
with  distance  from  the  boundary.  If  a  wave  number  spectrum  of  any  turbulent 

property  is  measured,  it  is  expected  that  the  wave  number  axis  can  be  nor¬ 

malized  by  the  reference  elevation  z  to  obtain  similarity.  Given  the  wave 
number  k  =  2tt/A  ,  where  A  is  the  wavelength,  then  spectra  should  be  similar 
when  plotted  as  a  funtion  of  z^/A  .  Most  measurements  are  made  in  the  fre¬ 
quency  domain,  however.  The  relation  between  frequency  and  wave  number  is 
generally  assumed  to  follow  Taylor's  hypothesis  (Tennekes  and  Lumley  1972). 
This  asserts  that  turbulence  properties  do  not  change  significantly  when  ad- 
vected  past  a  fixed  point  at  reference  wind  speed  U  .  The  cyclic  frequency 

f  corresponding  to  A  is  then  f  =  U^/A  .  Using  this,  the  quantity  z^/A 
becomes  fz^/U^  ,  which  is  sometimes  called  the  natural  frequency. 

156.  Spectra  of  turbulence  properties  have  been  measured  by  several 

groups  of  investigators.  Among  recent  papers  on  this  subject  are  those  by 
Schmitt,  Friehe,  and  Gibson  (1979),  Large  and  Pond  (1981),  and  Smith  and 

Anderson  (1984).  These  measurements  indicate  that  the  bulk  of  turbulent 

fluctuations  occur  for  fz^/U^  £  0.01  .  The  time  scale  corresponding  to  this 
frequency  is  1 / f  -  100  z^/U^  .  For  bulk  measurements  using  mean  properties, 
it  is  only  necessary  to  have  an  instrument  that  responds  to  frequencies  or  the 
order  of  0.01  U  /z  •  A  reasonably  stable  estimate  of  the  mean  is  obtained 
by  averaging  over  about  10  time  periods  corresponding  to  this  frequency,  i.e., 
the  averaging  time  T  is  T  =  1,000  z  /U 

157.  For  typical  values  of  wind  speed  and  reference  elevation  this 

gives  characteristic  averaging  times  necessary  to  resolve  mean  values.  For 

example,  if  z  ■  10  m  and  U  »  10  m/sec,  then  T  »  1,000  sec  *  16.7  min. 
r  r  3 

For  smaller  Uf  ,  a  longer  time  is  required.  For  example,  with  U  =  1  m/sec, 
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T  needs  to  be  10,000  sec  =  2.8  hr.  Since  common  averaging  times  seldom 
exceed  I  hr,  this  suggests  that  a  1-m/sec  wind  is  typically  undersampled  and 
so  may  have  substantial  variability.  In  much  of  the  older  data  used  by 
Kitaigorodskii  (1973)  and  Garratt  (1977),  averaging  times  were  around  10  to 
15  min.  By  the  scaling  given  here,  this  suggests  that  wind  speeds  less  than 
11  to  16  m/sec  were  undersampled  (for  z  =  10  m)  and  so  may  account  for  some 
of  the  reported  data  scatter.  For  wind  speeds  on  the  order  of  5  m/sec,  at 
which  increased  wave  generation  occurs,  the  averaging  time  is  about  30  min. 
Note  that  as  z^_  increases,  so  does  the  averaging  time.  For  z^.  =  20  m,  all 
averaging  times  double  over  those  required  at  z^_  =  10  m. 

158.  Uncertainties  of  stress  estimates  due  to  instrument  error  and 
platform-blockage  effects  were  considered  by  Blanc  (1986a).  In  this  study  he 
compared  four  bulk  schemes  using  weather-ship  data  from  the  open  North 
Atlantic  Ocean  and  allowed  the  data  to  vary  randomly  about  its  reported  values 
by  amounts  described  in  the  literature  as  being  characteristic  of  calibration 
errors  and  ship-influence  errors.  Error  estimates  used  were  fixed  values  and 
not  percentages,  so  relative  errors  tended  to  increase  at  smaller  values  of 
the  primary  variables.  For  sensor  calibration  error,  Blanc  (1986a)  assumed 
wind-speed  errors  of  ±0.5  m/sec  for  <  20  m/sec  and  ±1.0  m/sec  for 

U  >  20  m/sec,  errors  in  air  temperature  and  wet/dry  bulk  temperatures  (for 
humidity  determination)  of  ±0.3°  C  and  errors  in  sea-surface  temperature  of 
±0.5°  C.  Based  on  random  variations  within  these  limits,  Blanc  found  that  all 
four  of  his  test  models  agreed  very  well  with  each  other  and  were  all  within 
10  to  15  percent  of  the  error-free  estimates  of  wind  stress  for 
u^  >  0.25  m/sec.  This  is  a  tolerable  error  for  modeling  purposes. 

159.  For  ship-induced  errors,  Blanc  (1986a)  estimated  variations  in 
wind  speed  due  to  flow  distortion  around  the  ship,  temperature  errors  due  to 
variable  solar  heating  near  and  on  the  ship,  variation  in  humidity  due  to 
ship-induced  spray,  and  variation  in  sea-surface  temperature  due  to  ship 
engine  cooling  water  discharge.  These  variations  were  of  the  same  order  or 
slightly  larger  than  the  calibration  errors  discussed  above.  Ship-induced 
errors  resulted  in  25-  to  35-percent  variations  in  stress  estimates  relative 
to  those  from  uncontaminated  data.  Again,  this  result  is  marginally  accept¬ 
able  for  modeling  purposes. 

160.  In  the  study  by  Blanc  (1986a),  the  wind  blockage  effect  was  esti¬ 
mated  to  be  on  the  order  of  10  percent  of  mean  wind  speed.  This  alone  would 
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give  variations  on  the  order  of  20  percent  in  wind  stress.  This  estimate  of 
error  in  wind  speed  may  be  conservative  for  some  ship-mounted  anemometers, 
however.  In  a  subsequent  publication,  Blanc  11986b)  reported  results  of  wind- 
tunnel  tests  to  determine  wind  speed  and  direction  deviations  at  standard 
anemometer  locations  on  a  US  Navy  vessel.  He  found  wind-speed  errors  ranging 
from  40  percent  low  to  50  percent  high,  depending  on  direction  of  wind  ap¬ 
proach  to  the  ship.  If  used  in  a  quadratic  drag  law,  these  results  would  give 
errors  in  stress  ranging  from  60  percent  low  to  over  100  percent  high.  This 
is  not  acceptable  for  use  in  an  otherwise  accurate  oceanic  model. 

161.  The  studies  by  Blanc  (1986a, b)  are  mentioned  because  they  empha¬ 
size  chat  poorly  sampled  mean  data  will  vield  poor  results  even  if  the  drag- 
and  mass-transfer  coefficients  are  perfectly  known.  In  application  of  the 
bulk  method,  careful  consideration  of  sampling  scheme,  calibration  method,  and 
platform  effects  are  required  to  ensure  results  of  reasonable  quality.  The 
same  considerations  are  even  more  critical  in  experiments  designed  to  help 
refine  drag  coefficient  models.  The  purpose  of  this  report  is  to  justify  and 
describe  such  an  experiment.  Measurement  requirements  are  discussed  in  the 
remainder  cf  Part  IV,  which  deals  with  alternate  methods  of  wind-stress  mea¬ 
surement,  and  in  Part  V  which  describes  an  experiment  scheme. 

Profile  Method 


162.  In  the  profile  method,  measurements  of  mean  horizontal  wind  speed, 

mean  temperature,  and  mean  humidity  are  made  within  the  atmospheric  surface 

layer  (but  above  the  sublayer)  at  a  minimum  of  three  elevations.  The  data  are 

then  fitted  to  the  set  of  Equations  34  through  41  to  find  the  three  basic  flow 

parameters  u.  ,  z  ,  and  L  .  Wind  stress  is  then  determined  from 
*  o 


163.  This  method  is  a  slightly  more  objective  application  of  the  bulk 
method.  In  the  bulk  method,  the  flow  profile  relationships  given  by  Equa¬ 
tions  29  to  33  are  integrated  from  the  surface  (z  =  z  )  to  some  reference 

o 

elevation  (z  =  z^_) .  In  the  profile  method,  the  same  equations  are  integrated 
between  levels  in  the  surface  layer  only.  This  has  the  advantages  that  a 

model  for  z  is  not  necessary,  sea-surface  temperature  does  not  have  to  be 

o 

measured,  and  sca-surface  humidity  does  not  have  to  be  assumed  to  exist  at 
100-percent  saturation.  Requirement  ■  for  sensor  accuracy  and  frequency 
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response  are  the  same  as  for  instruments  used  in  the  bulk  method. 

164.  The  disadvantages  of  the  profile  method  are  that  it  requires  the 
same  assumptions  of  steadiness  and  horizontal  uniformity  as  are  required  for 
the  bulk  method,  more  sensors  are  required,  and  the  platform  used  must  be  un¬ 
obstructive  at  all  measurement  levels.  Furthermore,  for  optimum  accuracy,  the 
measurement  levels  should  be  spaced  logarithmically  in  z  due  to  the  dominant 
logarithmic  nature  of  the  profiles  given  by  Equations  34,  35,  and  36.  Wide 
vertical  separation  of  sensors  at  upper  levels  is  necessary.  This  places 
added  structural  constraints  on  the  platform  and  runs  the  risk  that  upper  sen¬ 
sors  are  outside  the  surface  layer  during  times  of  low  wind  speeds  or  high 
stable  stratification. 

165.  In  the  bulk  method  only  four  sensors  are  required  (sea-surface 
temperature  along  with  mean  wind  speed,  air  temperature,  and  humidity  at  the 
single  reference  elevation).  In  the  profile  method,  at  least  nine  sensors  are 
required  (mean  wind  speed,  air  temperature,  and  humidity  at  a  minimum  of  three 
elevations)  so  there  are  added  logistical  burdens  of  obtaining  and  maintaining 
a  greater  number  of  sensors. 

166.  Blanc  (1983),  in  an  analysis  of  the  mechanical,  logiscical,  and 
climatological  constraints  of  atmospheric  measurements  in  the  marine  environ¬ 
ment,  suggested  that  the  profile  method  is  best  suited  for  flux  estimation 
over  long  durations  and  recommended  numbers  and  spacings  of  instruments  to 
accomplish  this.  The  profile  method  has  been  employed  over  water  by  Hsu 
(1972),  SethuRaman  and  Raynor  (1975),  and  Frenzen  and  Hart  (1975). 

167.  However,  no  experiment  over  water  has  been  conducted  wherein  the 
flux-profile  relationships  used  in  any  formulation  are  tested  for  validity. 
This  is  extremely  difficult  in  the  open  ocean  due  to  mechanical  constraints  of 
sea-gcing  platforms.  Near  land,  such  an  experiment  is  more  feasible  but  the 
higher  probability  (in  general)  of  horizontal  nonuniformity  would  make  the 
profile  method  as  suspect  as  the  bulk  method.  Since  the  bulk  method  is  the 
easiest  to  employ  in  large-scale  climatological  studies,  this  suggests  that 
experiments  testing  flux-bulk  formulations  may  be  preferable.  The  correctness 
of  the  profile  relationships  would  then  be  implicit  in  the  results  of  such 
experiments . 
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Direct  Method 


168.  In  the  direct  method,  the  total  (three  component)  velocity  field 
is  measured  at  a  point  at  frequencies  high  enough  to  resolve  the  fluctuating 
horizontal  and  vertical  velocities  in  the  bandwidth  wherein  the  correlation 
u'w'  ,  and,  by  Equation  3,  the  stress  can  be  computed  directly.  This  method 
requires  the  fewest  assumptions  about  the  flow  field  because  it  is  a  direct 
estimator  of  the  vertical  Reynolds  momentum  flux  at  the  level  of  measurement. 
The  only  assumption  is  that  a  direct  measurement  of  u'w'  at  some  elevation 
above  the  sea  surface  represents  the  shearing  stress  at  the  sea  surface.  If 
this  assumption  is  violated  due  to,  say,  strong  horizontal  nonuniformity,  then 
a  more  complete  treatment  of  the  problem  is  necessary.  In  such  a  treatment  it 
is  still  necessary  to  know  the  turbulent  flux  of  momentum  through  a  control 
surface  parallel  to  the  mean  sea  surface.  A  direct  measurement  of  this  flux 
is  the  only  reliable  estimator  in  this  case,  so  the  need  remains  for  accurate 
sensing . 

169.  Devices  to  measure  all  three  components  of  velocity  at  frequencies 
necessary  to  resolve  the  Reynolds  stresses  are  inherently  more  sensitive  than 
the  typically  robust  instruments  used  in  the  bulk  and  profile  methods.  This 
means  they  are  more  susceptible  to  calibration  error,  platform  problems,  and 
environmental  degradation.  Futhermore,  a  greater  volume  of  data  needs  to  be 
accumulated  to  compute  the  cross  correlations.  These  conditions  increase  sig¬ 
nificantly  the  expense  of  obtaining,  deploying,  and  maintaining  such  instru¬ 
ments.  As  a  result,  experiments  using  such  instrumentation  tend  to  be  of 
short  duration  (usually  less  than  I  month)  and  total  historic  data  are  some¬ 
what  sparse.  Blanc  (1985)  estimated  that  less  than  400  hr  of  direct  measure¬ 
ments  over  water  had  been  reported  in  the  literature  at  the  time  of  his 
search.  Any  additions  to  this  data  base  are  of  great  value. 

170.  In  making  direct  measurements,  the  averaging  time  is  the  same  as 
for  the  bulk  and  profile  methods.  The  necessary  frequency  response  of  the 
sensors  is  higher.  Measurements  of  cospectra  of  the  product  u'w'  by 
Schmitt,  Friehe,  and  Gibson  (1979),  Large  and  Pond  (1981),  and  Smith  and 
Anderson  (1984)  indicate  that  the  bulk  of  the  stress  fluctuations  exist  for 
natural  frequency  fz^/U^  <,  10  or  f  <,  10  U  / z^_  .  If  this  represents  the 
Nyquist  frequency  of  spectral  resolution,  then  the  sampling  frequency  must  be 
twice  this  or  f  -  20  U  /z^  .  Sampling  requirements  are  higher  for  high 
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20  Hz 


reference  wind  speeds  and  low  measurement  elevations.  For  example,  f  - 

for  U  -  10  m/sec  with  z  -  10  m  and  f  -  10  Hz  for  0  =  10  m/sec  with 
r  r  r 

z  -  20  m  but  f  -  40  Hz  for  U  -  20  m/sec  and  z  -  10  m. 
r  r  r 

171.  Devices  which  can  resolve  these  frequencies  are  tvpicallv  heat 
transfer  (hot  films  and  hot  wires)  or  acoustic  devices.  They  tend  to  be  frag¬ 
ile,  expensive,  and  short  lived  in  the  field.  More  robust  and  less  expensive 
mechanical  devices  (wind  vanes  and  cup  or  impellor  anemometers)  tend  to  have  a 
lower  frequency  response  than  that  required  above.  They  can  be  used  for 
direct  stress  measurements  if  their  signals  are  corrected  through  a  measured 
transfer  function  for  decreased  high  frequency  response.  Finkelstein  et  al . 
(1986)  made  intercomparison  measurements  between  a  sonic  anemometer  (with  good 
high-frequency  response)  and  six  commonly  employed  configurations  of  lower 
frequency  response  mechanical  devices.  Their  results  suggest  that  transfer 
functions  can  be  found  to  improve  frequency  response  up  to  fz^/U^  -  1  .  This 
is  still  a  decade  lower  than  desired.  However,  a  further  correction  can  be 
employed  based  on  the  observed  u'w'  cospectra  mentioned  above.  For  a  given 
stability,  those  spectra  appear  to  be  similar  in  shape  so  that  the  missing 
fractional  variance  (on  the  order  of  10  to  30  percent)  can  be  estimated  from 
the  part  of  the  signal  that  is  resolved  with  the  low-frequency  devices.  This 
procedure  allows  longer  term  measurements  of  reasonable  accuracy  to  be  per¬ 
formed  with  less  delicate  instrumentation. 

172.  A  second  consideration  in  the  direct  method  is  leveling  of  the 
sensors  used  to  measure  the  three  components  of  velocity.  The  reason  for  this 
is  that  u'w'  is  but  one  element  of  a  nine-element  correlation  tensor.  Small 
misalignments  effectively  rotate  the  tensor  so  the  desired  quantity  is  con¬ 
taminated  by  contributions  from  other  terms.  The  essential  argument  is  that 
if  fluctuating  u'  and  w'  are  measured  in  a  coordinate  system  that  is  out 

of  level  bv  an  angle  X  ,  then  the  measured  velocities  u'  and  w'  are  re¬ 
in  m 

lated  to  the  true  velocities  u'  and  w'  by  u'  =  u'  cos  X  -  w'  sin  X  and 

m  m 

w'  =  u'  sin  X  +  w'  c^s  X  ,  where  X  is  the  angle  in  the  vertical  plane  from 

iii  TTl 

the  true  x-axis  to  the  x-axis  of  the  measuring  system.  Note  that  here  the 

coordinate  system  is  aligned  with  the  x-axis  in  the  mean  wind  direction  so  the 

cross-wind  fluctuating  velocity  does  not  enter  this  argument.  It  should  be 

included  in  the  general  case.  The  correlation  of  u'  and  w'  is  simply  the 

average  of  the  product  of  the  two  expressions.  With  some  manipulation  this 

vields  u'w'  =  (u'u'  -  w'w')  cos  X  sin  X  +  u'w'  (cos  X  -  sin“  X)  .  This 
m  tn  Tn  m  mm 
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expression  can  be  approximated  in  terms  of  the  true  friction  velocity 

u.  =  -u'w'  and  the  measured  friction  velocity  u.  =  -u' w'  .  Observations  of 
*  *m  m  m 

variances  of  fluctuating  velocities  (Hinze  1959)  in  the  surface  layer  indicate 

that  w'w'  -  u-  and  u'u'  -  5  u  ,  to  within  50  percent  or  so.  Using  these 
m  in  m  m  z..  ° 

approximations  in  the  expression  for  uTwT  shows  the  true  friction  velocity 

2  -?  2  2 

squared  to  be  u“  -  u'  (cos“  X  -  sin  X  -  5  sin  X  c os  x)  .  If  measured 
values  are  assumed  to  be  relative  to  a  level  coordinate  system,  the  error  is 
the  term  in  parentheses.  Table  4  shows  this  term  as  a  function  of  X  . 

Table  4 

Estimated  Error  in  Wind  Stress  Due 
to  Measurement  System  Being 
Out  of  Level  by  Angle  X 


X  ,  deg 

2,  2 

U*/U*m 

-3 

1.26 

-2 

1.17 

-1 

1.09 

0 

1.00 

1 

0.91 

2 

0.82 

3 

0.73 

173.  According  to  Table  4,  the  measurement  system  must  be  level  to 
within  about  ±2  deg  for  accuracy  to  within  ±20  percent.  This  is  extremely 
difficult  on  a  floating  platform  such  as  a  ship  or  buoy.  It  requires  that 
direction  estimates  be  very  accurate  and  that  some  type  of  inertial  displace¬ 
ment  measuring  system  be  employed  to  correct  for  all  platform  motion.  This 
effect  and  flow  distortion  by  large  platforms  (such  as  ships  and  offshore  oil 
platforms)  give  cause  why  little  direct  wind-stress  data  have  been  obtained  in 
the  open  ocean.  These  problems  are  more  readily  controlled  on  a  fixed  tower 
of  small  cross  section.  Such  towers  are  used  in  shallow  and  nearshore  waters 
where  stays  and  bracing  can  be  fixed  to  the  seafloor.  Geernaert,  Katsaros  and 
Richter  (1986)  used  such  a  tower  in  16  m  of  water  to  deploy  a  three-axis 
impellor  anemometer.  Smith  and  Anderson  (1984)  obtained  very  good  results 
using  a  tower  placed  near  the  waterline  on  a  beach  and  upon  which  was  mounted 
an  acoustic  anemometer. 
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174.  When  Che  direct  method  is  used  in  stratified  flow,  measurements 


must  also  be  made  of  the  correlations  w'T'  and  w' q '  ,  which  are  the  turbu¬ 
lent  fluxes  of  sensible  heat  and  water  vapor,  respectively.  Fluctuating  tem¬ 
peratures  are  usually  made  with  temperature-sensitive  resistance  elements  of 
low  thermal  inertia  in  a  bridge  circuit.  Humidity  fluctuations  are  often 
quantified  through  the  absorption  of  light  at  a  fixed  wavelength  by  water 
vapor  in  air.  Variations  in  light  intensity  are  related  in  a  known  way  to 
water  vapor  concentration.  The  required  frequency  bandwidth  for  measurement 
of  these  variables  is  the  same  as  for  velocity  fluctuations. 

175.  To  obtain  the  necessary  frequency  response  for  temperature  and 
humidity  sensors,  the  sampling  volumes  are  on  the  order  of  centimetres  in 
size.  This  makes  such  instruments  somewhat  fragile.  In  addition,  problems 
have  been  noted  in  field  use  of  these  instruments  due  to  the  relatively  high 
concentrations  of  salt  in  the  marine  atmospheric  surface  layer.  Phelps  and 
Pond  (1971)  noted  curious  differences  between  temperature  spectra  and  humidity 
spectra  in  measurements  in  the  Pacific  Ocean  off  San  Diego,  California.  Such 
differences  were  not  expected  in  Monin-Obukhov  similarity  theory.  The  differ¬ 
ences  were  explained  by  Schmitt,  Friehe,  and  Gibson  (1978)  as  being  due  to 
saltwater  accumulation  on  temperature  sensors.  Condensation  and  evaporation 
of  water  on  the  salt  droplets  in  conditions  of  high  humidity  yielded  a  spuri¬ 
ous  temperature  signal  because  of  latent  heat  released  or  absorbed.  As  noted 
by  Blanc  (1983),  salt  accumulation  can  also  alter  the  optical  transmissivity 
or  cause  corrosion  of  windows  in  light-intensity  sensitive  devices  such  as 
those  used  to  detect  humidity. 

176.  These  problems  indicate  that  collecting  information  from  these 
devices  is  labor-intensive,  requiring  virtually  constant  attention  and  fre¬ 
quent  cleaning  and  exchange  of  sensing  elements.  On  the  other  hand,  knowledge 
of  these  problems  aids  in  proper  experiment  design.  Furthermore,  the  added 
expense  is  compensated  for  by  the  acquisition  of  direct  Reynolds  flux  measure¬ 
ments,  which  is  the  only  proper  way  to  calibrate  or  test  model  equations. 

177.  Another  constraint  on  experiment  design  is  in  the  spacing  of  sen¬ 
sors.  The  natural  frequency  fz^/U^  is  the  Taylor  transform  of  wave  number 
dependency  z^/A  .  For  fz^/U^  *  10  ,  the  corresponding  high-frequency  wave¬ 
length  is  found  from  z  /A  *  10  or  A  *  z  / 1 0  .  For  z  =  10  m,  this  means 

r  r  r 

that  correlations  are  needed  for  spacings  on  the  order  of  A  =  1  m.  To  cor¬ 
relate  temperature  or  humidity  with  vertical  velocity,  the  respective  sensors 
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should  reside  in  a  volume  of  length  scale  1  m  or  less.  For  =  20  m,  the 
spacing  should  be  less  than  2  m. 


Inertial  Subrange  Method 

178.  A  fourth  method  of  estimating  turbulent  fluxes  is  related  to  the 
functional  shape  of  turbulence  spectra  in  a  particular  band  of  wave  numbers. 

As  with  the  bulk  and  profile  method,  it  is  indirect  and  relies  on  the  assump¬ 
tions  of  steadiness  and  horizontal  uniformity. 

179.  The  method  is  derived  from  Kolmogorov’s  hypothesis  (Tennekes  and 
Lumlev  1972)  which  argues  that,  for  large  Reynolds  number  flows,  there  exists 
a  band  in  the  wave  number  spectrum  of  turbulent  kinetic  energy  the  structure 
of  which  depends  only  on  wave  number  and  turbulent  energy  dissipation  rate. 
That  is,  if  S(k)  is  velocity  variance  per  wave  number  increment  at  wave  num¬ 
ber  k  and  e  is  the  rate  of  dissipation  from  Equation  20,  with  units  of 
velocity  variance  per  time,  then  a  closed  set  of  variables  is  S(k)  ,  k  ,  and 
£  .  Dimensional  analysis  dictates  the  relationship  as 


S (k)  =  Ae 


2/3k-5/2 


(45) 


where  A  is  the  universal  constant,  often  called  Kolmogorov's  constant. 

180,  If  the  flow  is  unstratified,  then  £  in  Equation  45  can  be  re¬ 
placed  by  Equation  28  which  yields 


,2/3 


S(k) 


-5/3 


\KZ 


(46) 


The  next  step  is  to  convert  Equation  46  from  wave  number  dependence  to  fre¬ 
quency  dependence  through  Taylor's  hypothesis.  The  relation  between  cyclic 
frequency  f  (in  Hertz)  and  radian  wave  number  k  in  a  Taylor  transformation 
is  k  =  2 n f / 1'  where  U  is  local  mean  wind  speed.  Since  S(k)  is  velocity 
variance  per  wave  number  increment,  it  is  converted  to  frequency  dependence 
by  S(f)  =  U  S(k)/2n  .  Using  these  transformations  in  Equation  46  and  combin¬ 
ing  terms  gives 
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S(f) 


(47) 


18L.  The  range  of  frequencies  for  which  Equation  47  is  valid  is  called 
the  inertial  subrange.  Since  there  is  no  direct  dependence  on  production  of 
turbulence  in  this  high-frequency  range,  it  is  generally  assumed  that  the 
motions  follow  the  theory  of  isotropic  turbulence.  In  this  case,  the  spectrum 
of  u'  follows  Equation  47  and  the  spectra  o f  v'  and  w'  (the  cross-wind 
components)  are  each  four-thirds  times  this.  Thus,  if  A  is  the  Kolmogorov 
constant  for  the  downwind  spectrum,  then  the  coefficient  of  Equation  47  for 
the  spectrum  of  horizontal  winds  should  be  7/3A  and  the  coefficient  for  the 
spectrum  of  winds  in  all  directions  should  be  11/3A  .  Frenzen  (1983)  lists 
some  experimentally  determined  values  of  A  ,  which  cover  the  range  0.50  <  A 
<  0.59  ,  and  recommends  A  =  0.52  from  his  own  work.  Schmitt,,  Friehe,  and 
Gibson  (1979)  use  A  =  0.50  in  their  work.  Large  and  Pond  (1981,  1982)  and 
Smith  and  Anderson  (1984)  use  A  =  0.55  .  This  last  value  is  a  reasonable 
estimator  and  can  be  assigned  an  uncertainty  of  ±10  percent. 

2 

182.  The  way  in  which  Equation  47  is  used  to  estimate  u*  is  to  mea¬ 
sure  and  compute  the  one-,  two-,  or  three-dimensional  spectrum  of  fluctuating 
velocities  and  then  fit  Equation  47  to  these  measurements  by  varying  u^  . 

This  is  subject  to  two  constraints.  The  mean  flow  Reynolds  number  (Uz/v)  must 
be  large  enough  for  an  inertial  subrange  to  exist  and  the  curve-fitting  pro¬ 
cess  must  only  be  done  in  that  subrange  of  frequencies.  Tennekes  and  Lumlev 

(1972)  estimate  that  an  inertial  subrange  should  exist  for  Uz/v  >  10^  . 

-5  2 

Since  v  =  1.5  x  10  m  /sec  for  air,  this  constraint  is  satisfied  for  mean 
wind  speeds  as  low  as  1  m/sec  at  elevations  as  low  as  1  m.  Higher  wind  speeds 
and  higher  observation  elevations  serve  to  increase  the  Reynolds  number  so 
this  constraint  is  not  a  problem  in  most  flows  of  interest. 

183.  Spectral  measurements  over  water  by  Schmitt,  Friehe,  and  Gibson 
(1979),  Large  and  Pond  (1981)  and  Smith  and  Anderson  (1984)  indicate  that  the 
inertial  subrange  has  a  low  frequency  bound  at  a  natural  frequency  of  about 
fz/U  =  1  ,  Since  the  upper  frequency  bound  in  the  direct  method  of  stress 
estimation  was  fntmH  no  be  fz/U  «  10  ,  this  means  that  sensors  of  the  same 
type  can  be  used  in  the  inertial  subrange  method.  That  is,  if  the  spectrum 
is  resolved  for  the  range  1  <  fz/U  <  10  ,  then  Equation  47  can  be  used  to 
find  u*  . 
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184.  The  advantages  of  this  method  are  several.  First,  only  a  one¬ 
dimensional  spectrum  is  necessary  to  estimate  stress  so  it  is  not  necessary  to 
measure  all  three  velocity  components.  Second,  stress  can  be  estimated  from 
spectra  of  just  horizontal  velocities,  as  could  be  obtained  from  a  cup 
anemometer  with  a  high  enough  frequency  response.  Frenzen  (1983)  describes 
results  from  measurements  using  such  a  device.  Third,  because  the  range  of 
frequencies  in  the  inertial  subrange  is  relatively  high,  measurements  are  less 
sensitive  to  platform  motion  than  are  measurements  using  the  direct  method. 
This  suggests  that  such  devices  can  be  used  on  shipboard  if  adequate  care  is 
taken  to  avoid  flow  distortion  by  the  ship's  structure.  Since  sensors  which 
work  at  the  frequencies  cited  tend  to  be  fragile,  long  term,  unattended 
deployment  on  buoys  or  remote  stations  is  not  recommended. 

185.  If  the  flow  is  stratified,  then  the  inertial  subrange  method  must 
be  adapted  to  account  for  the  fluxes  of  temperature  and  humidity.  The  conven¬ 
tional  derivation  is  given  here  for  temperature,  following  Tennekes  and  Lumley 
(19721.  The  arguments  for  humidity  are  identical.  When  Kolmogorov's  hypothe¬ 
sis  is  applied  to  a  state  variable  (temperature  or  humidity)  it  is  assumed 
that  the  variable  acts  as  a  passive  contaminant  at  wave  numbers  in  the  iner¬ 
tial  subrange.  That  is,  density  variations  have  a  primary  effect  at  large 
(production)  length  scales  but  not  directly  at  shorter  length  scales.  It  is 
then  assumed  that  the  temperature  variance  spectral  density  ST(k)  is  a  func¬ 
tion  of  wave  number  k  and  the  molecular  loss  rate  of  temperature  variance 

E.J,  .  One  more  variable  is  needed  because  has  dimensions  of  squared  tem¬ 

perature  per  time  and  time  is  not  a  dimension  of  k  or  S^(k)  .  The  added 
variable  is  assumed  to  be  e  ,  the  energy  dissipation  rate. 

186.  The  set  of  variables  ST(k)  ,  k  ,  eT  ,  and  e  forms  a  dimen¬ 
sionally  unique  grouping  which  can  be  written  as 

ST(k)  -  A^"173  et  k"5/3  (48) 

where  A^  is  a  universal  constant.  The  same  argument  applied  to  humidity 
yields 


S  (k) 


-5/3 


q 


(49) 
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where  is  a  universal  constant.  Equations  45,  48,  and  49  are  used  to  find 

the  fluxes  once  e  ,  e_  ,  and  z  are  modeled. 

187 .  In  the  case  of  stratified  flow,  an  estimate  of  z  is  obtained 
from  Equation  27  under  the  assumption  that  D  =  0  by  using  the  Monin-Obukhov 
length  L  as  notation  in  the  second  term  on  the  right  side  and  then  replacing 
the  term  on  the  left  side  with  Equation  29.  The  expression  for  z  is  then 

3 

u* 

e  =  7? 

This  expression  is  subject  to  argument  because  it  has  not  been  shown  that  D 
is  negligible  in  stratified  flow.  Frenzen  (1983)  indicates  that  D  mav  be  on 
the  order  of  20  to  40  percent  of  the  dissipation  rate.  Long  (1981)  argues 
that  D  is  on  the  order  of  4  times  the  buoyant  sink  term  in  Equation  27  for 
stratified  flow  which  is  stable.  This  is  a  matter  for  future  research.  Equa¬ 
tion  50  has  been  used  in  several  recent  investigations  (Schmitt,  Friehe,  and 
Gibson  1979;  Large  and  Pond  1981;  Smith  and  Anderson  1984;  Boyle,  Davidson, 
and  Spiel  1987)  and  so  is  worth  consideration  as  an  interim  model. 

188.  The  dissipation  rate  for  temperature  variance  is  estimated  from 
the  temperature-variance  equation  (Tennekes  and  Lumley  1972).  For  steady, 
uniform  flow  and  by  assuming  negligible  gradient  diffusion  of  temperature 
variance  (like  assuming  D  =  0  in  the  energy  equation) ,  this  is 


at)  - 1 


f50) 


w'T'  — 
3z 


(51) 


When  combined  with  Equation  30,  which  was  deduced  from  Monin-Obukhov  similar¬ 
ity,  the  temperature  gradient  3T/3z  can  be  eliminated.  This  yields 


eT 


(w’T’)2 


The  equivalent  expression  for  humidity  is 


(52) 
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n 


(53) 


Using  Equations  50  and  52  in  Equation  48,  converting  S^Ck)  and  k  to 
(l72n)S^,(f )  and  (2'r/L’)f  ,  respectively,  by  the  Taylor  transformation  and 
rearranging  terms  gives 


(54) 


The  expression  for  the  spectrum  of  humidity  is 


(55) 


The  one-dimensional  velocity  variance  spectrum  is  obtained  by  using  Equa¬ 
tion  50  in  Equation  45,  transforming  Taylor's  hypothesis,  and  rearranging 
terms.  This  yields 

»<»  -  *4  [uD  - 

189.  Flux  estimates  are  obtained  through  spectral  measurement  of  fluc¬ 
tuating  velocity,  temperature,  and  humidity  and  then  fitting  Equations  54,  55, 

and  56  to  these  measurements  using  Equations  32  and  33  for  <f>  ,  f>  ,  and  t 

m  q  T 

along  with  the  definition  of  L  following  Equation  28.  Values  for  mean  wind 
velocity  U  and  measurement  elevation  z  are  also  needed.  Large  and  Pond 
(1982)  and  Smith  and  Anderson  (1984)  use  =  A  =  0.80  for  the  Kolmogorov 
constants,  citing  an  uncertainty  on  the  order  of  10  percent  based  on  previous 
work.  Where  wind  stress  is  the  primary  variable  sought,  uncertainty  in  these 
constants  is  not  too  harmful  since  the  mass  fluxes  are  determined  only  for 
corrections  to  stress  estimates.  The  equation  set  is  transcendental  in  z/L 
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and  so  the  curve-fitting  process  must  be  done  in  an  iterative  fashion.  This 
can  be  done  readily  on  any  common  computer. 

L90.  Tests  of  the  hypotheses  made  in  the  derivation  of  this  method  have 

been  carried  out  by  Schmitt,  Friehe,  and  Gibson  (1979),  Large  and  Pond  (1981, 

1982),  Smith  and  Anderson  (1984),  and  others  in  measurements  over  water. 

Results  were  very  good  where  reasonable  assurance  of  negligible  platform 

blockage,  platform  motion,  nonuniformity,  and  instrument  error  effects  could 

-5/3" 

be  obtained.  All  spectra  had  an  f  dependence  to  within  measurement  ur.- 

certainity.  Spectral  amplitudes  in  the  inertial  subrange  agreed  well  with 
measurements  made  over  land.  A  direct  comparison  of  u*  determined  from  the 
inertial  subrange  method  with  u^  determined  from  the  direct  method  was  made 
by  Large  and  Pond  (1981).  In  the  mean,  the  two  methods  agreed  to  within  about 
3  percent.  The  scatter  was  generally  within  ±20  percent,  which  is  reasonable 
since  the  uncertainty  from  each  method  is  about  ±10  percent  for  u*  . 

191.  An  application  of  the  method  was  reported  by  Boyle,  Davidson,  and 
Spiel  11987)  in  measurements  over  the  open  Pacific  Ocean.  In  steady,  uniform 
conditions  the  inertial  subrange  method  gave  results  that  compared  very  fav¬ 
orably  with  the  bulk  method,  based  on  results  given  by  Large  and  Pond  (1982). 
As  mentioned  previously,  this  experiment  showed  that  the  bulk  method  severely 
underestimated  wind  stress  when  wind  and  wind-driven  waves  were  not  in  equi¬ 
librium.  Since  wind-driven  waves  tend  to  be  long  lived,  relative  to  time 
scales  of  Reynolds  stresses  and  inertial  subrange  variances,  this  suggests 
that  the  inertial  subrange  method  is  a  better  indicator  of  instantaneous  wind 
stress  than  the  bulk  method.  That  is,  corrections  to  or  further  modeling  of 
the  bulk  method  to  account  for  time  dependence,  variations  in  sea  state, 
finite  water  depth,  etc.,  can  be  achieved  using  the  inertial  subrange  method 
to  estimate  wind  stress. 


Summary 


192.  Four  methods  are  described  here  for  direct  or  indirect  methods  of 
wind-stress  measurement.  General  instrumentation  requirements,  averaging 
times,  and  frequency  responses  are  given.  To  effect  an  improved  parametric 
dependence  of  the  bulk  method,  alternate  stress  estimates  need  to  be  made 
simultaneously.  In  order  of  preference,  the  best  method  is  the  direct  method, 
followed  by  the  inertial  subrange  and  profile  methods.  In  order  of  expense  of 
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instrument  acquisition  and  maintenance,  the  list  is  reversed.  In  all  of  the 
methods,  the  theory,  instrument  capability,  and  deployment  constraints  ar» 
being  developed  and  improved  concurrently.  For  this  reason,  any  experiment 
design  should  include  as  many  of  these  methods  as  is  possible  to  intercompare 
results  and  ensure  reliability.  This  is  especially  true  in  shallow  water  and 
nearshore  environments  where  the  physics  of  air-sea-land  interaction  is  more 
complex  than  in  steady,  uniform  conditions  over  land  or  over  the  open  ocean. 
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PART  V:  AN’  EXPERIMENT  PROGRAM 


193.  Ar  experiment  program  has  been  designed  to  test  and  evaluate  some 
of  the  hypotheses  expressed  above  concerning  modification  of  drag  coefficients 
in  shallow  and  nearshore  waters.  The  intent  is  to  measure  over-water  mean 
wind  speed  and  wind  stress,  so  as  to  compute  a  drag  coefficient,  and  a  number 
of  auxiliary  parameters  to  determine  drag  coefficient  dependence  on  them. 

These  parameters  include  wind  direction,  air-sea-land  temperatures,  fluxes  cf 
heat  and  humidity,  water  depth,  distance  to  shore,  and  detailed  sea  state  pa¬ 
rameters  derived  from  currents  and  directional  wave  spectra. 

194.  For  logistical  and  physical  reasons,  the  site  chosen  for  the  ex¬ 
periment  program  was  the  Field  Research  Facility  (FRF)  of  the  l'S  Army  Engineer 
Waterways  Experiment  Station,  Coastal  Engineering  Research  Center.  This  site 
is  on  the  Atlantic  coastal  barrier  islands  of  North  Carolina  just  north  of  the 
village  of  Duck  (Figure  l) .  A  pier  extending  about  500  m  seaward  of  the 
waterline  provides  dry  access  to  an  over-water  location,  ideal  for  nearshore 
wind  stress  measurements.  The  site  is  also  very  uniform  for  several  kilo¬ 
metres  in  the  alongshore  direction.  This  allows  reasonable  correlation  of 
drag  coefficient  with  the  parameters  listed,  without  added  complication  due  to 
a  convoluted  coastline. 

195.  Because  of  the  wide  variety  of  parameters  needed,  this  program  was 
conducted  as  part  of  a  large-scale  joint  experiment  known  as  SUPERDUCK.  It 
was  coordinated  by  the  Corps  of  Engineers  and  had  added  participation  by  uni¬ 
versity  scientists  and  other  government  research  groups.  Overview  papers  de¬ 
scribing  Sl’PERDUCK  are  expected  to  appear  soon  in  the  literature.  This  report 
describes  the  wind-related  measurements  of  this  experiment. 

196.  The  atmospheric  measurement  program  took  place  in  two  parts.  The 
first  part  took  place  in  the  fall  of  1985  in  a  pre-SUPERDUCK  experiment  known 
as  DUCK85.  This  work  was  intended  to  provide  a  preliminary  survey  of  horizon¬ 
tal  variability  of  winds  at  the  FRF  and  a  measure  of  the  adequacy  of  various 
existing  platforms  for  accurate  wind-stress  measurements.  Based  in  part  on 
the  results  of  DUCK85,  the  main  experiment  was  designed  to  make  optimum  use  of 
existing  platforms  yet  minimize  data  contamination  due  to  platform  effects. 
This  report  describes  the  DUCK85  experiment  and  results  in  detail  and  gives  a 
cursory  description  of  the  SUPERDUCK  experiment.  A  detailed  description  of 
the  latter  experiment  and  its  results  will  be  given  in  a  subsequent  report. 
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Figure  I.  Site  locations  of  the  CERC  Field  Research  Facility  (from  Miller  1982) 


DUCK85 


I°~.  The  preliminary  wind  experiment  had  four  main  objectives.  The 
first  was  to  provide  over-water  environmental  wind  conditions  for  other  par¬ 
ticipants  in  the  DUCK85  experiment  (Hubertz  1987) .  The  second  was  to  deter¬ 
mine  the  magnitude  of  cross-shore  horizontal  nonuniformity  of  winds  from 
various  directions  at  the  FRF  site.  The  third  was  to  compare  measurements  by 
the  permanently  installed  anemometer,  used  by  FRF  personnel  for  routine 
weather  observations,  to  measurements  over  water  to  see  how  well  it  repre¬ 
sented  over-water  winds.  The  fourth  objective  was  to  test  the  adequacv  of  a 
tower  at  the  seaward  end  of  the  pier  for  making  wind-stress  measurements  using 
the  profile  method  described  in  Part  IV. 

Instrumentation 

198.  Wind  measurements  were  made  from  three  platforms,  two  on  land  near 
the  landward  end  of  the  pier  and  one  at  the  seaward  end  of  the  pier.  Figure  2 
is  a  schematic  plan  view  of  the  FRF  site,  showing  locations  of  some  of  the 
instrumentation  from  the  DUCK85  experiment.  Figure  3  is  a  schematic  diagram 
showing,  in  three  pieces,  shore-normal  sections  of  the  three  wind  measurement 
platforms . 

199.  The  pier-end  tower  (Figure  2,  Item  12;  and  Figure  3)  is  an  open- 
frame,  three-legged  tower  with  legs  0.61  m  apart,  and  with  a  maximum  elevation 
of  22.1  n  above  mean  sea  level  (msl) .  Its  base  is  on  the  southern  end  of  the 
terminal  bent  of  the  pier  at  6.8  m  above  msl.  The  pier  deck  is  at  7.7  m  above 
msl.  One  of  the  land-based  platforms  is  a  6.2-m  tower  atop  the  FRF  laboratory 
building  (Figure  2,  Item  13;  Figure  3).  The  top  of  this  tower  is  at  18.7  m 
above  msl.  It  is  called  the  Skyvane  platform  in  this  report  because  of  the 
manufacturer’s  name  for  the  instrument  mounted  there.  The  Skyvane  tower  is 
about  600  m  from  the  pier-end  tower.  The  other  land-based  platform  is  an 
antenna  tower  located  about  50  m  north  of  the  Skyvane  tower  and  equidistant 
from  the  beach.  Called  here  the  north  tower,  it  is  even  in  the  alongshore 
direction  with  Item  18  in  Figure  2  and  its  section  view  is  the  center  sketch 
in  Figure  3.  It  is  an  open-frame,  three-legged  telescoping  tower  with  legs 
varying  between  0.61  m  and  0.33  m  apart.  It  reaches  a  maximum  elevation  of 
about  22  m  above  msl. 

200.  Basic  wind  instrumentation  consisted  of  four  cup  anemometers,  four 
wind-vane  direction  sensors  and  one  impeller  vane,  and  a  combination  speed  and 
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horizontal  direction  sersor  (the  Skvvane).  Table  5  lists  some  of  the  attri¬ 
butes  of  these  devices.  The  cup-and-vane  sensor  sets  were  installed  at  two 
elevations  the  pier-end  tower  and  at  two  elevations  on  the  north  tower. 

The  Skvvane  was  at  an  elevation  of  18.7  m  above  msl,  so  cup  anemometers  were 
also  placed  at  this  elevation  on  the  pier-end  and  north  towers.  Wind  vanes 
were  installed  I  m  below  this  level,  at  17.7  m  above  msl.  Cup  anemometers 
were  also  placed  at  14.0  m  above  msl  on  the  pier-end  and  north  towers  to  in¬ 
vestigate  vertical  variations  in  mean  horizontal  wind  speed.  Wind  vanes  were 
also  mounted  1  m  below  these  at  13.0  tn  above  msl.  Lines  protruding  normal  to 
the  pier-end  and  north  towers  in  Figure  3  represent  the  mounting  booms  and 
show  these  sensor  locations.  Mounting  booms  were  at  least  three  characteris¬ 
tic  tower  diameters  long  to  minimize  flow  distortion  due  to  the  towers  proper. 
The  Skvvane  was  at  the  top  of  the  Skvvane  tower  shown  in  Figure  3. 

201.  For  gross  climatology,  the  air  temperature  over  land  was  measured 
at  site  18  in  figure  2,  which  represents  the  permanently  installed  FRF  weather 
station.  Sea  temperature  measurements  were  made  in  conjunction  with  deploy- 
ment  of  a  Sea  Data  PUV  wave  gage  at  site  11  in  Figure  2.  The  temperature 
sensor  of  this  gage  was  located  near  the  seabed  and  so  does  not  necessarily 
represent  sea-surface  temperature.  Air  temperature  over  water  was  not  mea¬ 
sured.  Temperature  differences  AT  (air  temperature  minus  water  temperature) 
give  a  crude  signature  of  possible  nearshore  thermal  forcing.  Some  of  the 
thermometer  characteristics  are  given  in  Table  5.  Characteristic  wave  heights 
were  obtained  through  standard  spectral  analysis  of  data  from  a  pressure  gage 
located  at  site  5  in  Figure  2  and  described  in  Table  5. 

Data  collection  and  editing 

202.  With  the  exception  of  the  Sea  Data  temperature  gage,  all  sensors 
were  sampled  at  2  Hz  for  40  min,  beginning  on  the  hour.  Data  collection  and 
processing  was  executed  using  a  Data  General  Nova  Model  1200  computer.  Mean 
values  and  standard  deviations  were  printed  soon  after  each  40-min  sampling 
interval.  Raw  time  series  were  stored  on  magnetic  tape.  For  this  study,  only 
mean  values  were  used  for  analysis.  Variances  were  used  for  error  checking. 
The  Sea  Data  gage  was  internally  recording.  Temperature  data  from  this  gage 
were  entered  in  a  master  data  file  after  the  experiment. 

203.  Prior  to  deployment,  all  wind-speed  sensors  were  calibrated  at  the 
National  Bureau  of  Standards  test  facility  in  Washington,  DC.  All  gages  were 
found  to  perform  t  •  within  manufacturer's  specifications.  Direction  sensors 
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Table  5 

Instrumental  Ion  llseil  in  the  Wind  Measurement  Study  During  the  DUCK81)  Experiment 
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Model  CP973C  for  40  min  teristic  wave  height 


were  assumed  accurate  to  manufacturer's  specifications.  They  were  aligned 
after  mounting  bv  locking  the  vanes  parallel  to  the  mounting  booms  and  aver¬ 
aging  the  measurements  for  40  min  in  this  configuration.  The  boom  azimuths 
were  found  through  triangulation  to  fixed  points  about  the  FRF  site,  the  loca¬ 
tions  of  which  were  determined  with  a  Zeiss  Elta-2  three-axis  theodolite.  In 
this  way,  absolute  directions  were  determined  to  within  il  deg.  which  is  better 
than  instrument  accuracy. 

204.  Data  were  collected  round  the  clock  with  the  exception  of  one  hour 
daily  for  computer  maintenance.  With  the  further  exception  of  infrequent 
problems  due  to  poor  connections  and  cross-talk,  data  return  was  excellent. 
Problems  were  detected  through  observed  radical  changes  in  mean  values  or 
variance.  To  ensure  that  all  data  were  above  instrument  threshold,  no  data 
were  retained  for  a  speed  and  direction  sensor  pair  where  the  speed  was  less 
than  1  m/sec. 

205.  For  the  general  climatology,  data  determined  to  be  satisfactory  in 
the  sense  of  instrument  response  were  used  in  toto.  However,  for  the  other 
objectives  of  the  experiment,  an  intercomparison  among  sensors  was  conducted. 
To  minimize  some  of  the  problems  of  platform  blockage  on  these  comparisons, 
data  were  edited  further  by  eliminating  from  analysis  all  speeds  and  direc¬ 
tions  for  which  the  sensors  were  directly  downwind  of  a  tower  or  large 
building . 

206.  Table  6  shows  the  ranges  of  azimuths  excluded  from  intercomparison 
analyses.  The  coordinate  system  used  in  this  report  is  a  compass  system 
(direction  increasing  clockwise)  with  the  effective  North  azimuth  directed 
northward  parallel  to  the  beach.  True  North  is  at  20  deg  in  this  coordinate 

Table  6 

Wind  Azimuths  for  Which  Data  Have  Been  Eliminated  from  Analysis 


Sensor  Location 

Range 

of  Azimuths,  deg 

Pier  end  (18.7  m) 

101  to  118 

Pier  end  (14.0  m) 

95  to  126 

North  tow’or  (18.7  m) 

114  to  178, 

0  to  8,  and  343  to 

360 

North  tower  (14.0  m) 

114  to  178, 

0  to  8,  and  343  to 

360 

Skyvane 

No  blockage 

assumed 
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s vs  tern. 


True  East  is  at  110  deg,  and  the  direction  seaward  along  the  pier  is 
90  deg.  Winds  are  considered  in  their  vector  sense  so  directions  given  are 
those  tow'ard  which  the  wind  blows.  Thus,  a  wind  blowing  southward  along  the 
beach  has  an  azimuth  of  180  deg.  A  wind  blowing  directly  onshore  (toward 
land)  has  an  azimuth  of  270  deg. 

Experiment  duration  and  climatology 

207.  Figure  4  gives  a  summary  of  some  of  the  general  climatology  of  the 
DUCK85  experiment.  Data  from  the  wind  and  other  sensors  have  been  tabulated 
by  Hubertz  et  al.  (1987).  Data  were  collected  more  or  less  continuouslv  from 
2100  (all  times  here  are  Eastern  Standard  Time)  on  5  September  1985  to  0700  on 
21  September  1985  when  the  experiment  was  scheduled  to  terminate.  However, 
during  the  cleanup  and  equipment  disassembly  period  following  the  experiment, 
it  became  obvious  that  the  path  of  Hurricane  Gloria  would  pass  near,  if  not 
over,  the  FRF.  To  sample  this  unique  event,  some  of  the  sensors  were  re¬ 
mounted  and  reconnected.  Additional  data  were  then  collected  from  1100  on 

26  September  1985  to  0700  on  28  September  1985. 

208.  Wind  data  are  shown  as  five  sets  of  stick  diagrams  in  Figure  4. 
From  the  beginning  of  the  experiment  until  11  September  the  winds  were  rela¬ 
tively  light,  generally  less  than  7  m/sec.  About  midday  on  11  September,  a 
front  passed  the  FRF  and  became  stationary  just  offshore.  Wind  speeds  were 
typically  in  the  range  of  10  to  15  m/sec  until  the  evening  of  14  '  ptember. 

The  influence  of  the  front  continued  until  18  September.  During  this  time  the 
winds  continued  to  blow  from  the  northeast  at  5  to  9  m/sec.  The  winds  became 
more  variable  from  18  to  21  September. 

209.  During  Hurricane  Gloria,  the  wind  followed  the  classic  pattern  for 
a  low-pressure  center  approaching  from  the  south.  The  winds  were  generally 
from  the  east  and  intensified  as  the  storm  center  neared  the  FRF.  Gloria 
passed  the  FRF  at  about  midnight  on  26  September.  The  winds  shifted  by  almost 
180  deg  in  5  hr  and  reached  a  peak  speed  of  about  23  m/sec.  It  is  noted  that 
the  stick  vectors  representing  the  Skyvane  appear  different  from  those  repre¬ 
senting  the  pier-end  tower  and  north  tower.  This  is  examined  in  detail 
shortly. 

210.  The  axis  in  Figure  4  marked  AT  shows  the  temperature  difference 
between  the  air  over  land  (near  the  FRF  laboratory  building)  and  the  water  at 
about  5-ra  depth  and  roughly  700  m  away.  The  record  is  incomplete  because  the 
Sea  Data  device  had  a  limited  duty  cycle.  The  measurements  are  somewhat  crude 


85 


86 


Figure  4.  Summary  of  datt  from  the  DUCK85  wind  measurement  experiment 


for  quantitative  interpretation  but,  nonetheless,  show  a  clear  diurnal  cycle. 
This  daily  variation  in  temperature  difference  can  be  visually  correlated  with 
daily  variations  in  the  wind  patterns  from  all  five  gages  and  under  all  wind 
conditions.  The  most  noticeable  effect  is  the  ordered  rotation  of  the  wind 
vectors  as  the  temperature  difference  intensifies  or  weakens.  This  suggests 
that  the  sea-breeze  forcing,  brought  about  by  land-sea  temperature  differ¬ 
ences,  acts  to  modulate  larger  scale  forcing  mechanisms.  The  effect  seems  to 
dominate  at  lower  wind  speeds.  During  the  northeast  storm  of  11-15  September, 
a  diurnal  variation  is  present  but  does  not  appear  to  be  a  major  modulating 
force . 

211.  Observed  wave  heights  are  shown  on  the  bottom  axis  of  Figure  4, 
They  are  characterized  as  four  times  the  standard  deviation  of  sea-surface 
elevation  where  the  latter  estimate  was  obtained  from  the  Sensa-Metric  pres- 
sur-  gage  described  above.  The  sea  was  relatively  calm  for  the  first  few  days 
of  the  experiment  but  wave  heights  grew  rapidly  at  the  onset  of  higher  winds 
from  the  northeast  storm.  Characteristic  wave  heights  reached  about  1.5  m  (at 
the  location  of  this  sensor)  at  the  height  of  the  storm  as  compared  with  0.4  m 
prior  to  the  storm.  Wave  heights  stayed  high  after  the  main  part  of  the  storm 
passed  on  15  September,  due  partly  to  forcing  by  6-  to  8-m/sec  wind?  and  prob¬ 
ably  due  to  added  storm  waves  propagating  landward  from  the  somewhat  distant 
front.  Wave  heights  increased  again  during  the  passage  of  Hurricane  Gloria, 
approaching  saturation  for  the  water  depth  in  which  the  gage  was  located. 
However,  Gloria  passed  the  FRF  site  relatively  quickly  and  the  wind  blew  sea¬ 
ward  (opposite  the  propagation  direction  of  landward-traveling  storm  waves) 
with  reasonably  high  intensity  after  the  storm  passage  so  wave  heights  dropped 
off  rapidly. 

212.  These  observations  indicate  two  important  features  of  nearshore 
wind  coupling  to  the  sea.  First,  the  wave  field  tends  to  grow  rapidly  at  the 
onset  of  high  winds.  This  suggests  that  instantaneous  sea  state  conditions 
are  in  equilibrium  with  the  wind  when  seas  are  building.  A  wave-growth  model 
could  thus  be  designed  using  this  concept.  Second,  in  decaying  winds  the  seas 
can  remain  high  due  to  the  relatively  long-lived  nature  of  waves  and  their 
ability  to  radiate  from  the  region  of  generation.  In  this  case,  winds  and 
waves  are  not  in  equilibrium  so  that  coupling  (drag)  coefficients  can  be  ex¬ 
pected  to  vary  greatly,  i.e.,  not  be  coupled  uniquely  to  instantaneous  sea 
state.  Accurate  modeling  of  a  drag  coefficient  must  rely  on  observations 
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under  a  wide  variety  of  conditions  like  these  because  the  physics  of  these 
interactions  is  quite  complex  and  is  not  understood  fully. 

213.  The  sequence  of  wind  vectors  denoted  by  the  stick  diagrams  in 
Figure  4  indicates  that,  for  this  set  of  observations,  the  wind  field  gener¬ 
ally  satisfies  the  conditions  of  quasi-steadiness  given  by  scaling  arguments 
in  the  last  section.  Two  minor  exceptions  to  this  are  noted.  At  low  wind 
speeds,  less  than  5  m/sec  and  for  which  little  wave  generation  occurs,  the 
speeds  often  change  by  more  than  10  percent  per  hour  and  directions  change  bv 
more  than  20  deg  per  hour.  For  low-energy  conditions  this  is  not  critical. 
Also,  at  the  onset  of  the  northeast  storm  on  11  September  the  wind  increased 
and  changed  directions  somewhat  rapidly.  However,  since  it  was  a  transition 
from  low-  to  high-wave  generation  conditions,  it  is  only  the  onset  of  high 
winds  that  is  important.  The  nature  of  the  transient,  low-wind  conditions 
preceding  this  is  not  too  critical. 

214.  On  the  other  hand,  a  major  exception  to  the  quasi-steady  scaling 
occurs  during  the  passage  of  Hurricane  Gloria  on  26  to  28  September.  Durir.s 
the  storm  passage,  the  winds  shifted  at  rates  approaching  90  deg  per  hour.  In 
3  hr  the  winds  shifted  by  almost  180  deg.  Thus,  while  wind  speeds  were  high 
enough  ( > 20  m/sec)  for  efficient  wave  generation,  the  rapid  change  in  wind 
direction  resulted  in  either  the  wind  blowing  nearly  perpendicular  to  the 
direction  of  propagation  of  waves  generated  the  previous  hour  or  the  wind 
blowing  opposite  the  direction  of  propagation  of  waves  generated  several  hours 
previously.  In  the  latter  case,  momentum  transfer  is  likely  to  be  from  waves 
to  wind  (or  negative  momentum  transfer  from  wind  to  waves)  with  a  consequent 
loss  of  wave  energy.  A  quantification  of  the  rate  of  momentum  transfer  (i.e., 
a  drag  coefficient)  is  likely  to  be  quite  different  from  that  for  growing 
waves,  but  its  magnitude  is  not  known.  Clearly,  this  is  an  important  topic 
for  further  research. 

Correlations  between  sensors 

215.  The  primary  analysis  of  this  data  set  is  done  through  correlation 
plots  comparing  various  pairs  of  sensors.  Relative  sensor  behavior  and  wind- 
field  differences  are  seen  readily  in  this  way.  Of  major  importance  are  com¬ 
parisons  of  wind  speeds.  Since  wind-speed  differences  are  likely  to  be  strong 
functions  of  wind  direction,  it  is  useful  to  classify  the  data  by  wind 
direction. 

216.  This  was  done  by  dividing  the  azimuth  circle  into  12  equal  parts 
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as  shown  in  Figure  5.  Each  sector  is  30  deg  wide.  Four  of  the  sectors  are 
centered  on  cardinal  points  in  the  pier-aligned  coordinate  system.  Symbols 
are  shown  in  Figure  5  to  define  uniquely  each  sector.  Triangles  designate 
cardinal  points  (in  the  pier  system)  and  can  be  thought  of  as  arrowheads 
pointing  in  the  direction  toward  which  the  wind  blows.  Intervening  symbols 
also  represent  directions  toward  which  the  wind  blows,  e.g.,  the  plus  symbol 
represents  winds  blowing  toward  15  to  45  deg  clockwise  from  the  y-axis.  The 
y-axis  is  shore-parallel  and  positive  northward.  The  x-axis  is  pier-parallel 
and  positive  seaward. 


LEGEND 

Figure  5.  Symbol  legend  showing  wind-direction 
bands  and  representative  symbols 

217.  Correlations  between  the  three  wind-speed  and  wind-direction 
sensor  sets  at  the  18.7-m  level  are  shown  in  Figures  6-8.  Symbols  denote  wind 
direction,  and  axes  denote  wind  speeds.  Figure  6  shows  pier-end  winds  com¬ 
pared  with  north  tower  winds.  There  is  some  scatter  in  the  data,  and  some 
trends  can  be  seen  within  symbol  groupings.  This  will  be  discussed  in  detail 
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Skyvane  Wind  Speed  (m/sec) 


Figure  6.  Correlation  plot  of  north  tower  18.7-m 
wind  speed  compared  with  pier-end  18.7-m  wind  speed 
using  data  from  all  direction  bands 


Pier-End  18.7-m  Wind  Speed  (m/sec) 


Figure  7.  Correlation  plot  of  Skyvane  wind 
speed  compared  with  pier-end  18.7-m  wind 
speed  using  data  from  all  direction  bands 


Figure  8.  Correlation  plot  of  Skvvane  wind 
speed  compared  with  north  tower  18.7-m  wind 
speed  using  data  from  all  direction  bands 


below.  The  general  observation  in  Figure  6  is  that  the  bulk  of  the  data  agree 
to  within  about  10  percent.  If  these  sensors  were  over  the  open  ocean,  one 
would  expect  a  perfect  correlation.  The  scatter  indicates  that  the  two 
sensors  are  not  detecting  the  same  wind  pattern.  The  difference,  of  course, 
is  that  one  sensor  is  over  land  and  the  other  is  over  water;  the  boundary 
layer  structure  in  and  between  the  two  regimes  is  different. 

218.  The  general  difference  shown  in  Figure  6  is  less  than  about 

10  percent  over  a  horizontal  separation  of  about  600  m.  Based  on  the  scaling 
arguments  given  in  the  section  on  horizontal  nonuniformity,  this  suggests  that 
the  flow  is  uniform  to  a  roughly  acceptable  degree  between  the  two  stations 
for  most  of  the  data  at  the  18.7-m  level.  That  is,  wind  measurements  from  the 
north  tower  are  roughly  representative  of  winds  over  water  600  m  away.  The 
more  detailed  analysis  below  indicates  that  improved  representations  are 
possible . 

219.  Figure  7  shows  a  correlation  of  pier-end  winds  with  those  detected 
by  the  Skyvane.  At  wind  speeds  less  than  about  10  m/sec  the  general  data 
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pattern  in  Fieure  '  is  about  the  same  as  in  Figure  6.  However,  for  speeds  be¬ 
tween  '0  an;!  !5  m/sec  the  Skvvane  tends  to  register  slightly  higher  wind 
speeds  than  sensors  at  the  pier  end.  This  should  be  compared  with  the  north 
tower  which  tends  to  register  sligitly  low  relative  to  pier-end  sensors. 
Alarmingly,  the  Skvvane  does  not  register  speeds  i”  excess  cf  13  m/sec  even 
though  the  pier-end  anemometer  senses  speeds  exceeding  22  m./sec.  In  fact,  for 
the  higher  pier-en.d  speeds  the  Skvvane  does  not  read  above  11  m/sec;  this 
behavior  is  not  seen  in  data  from  the  north  tower. 

220.  Figure  8  shows  the  correlation  between  the  north  tower  18.7-m  wind 
speed  and  the  Skvvane.  Since  both  towers  are  equidistant  from  the  beach  their 
readings  should  be  identical.  The  correlation  is  quite  good  for  speeds  up  to 
about  9  m/s ec.  From  this  to  speeds  of  about  15  m/sec  the  Skvvane  appears  to 
overspeed  slightly  and,  as  with  the  pier-end  comparison,  drop  considerably  for 
north  tower  winds  in  excess  of  15  m/sec. 

221.  There  are  two  possible  explanations  for  this  behavior.  One  is 
that  the  sensor  malfunctioned.  Though  it  calibrated  to  within  manu'ac turer ’ s 
specifications  prior  to  the  experiment  and  behaved  properly  both  before  and 
after  Hurricane  Gloria  (compare  the  stick  vectors  in  Figure  4  for  27  and 

28  September) ,  it  is  possible  that  the  sensor  signal  was  degraded  during  the 
high  winds  of  the  hurricane.  A  more  likely  explanation  is  that  the  wind  was 
strongly  affected  by  the  bluff  upwind  faces  of  the  FRF  building  such  that  the 
Skvvane  sensor  was  in  a  low-velocity  wake  region  of  flow  at  high  wind  speeds. 
The  svmbols  for  the  highest  wind  speeds  indicate  air  flow  directly  onshore. 

The  Skvvane  platform  section  shown  in  Figure  3  indicates  that,  for  onshore 
winds,  the  air  passed  up  the  seaward  dune  face  and  immediately  had  to  pass 
over  or  around  the  FRF  building. 

222.  There  is  no  simple  theory  for  flow  separation  around  angular  ob¬ 
structions  (like  the  FRF  building).  However,  Batchelor  (1970)  notes  that 
fluids  tend  to  have  a  stagnation  zone  downstream  of  sharp  edges,  and  that  this 
zone  increases  with  increasing  Reynolds  number  (governed  by  wind  speed  in  this 
case).  He  also  notes  that  fluid  flow  tends  to  follow  a  tangent  to  the  upwind 
surface  of  an  obstruction.  In  the  case  of  the  FRF  building,  the  upwind  sur¬ 
face  is  vertical  so  the  effect  would  be  expected  to  be  dramatic.  The  apparent 
slight  overspeeding  of  the  Skvvane  at  speeds  of  9  to  15  m/sec  is  also  indica¬ 
tive  of  flow-acceleration  effects  around  an  obstacle. 

223.  While  more  detailed  mapping  of  the  wind  field  over  and  around  the 
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FRF  building  at  high  wind  speeds  is  necessary  to  resolve  this  problem,  it  is 
evident  thac  the  Skvvane  sensor  did  not  detect  the  same  wind  field  as  the 
north  tower  sensor.  Since  instrument  malfunction  or  platform  distortion  are 
suspected  to  be  the  cause  of  this,  the  Skvvane  data  are  not  considered  further 
in  this  report.  It  is  recommended  that  the  Skvvane  platform  and  platforms 
like  it  be  studied  in  detail  to  ensure  that  they  do  not  induce  grievous  errors 
in  field  observations. 

224.  Figure  9  is  a  correlation  plot  of  wind  direction  at  the  pier-end 


Figure  9.  Correlation  plot  of  north  tower 
17.7-m  wind  direction  compared  with  pier- 
end  17.7-m  wind  direction 


tower  with  that  at  the  north  tower.  Most  of  the  data  is  within  one  symbol 
width  (about  6  deg)  of  the  perfect  correlation  line.  The  only  clear  deviate 
trend  occurs  for  winds  with  an  azimuth  between  270  and  360  deg,  i.e.,  winds 
heading  northwest  or  coming  from  the  southeast.  Winds  in  this  direction  pass 
over  the  pier  structure  itself  (see  Figure  2)  prior  to  reaching  the  north 
tower.  It  is  possible  that  the  flow  in  the  wake  of  the  pier  may  have  had  some 
effect  on  wind  direction  at  the  north  tower.  Another  possibility  is  that  the 
vane  on  the  north  tower  was  deflected  by  flow  around  the  tower.  The  vane 
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would  be  nearest  the  tower  at  the  instrument  boom  azimuth  of  336  deg.  Flow 
deflection  around  the  tower  would  tend  to  result  in  a  reduced  reading  for  the 
north  tower.  This  is  consistent  with  deviations  seen  in  Figure  9.  Aside  from 
this  possible  platform  effect,  no  major  direction  deviation  trends  are  evident 
in  this  data  set. 

225.  The  scatter  of  points  shown  in  Figure  6  is  due,  at  least  in  part, 
to  the  combining  of  wind  data  from  all  directions  on  a  single  graph.  The  wind 
system  at  the  FRF  site  is  more  orderly  than  is  suggested  by  Figure  6.  This  is 
shown  in  Figures  10  to  17  wh  re  pier-end  and  north  tower  18.7-m  wind  speeds 
are  correlated  for  each  of  eight  azimuthal  bands.  The  four  azimuthal  bands 
not  shown  had  too  little  data  to  be  instructive.  Figures  10  to  17  are  ordered 
by  increasing  wind  azimuth  such  that  Figures  10  to  12  contain  winds  with  an 
offshore-directed  component  and  Figures  13  to  17  contain  winds  with  an 
onshore-directed  component  (see  Figure  5  for  symbol  legend). 

226.  The  greatest  fractional  wind  deviation  occurs  in  Figure  10.  Here 
the  wind  is  headed  northeasterly.  Examination  of  Figure  1  shows  that  these 
winds  have  traveled  the  greatest  distance  over  land  before  reaching  either 
tower  and,  further,  wind  at  the  pier-end  tower  has  passed  the  greatest  dis¬ 
tance  over  water  (for  winds  with  an  offshore  component)  prior  to  reaching  the 
pier-end  tower.  Thus,  one  would  expect  the  highest  contrast  between  the 
sensors  in  this  direction  (and  its  symmetric  counterpart  about  the  x-axis, 
i.e.,  the  circle-x  symbols,  not  shown  because  sensors  were  in  the  wakes  of 
towers) .  The  over-land  boundary  layer  would  be  most  fully  developed  and 
register  maximum  land-roughness  effects  in  the  retardation  of  wind  flow.  Wind 
at  the  pier  end  would  have  passed  generally  over  the  most  quiescent,  and 
therefore  smoothest,  water  surface  since  fetch  limitations  would  inhibit  any 
significant  local  wave  growth.  This  sensor  would  then  show  the  least  retarda¬ 
tion.  The  contrast,  about  25  percent  for  an  8-m/sec  pier-end  wind  speed,  is 
significantly  large.  A  wind  stress  based  on  the  land-derived  (north  tower) 
wind  would  be  about  50  percent  low  using  a  quadratic  stress  law. 

227.  In  the  next  azimuthal  band  (x  symbols,  Figure  11)  the  north  tower 
wind  also  tends  to  read  systematically  low  relative  to  the  pier  end.  Here  the 
maximum  deviation  is  only  about  15  percent.  This  is  because  the  north  tower 
air  has  passed  less  distance  over  land  than  for  the  azimuthal  band  discussed 
in  the  last  paragraph.  Also,  the  pier-end  wind  has  passed  less  distance  over 


94 


Pier-End  18  7-m  //ind  Speed  im/sec.) 

Figure  10.  Correlation  plot  of  north  tower  18.7-m  wind 
speed  compared  with  pier-end  18.7-m  wind  speed  for  winds 
blowing  toward  15  to  45  deg  in  the  pier  coordinate  system 


Pier-End  18.7-m  Wind  Speed  (m/sec) 

Figure  11.  Correlation  plot  of  north  tower  18.7-m  wind 
-speed  compared  with  pier-end  18.7-m  wind  speed  for  winds 
blowing  toward  45  to  75  deg  i..  the  pier  coordinate  system 
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Figure  12.  Correlation  plot  of  north  tower  18.7-m  wind 
speed  compared  with  pier-end  18.7-m  wind  speed  for  winds 
blowing  toward  75  to  105  deg  (generally  offshore)  in  the 
pier  coordinate  system 


0.0  5.0  10.0  15.0  20.0 


Pier-End  18.7-m  Wind  Speed  (m/sec) 

Figure  13.  Correlation  plot  of  north  tower  18.7-m  wind 
speed  compared  with  pier-end  18.7-m  wind  speed  for  winds 
blowing  toward  195  to  225  deg  in  the  pier  coordinate  system 


0.0  5.0  10.0  15  0  20.0 

Pier-End  18.7-m  Wind  Speed  (m/sec) 


Figure  .14.  Correlation  plot  of  north  tower  18.7-m  wind 
speed  compared  with  pier-end  18.7-m  wind  speed  for  winds 
blowing  toward  225  to  255  deg  in  the  pier  coordinate  system 


figure  15,  Correlation  plot  of  north  tower  18.7-m  wind 
speed  compared  with  pier-end  18.7-m  wind  speed  for  winds 
blowing  toward  725  to  285  deg  (.generally  onshore)  in  Che 
pier  coordinate  system 
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Figure  16.  Correlation  plot  of  north  tower  18.7-m  wind 
speed  compared  with  pier-end  18.7-m  wind  speed  for  winds 
blowing  toward  285  to  315  deg  in  the  pier  coordinate  system 


Pier-End  18.7-m  Wind  Speed  (m/sec) 

Figure  17.  Correlation  plot  of  north  tower  18.7-m  wind 
speed  compared  with  pier-end  18.7-m  wind  speed  for  winds 
blowing  toward  315  to  345  deg  in  the  pier  coordinate  system 
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water  and  so  still  retains  some  land-derived  turbulence  (see  the  section  on 
horizontal  nonuniformity  in  Part  III). 

228.  Figure  12  shows  the  correlations  for  the  azimuthal  band  centered 
on  the  offshore-directed  axis.  There  is  some  scatter,  but  the  systematic 
pattern  is  closer  to  the  perfect  correlation  line  than  in  either  of  the  two 
previous  cases.  Here  the  wind  has  passed  a  minimum  distance  over  land  at  the 
north  tower  and  a  minimum  distance  over  water  at  the  pier  end.  Evidently  the 
land-influenced  flow  has  not  adjusted  to  over-water  conditions  in  this  short 
(600  m)  distance  so  that  for  speeds  in  excess  of  about  4  m/sec  the  two  wind 
measurements  are  about  the  same.  The  two  deviate  readings  at  speeds  less  than 
2  m/sec  were  obtained  at  times  of  light  and  variable  winds.  Threshold  prob¬ 
lems  and  high  variances  in  instantaneous  speeds  and  directions  used  to  con¬ 
stitute  mean  values  may  have  caused  the  deviations  seen. 

229.  Figures  13  to  17  show  correlations  for  winds  with  an  offshore 
component.  Symmetry  about  the  x-axis  in  Figure  5  suggests  that  Figures  13 
and  17  should  be  identical  and  that  Figures  14  and  16  should  also  be  identi¬ 
cal.  To  within  the  data  scatter  at  wind  speeds  less  than  about  6  m/sec  (the 
maximum  speeds  attained  in  Figures  16  and  17),  this  is  seen  to  be  the  case. 

For  this  reason,  references  to  Figures  13  and  14  will  imply  references  to 
Figures  16  and  17,  respectively,  in  the  following  discussion. 

230.  Figure  13  shows  that  for  winds  with  a  slight  onshore  component  the 
correlation  is  very  good  to  speeds  of  about  10  m/sec.  For  higher  speeds,  the 
north  tower  winds  are  retarded  somewhat,  in  general.  The  point  of  maximum 
deviation  (about  16  m/sec  at  the  pier  end)  occurred  during  Hurricane  Gloria  at 
a  time  of  rapid  temporal  change  and  so  could  be  due,  in  an  unknown  way,  to 
flow  acceleration  effects.  The  lesser  but  systematic  deviations  for  the 
higher  wind  speeds  are  likely  due  to  the  influence  of  air  flow  over  the  dune 
shown  in  Figure  3.  A  line  tangent  to  the  maximum  slope  of  the  seaward  dune 
face  passes  above  the  1  .7-m  level  of  the  north  tower.  The  argument  given  for 
the  Skyvane  platform  also  applies  here.  At  some  wind  speed,  it  is  expected 
that  air  flow  due  to  the  wake  of  the  dune  will  reach  the  18. 7-m  elevation.  In 
Figure  3  this  appears  to  occur  for  wind  speeds  between  10  and  13  m/sec.  The 
increased  scatter  at  a  speed  of  about  13  m/sec  may  also  be  due  to  this. 

231.  The  dune  effect  is  expected  to  become  evident  at  lower  wind  speeds 
as  winds  become  more  shore-normal  since  such  winds  have  a  stronger  onshore 
component.  This  is  evident  in  Figure  14  where  a  slight  overspeeding  of  the 
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north  tower  signal  at  speeds  around  7  m/sec  suggests  flow  acceleration  over 
the  dune  crest.  The  increase  in  scatter  occurs  at  speeds  around  10  m/sec  com¬ 
pared  with  !3  m/sec  in  the  last  case  discussed.  In  Figure  15  the  overspeeding 
of  the  north  tower  winds  occurs  at  speeds  around  5  m/sec  and  the  increase  in 
scatter  at  about  7  m/sec.  Winds  in  Figure  15  are  blowing  within  ±15  deg  of 
directly  onshore,  and  so  would  show  dune  wake  effects  at  the  lowest  wind 
speeds.  Note  in  Figure  15  that  the  higher  wind  speeds  encountered  during 
Hurricane  Gloria  resulted  in  about  10-percent  retardation  of  the  north  tower 
winds  relative  to  those  at  the  pier  end. 

232.  If  the  conjecture  is  true  in  the  foregoing  discussion  about  varia¬ 
tions  in  the  wind  pattern  at  the  18.7-m  level,  then  deviations  would  be  even 
more  pronounced  at  lower  levels.  Figure  18  shows  the  general  comparison  of 


Figure  18.  Correlation  plot  of  north  tower 
14.0-m  wind  speed  compared  with  pier-end 
14.0-m  wind  speed 

pier-end  and  north  tower  wind  speeds  at  the  14.0-m  level,  analogous  to  Fig¬ 
ure  6.  Comparison  of  Figures  6  and  18  shows  that  the  deviation  is,  in  fact, 
stronger  at  the  14.0-m  level.  Though  Figure  18  is  not  disassembled  by  wind 
azimuth  classes  as  was  done  for  Figure  6,  the  trends  of  some  of  the  classes 
are  clear.  For  instance,  wind  blowing  toward  15  to  45  deg  (plus  symbol) 
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deviates  at  the  north  tower  by  almost  40  percent  from  pier-end  winds  of  about 
8  m/sec.  This  is  somewhat  larger  than  the  25-percent  difference  found  at  the 
18.7-m  level.  Also,  winds  blowing  toward  195  to  225  deg  (boxed  plus  symbol) 
show  about  10-percent  attenuation  at  the  north  tower  in  Figure  18  as  compared 
with  roughly  5  percent  in  Figure  6. 

233.  One  result  of  this  analysis  is  that,  though  there  are  differences 
between  sensors  separated  by  600  m  across  the  nearshore  zone,  the  differences 
are  highly  systematic  in  wind  speed  and  direction.  If  other  potentionallv  im¬ 
portant  parameters,  such  as  air-sea  and  sea-land  temperature  differences  or 
sea-state  variables,  were  included,  wind  differences  may  be  even  more  highly 
systematized.  Such  additional  analysis  was  not  done  here  because  of  the 
paucity  of  data.  Nonetheless,  this  means  that  very  good  corrections  to  land- 
based  observations  are  possible  (for  this  site  and  between  the  two  observation 
platforms  considered)  by  way  of  direct  correlations.  Moreover,  wind  differ¬ 
ences  are  explainable  qualitatively  in  terms  of  recognized  boundary  layer 
behavior.  This  suggests  that  further  research  can  enable  good  numerical 
modeling  of  nearshore  winds.  Such  modeling  would  have  broader  application 
than  at  this  single  experiment  site  and  therefore  be  more  generally  useful. 

234.  The  last  objective  in  this  experiment  is  to  evaluate  the  pier-end 
tower  for  adequacy  in  making  wind-profile  measurements.  The  test  performed 
was  a  comparison  of  wind  speeds  at  two  elevations  (18.7  m  and  14.0  m)  on  the 
pier-end  tower  to  examine  variations  in  wind  shear.  Figure  19  is  a  correla¬ 
tion  plot  of  the  speeds  recorded  at  the  two  elevations.  This  result  is  very 
remarkable.  The  plot  shows  chat  the  14.0-m  wind  speeds  are  less  than  the 
18.7-m  wind  speeds  by  about  0.2  m/sec  for  all  wind  speeds  and  directions;  this 
is  not  an  expected  behavior.  If  the  flow  was  neutrally  stratified,  such  that 
the  logarithmic  profile  (Equation  5)  is  valid,  then  the  velocity  difference 

L’2  -  between  two  elevations  z^  and  z^  ,  would  be 

=  (u./<)ln(z  /z  )  .  Since  k  ,  z  ,  and  z,  are  fixed  in  this  case  at 
*  2  l  2  1 

0.41  m,  18.7  m,  and  14.0  m,  respectively,  then  the  result  *  0.72  u* 

should  hold.  If  -  Uj  is  constant,  then  u^  must  be  constant,  by  this 

derivation,  and  tne  wind  stress  must  also  be  constant. 

235.  For  Uj  -  Uj  *  0.2  m/s ec,  the  above  expression  gives  u* 

~0.3  m/sec.  This  is  a  reasonable  value  for  a  general  estimate  of  u^  . 
However,  according  to  Garratt  (1977),  uA  should  be  0.16  and  0.64  m/sec  for 
wind  speeds  of  5  and  15  m/sec,  respectively.  By  the  expression  given  above, 
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Figure  19.  Correlation  plot  of  pier-end  wind 
speeds  at  14.0-m  compared  with  those  at  18.7-m 


l\  -  U  should  then  be  0.12  and  0.46  m/sec.  respectively.  Figure  19  shows 
that  this  is  not  the  case,  so  something  is  amiss.  It  is  possible  that  the 
instruments  malfunctioned  or  were  recorded  incorrectly.  It  i6  also  possible 
that  the  sensor  at  the  14.0-m  elevation  was  detecting  a  wind  field  distorted 
by  structures  and  appendages  at  the  pier  deck  level.  Figure  3  shows  a  small 
building  of  about  2  m  in  height  near  the  base  of  the  pier-end  tower.  The  rule 
of  thumb  for  flow  around  isolated  obstacles  is  that  the  flow  is  distorted 
within  three  characteristic  length  scales  of  the  center  of  obstacle.  For  the 
size  and  location  of  the  small  building  at  the  end  of  the  pier,  this  indicates 
that  measurements  up  to  the  15-m  elevation  on  the  pier-end  tower  would  be 
affected.  Note  that,  except  in  the  wake  of  the  small  building,  the  flow  would 
tend  to  be  accelerated  in  the  vicinitv  of  the  obstruction.  This  would  reduce 
the  velocity  difference  between  the  two  measuring  heights  on  the  pier-end 
tower.  Measurements  in  the  wake  of  the  tower  were  edited  out  before  analysis. 
This  excluded  measurements  in  the  wake  of  the  small  building,  so  only  accel¬ 
erated  measurements  are  seen  in  Figure  19. 

236.  This  explanation  of  observations  leads  to  the  conclusion  that, 
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unless  some  par t icularlv  imaginative  mechanical  engineering  is  done  to  mount 
instruments  veil  away  from  pier  blockage  effects  at  lower  elevations,  it  would 
be  a  grave  mistake  to  attempt  profile  measurements  from  the  p^er-end 
structure . 

Summary  of  results 

237.  Results  of  the  preliminary  wind  experiment  in  the  fall  of  1985 
pr'vided  valuable  insight  into  the  mechanics  of  the  nearshore  wind  field  at 
the  experiment  site  and  guidance  for  the  design  of  the  main  experiment,  year- 
shore  winds  appear  to  vary  in  the  horizontal  in  a  highly  systematic  way,  de¬ 
pending  strongly  on  wind  speed  and  direction.  This  means  that  long-term 
correlations  between  wind  speeds  and  directions  at  two  points  at  a  given  site 
can  be  used  to  deduce  an  empirical  transfer  function  for  land-based  winds  to 
over-water  winds.  It  also  suggests  that,  with  further  study,  a  reasonably 
simple  model  could  be  constru  '.ted  to  predict  the  transformation  for  a  more 
arbitrary  site.  Caution  is  advised  in  arbitrary  placement  of  a  land-based 
anemometer.  Results  from  this  experiment  indicate  that  the  rule  of  thumb  of 
three  characteristic  length  scales  of  upwind  and  lateral  separation  of  instru¬ 
ments  from  possible  flow-distorting  obstructions  is  a  strict  rule.  The  data 
analysis  suggests  that  measurements  from  the  building  top  anemometer  and  the 
lower  pier-end  anemometer  were  both  affected  by  nearby  structures.  The  analy¬ 
sis  also  suggests  that  the  pier-end  tower  should  not  be  used  for  wind-stress 
estimation  by  the  profile  technique. 


SUPERDUCK 

238.  In  view  of  the  results  of  the  DUCK85  experiment,  it  was  decided  to 
use  the  direct  method  for  measuring  wind  stress  over  water  during  the  fall  of 
1986  experiment  known  as  SUPERDUCK.  This  endeavor  was  a  joint  venture  with 
CERC  and  ("with  CERC  funding)  the  Department  of  Marine,  Earth  and  Planetary 
Sciences  of  North  Carolina  State  University  (NCSU)  in  Raleigh,  North  Carolina. 
Instrumentation  type,  method  of  deployment,  data  collection,  and  data  editing 
were  mostly  selected,  designed,  or  performed  by  NCSU  personnel.  Instruments 
were  provided  by  NCSU.  Data  analysis  is  currently  in  progress  and  will  be 
described  in  a  future  report.  The  intent  of  this  section  of  this  report  is  to 
describe  the  experiment  design,  instrumentation,  and  data  collection  logistics 
associated  with  the  wind  measurements. 
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Measurement  requirements 


23Q.  In  direr1-  flux  measurements  where  stratification  is  important,  it 
is  necessarv  to  measure  fluctuating  temperature,  humidity,  and  three  compo¬ 
nents  cf  velocity.  Cross  correlations  of  these  quantities  provide  a  measure 
of  the  Revnolds  fluxes  of  heat,  humidity,  and  momentum.  Relating  these  ob¬ 
servations  to  the  bulk  method  requires  that  mean  horizontal  wind  speed,  wind 
direction,  air-sea  temperature  difference,  and  air-sea  humidity  difference  be 
measured  (or  estimated)  as  well.  These  fluxes  and  mean  parameters  are  the 
primary  variables  sought  in  this  experiment. 

240.  In  analysis  of  data  from  this  experiment  it  is  necessary  to  have 
measurements  of  air  temperature  over  land,  sea  state  variables  (at  least  char¬ 
acteristic  wave  height,  wave  period,  and  dominant  travel  direction)  and  sea- 
surface  current.  These  variables  are  necessary  to  test,  calibrate,  or  modify 
the  hypotheses  and  models  given  in  Part  II  of  this  report.  Measurements  of 

all  these  variables  were  made  during  Sl'PERDUCK  by  other  groups  of  invest iga-  < 

tors.  They  have  agreed  to  share  their  data  for  use  in  wind-stress  analysis. 

Platform 

241.  The  primary  measurements  in  the  wind  measurements  part  of 
Sl'PERDUCK  were  all  made  at  or  very  near  the  seaward  end  of  the  FRF  pier 
(station  12  in  Figure  2).  Three  instrument  mounting  booms  were  attached  to 
the  pier-end  tower.  A  2.4-m  retractable  main  boom  was  located  at  18.7  m  above 
msl,  corresponding  in  location  to  the  upper  boom  shown  for  the  pier-end  tower 
in  Figure  3.  A  2.0-m  fixed  boom,  called  the  top  boom,  was  located  at  22  m 
above  msl,  corresponding  to  the  top  of  the  pier  tower  in  Figure  3.  A  second 
2.0-m  fixed  boom,  called  the  lower  boom,  was  located  at  13.7  m  above  msl, 
corresponding  (about)  to  the  lower  of  the  lower  pair  of  booms  on  the  pier-end 
tower  in  Figure  3. 

242.  A  sea-temperature  gage  was  mounted  on  the  second  pier  piling  from 
the  pier  end  at  depths  varying  between  0.6  and  2.0  m  below  msl.  Depths  varied 
because  several  mounting  schemes  were  tried  before  a  satisfactory  arrangement 
was  found.  Measurements  from  this  gage  were  recorded  automatically  and  were 
checked  three  times  a  day  with  a  bucket  thermometer  lowered  from  the  pier  deck 
at  a  location  generally  midway  between  the  last  two  pier  pilings. 

243.  Hand-held  air  temperature  and  relative  humidity  instruments  were 

also  operated  three  times  daily  at  about  9  m  above  msl  (roughly  waist  high  l 

above  the  pier  deck  in  the  lee  of  the  instrument  trailer).  This  trailer  is 
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schematised  in  Figure  3  as  a  rectangle  centered  over  '  h 
the  pier  end. 


second  piling  from 


Instrumentat  lor. 

2-iA.  Thirteen  sensors  were  deployed  to  provide  primary  and  redundant 
measurements  of  mean  and  fluctuating  variables  of  the  wind  field.  Sensors 
were  divided  into  groups  depending  on  the  frequencies  at  which  their  outputs 
were  sampled.  Data  from  each  group  of  sensors  were  recorded  on  separate  data 
logging  devices.  The  two  sampling  frequencies  were  5  and  10  Hz. 

345.  The  low-frequency  (5  Hz)  sensors  were  intended  primarily  to  mea¬ 
sure  mean  quantities  for  use  both  in  bulk  parameterization  and  in  checking  for 
mean,  drift  in  the  high-frequency  sensors.  Table  7  lists  the  variables 
measured,  the  devices  used,  the  number  of  data  channels,  and  the  locations  of 
the  low-f requencv  sensors.  These  measurements  include  three  orthogonal 


Table  7 

Devices  Sampled  at  Low  Frequency  (5  Hz)  During  Wind 
Measurements  Component  of  SUPERDUCK 


Var iab le 

Number  of 
Channels 

Manufacturer/ 
Model/Sensor  Type 

Location 

Elevation  above  msl 

> 

Horizontal  wind 
speed  and 

>  direction 

2 

R.  M.  Young 

Model  05103 

Impellor  vane 
anemometer 

22.0 

m  (top  boom) 

Vertical 
wind  speed 

1 

R.  M.  Young 

Model  27106 

Gill  impellor 
anemometer 

22.0 

m  (top  boom) 

Air  temperature 
over  water 

1 

Yellow  Springs 

Model  205 

Thermistor 

18.7 

m  (main  boom) 

Humid  it v 

over  water 

1 

Campbell  Scientific 

Model  207 

Phys-Chemical  Research 
Relative  humidity  sensor 

18.7 

m  (main  boom) 

Horizontal 
wind  speed 

1 

Met-One 

Model  010B 

Cup  anemometer 

13.7 

m  (lower  boom) 

'  Sea  temperature 

1 

Yellow  Springs 

Model  205 

Thermistor 

-0.6 

ro  to  -2.0  m 
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components  of  wind  speed  at  the  top  (22-m  elevation)  of  the  pier-end  tower, 
mean  air  temperature  and  humidity  at  the  level  of  the  main  boom,  horizontal 
wind  speed  a:  the  lower  boom  (13.7-m  elevation),  and  sea  temperature. 

216.  The  high-frequency  (10  Hz)  devices  were  intended  to  resolve  the 
fluctiating  (turbulent)  components  of  temperature,  humidity,  and  all  three 
speed  components  of  the  wind.  These  sensors  were  all  mounted  on  the  main  boom 
within  a  volume  of  characteristic  dimension  1  m.  Cross  correlations  between 
sensors  could  then  be  used  to  estimate  Reynolds  fluxes  of  mass  and  momentum. 
Variables  measured,  devices  used,  and  numbers  of  data  channels  are  given  in 
Table  8.  Three  axes  of  wind  speed  were  measured  with  an  orthogonal  impellor 
anemometer.  A  redundant  vertical  velocity  was  obtained  with  a  sonic  anemom¬ 
eter,  and  a  redundant  horizontal  wind  speed  was  measured  with  a  hot  film 

Table  8 


Devices 

Sampled  at 

10  Hz 

on  Main  Instrument  Boom  During 

Wind  Measurement 

Component  of  SUPERDUCK 

Number 

of 

Variable 

Channels 

Manufacturer/Model/Sensor  Type 

Three  axis 
wind  speed 

3 

R.  M.  Young 

Model  27005 

Gill  UVW  impellor  anemometer 

Vertical 
wind  speed 

1 

Campbell  Scientific 

Model  CA27T 

Sonic  anemometer 

Horizontal 
wind  speed 

1 

TSI 

Model  1610 

Hot  film  anemometer 

Air  temperature 
over  water 

1 

Campbell  Scientific 

Model  CA27T 

Chromel  constantan  thermocouple 

Air  temperature 
over  water 

1 

Atmospheric  Instrumentation  Research 
Model  FT-1A 

Platinum  wire  resistance 

Humidity 
over  water 

1 

Weathermeasure 

Model  5120-B 

Thin  film  capacitor 

Humidity 
over  water 

1 

Krypton 

Model  KH-20 

Optical  absorption  hygrometer 
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anemometer.  Redundant  air  temperatures  were  measured  with  a  fine  wire  resis¬ 
tance  thermometer  and  a  fine  wire  thermocouple  thermometer.  Redundant  humid¬ 
ity  measurements  were  obtained  with  a  thin  film  capacitor,  a  relative  humiditv 
sensor,  and  an  optical  absorption  hygrometer. 

Data  collection 

247.  Data  were  collected  using  two  programmable,  internally  storing 
data  loggeis  (Campbell  Scientific  Model  21X).  These  made  mean  and  covariance 
calculations  in  real  time  to  save  raw  data  storage  media.  The  data  loggers 
were  internally  clocked  and  had  battery  back-up  power  sources  to  avoid  data 
loss  when  conventional  60-Hz  power  was  interrupted. 

248.  All  mean  and  cross-correlated  variables  were  sampled  for  1  hr  to 
ensure  adequate  low-frequency  resolution.  About  every  2  days,  data  in  the 
data  logger  memory  buffers  were  transferred  in  binary  format  to  standard  audio 
cassette  tape.  Data  on  these  tapes  were  read  by  a  microcomputer  and  trans¬ 
ferred  via  phone  to  a  large  computer  on  the  NCSU  campus  for  archiving  and 
analysis . 

249.  Since  high  accuracy  is  required  for  turbulent  flux  measurements, 
raw  data  from  the  high-frequency  data  logger  were  collected  once  daily  for 

30  min.  This  allowed  for  detailed  error  checking  and  for  computation  of  spec¬ 
tra  of  various  turbulence  quantities.  The  data  loggers  had  digital  displays 
to  monitor  input  on  individual  data  channels.  Several  of  the  sensors  had 
analog  displays.  These  were  checked  frequently  to  ensure  that  all  sensors 
were  functioning  properly. 

250.  The  primary  variable  to  be  measured  is  the  Reynolds  stress.  In 
general,  this  is  a  tensor  with  nine  elements.  In  terms  of  velocity  compo¬ 
nents,  it  has  six  unique  terms  of  which  the  two  representing  vertical  fluxes 
of  horizontal  momentum,  i.e. ,  uTwT  and  v'w'  ,  are  most  important.  In  terms 

_  2  2  1/2 

of  horizontal  wind  speed,  this  is  s'w'  where  s'  =  (u'  +  v'  )  .In  this 

case,  the  horizontal  stress  vector  is  assumed  to  align  with  the  mean  horizon¬ 
tal  wind  vector.  With  the  instrumentation  described  here,  the  Reynolds  stress 
can  be  estimated  in  five  ways.  With  the  three-axis  impellor  anemometer  (see 
Table  8)  u'w'  and  v'w'  can  be  computed  directly.  Using  the  sonic  anemom¬ 
eter  for  w'  and  the  hot  film  anemometer  for  s'  ,  s'w'  can  be  computed 
directly.  With  w'  from  the  impellor  anemometer  and  s'  from  the  hot  film 
anemometer,  s’w'  can  be  computed  directly.  With  u'  and  v'  from  the 
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impellor  anenoneter  and  w'  from  the  sonic  anemometer,  u'w'  and  v'w' 
be  computed  directly.  From  the  raw  data  tapes,  spectra  of  horizontal  wind 
speed  can  be  computed  from  the  hot  film  anemometer  data  and,  from  this,  the 
inertial  subrange  method  (see  Part  IV)  can  be  used  to  estimate  u^  . 

251.  Three  additional,  but  less  reliable,  estimates  of  stress  can  also 
be  made.  Correlation  of  horizontal  wind  speed  from  the  top-boom  impellor  vane 
anemometer  with  w'  from  the  top-boom  vertical  impellor  anemometer  (both 
sampled  at  5  Hz,  Table  7)  provides  an  estimate  of  s'w'  but  at  somewhat  re¬ 
duced  high-frequency  resolution.  Spectra  of  u'  and  v'  from  raw  measure¬ 
ments  bv  the  three-axis  impellor  anemometer  can  be  used  to  estimate  u^  by 
the  inertial  subrange  method  if  corrections  are  applied  to  compensate  for 
reduced  high-frequency  response.  If  allowance  is  made  for  platform  blockage 
effects  at  the  13.7-m  level  on  the  pier-end  tower  (as  deduced  in  the  DUCK85 
experiment),  then  the  profile  method  can  be  used  to  estimate  u^  from  mean 
wind  speeds  measured  with  the  cup  anemometer  on  the  lower  boom,  any  of  the 
speed  sensors  on  the  main  boom,  and  the  impellor  vane  anemometer  on  the  top 
boom. 

252.  In  a  like  manner,  various  combinations  of  temperature  and  vertical 
velocity  or  humidity  and  vertical  velocity  sensors  can  be  used  to  make  esti¬ 
mates  of  Reynolds  fluxes  of  heat  and  humidity.  These  can  be  used  in  combi¬ 
nation  with  stress  measurements  to  deduce  stratification  effects  in  drag 
coefficient  formulae  based  on  Monin-Obukhov  similarity  theory  (Part  III). 

Data  analysis 

253.  When  the  data  from  the  experiment  have  been  compiled  and  edited 
for  completeness  and  accuracy  and  when  data  from  the  phases  of  SUPERDUCK  re¬ 
lating  to  auxiliary  parameters,  such  as  ocean-surface  current,  sea  state,  and 
over-land  air  temperature,  have  been  compiled,  then  the  model  equations  given 
in  this  report  can  be  tested  for  similarity  with  results  obtained  in  open- 
ocean  conditions.  A  report  describing  tests,  results,  and  recommended  changes 
to  the  model,  if  any,  will  follow. 
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PART  VI:  CONCLUSIONS 


254.  Drae  coefficients,  relating  wind  stress  to  a  set  of  mean  environ¬ 
mental  parameters,  are  highly  useful  and  practical  forms  of  closure  for  char¬ 
acterizing  wind  forcing  of  waves  and  currents  in  large-scale  numerical  models. 
A  survey  of  recent  literature  on  the  subject  indicates  that  drag  coefficients 
are  neither  constant  nor  simply  dependent  on  wind  speed,  but  are  dependent  on 
a  broader  range  of  parameters  intimately  related  to  the  subjects  of  air-sea 
interaction  and  physics  of  the  atmospheric  boundary  layer. 

255.  Research  has  shown  that  bulk  formulations  depending  on  gravity, 
wind  speed,  air-sea  temperature  differences,  and  air-sea  humidity  difference 
provide  reasonable  estimates  of  wind  stress  over  the  open  ocean  under  steady, 
horizontally  uniform  conditions.  However,  errors  in  excess  of  100  percent 
occur  under  unsteady  conditions  when  wind  and  sea  state  are  not  in  dynamic 
equilibrium.  This  suggests  that  in  shallow  and  nearshore  waters,  where  waves 
are  influenced  by  the  sea  bottom  as  well  as  the  wind  and  where  wave-induced 
currents  are  different  than  in  deep  water,  bulk  formulae  based  on  deep-ocean 
work  may  be  substantially  in  error.  Evidence  to  this  effect  has  been  given  in 
recent  literature. 

256.  Since  the  bulk  of  Corps  interests  are  in  waters  of  intermediate  to 
shallow  depths,  it  is  imperative  that  improved  models  be  developed  so  that 
proper  modeling  of  nearshore  processes  can  be  done.  Based  on  physical  reason¬ 
ing  and  research  reported  in  the  literature,  it  is  hypothesized  that  nearshore 
drag  coefficients  depend  on  water  depth,  distance  from  shore,  instantaneous 
sea  state  (at  least  significant  height,  peak  period,  and  primary  propagation 
direction),  wind  direction,  local  water-surface  speed  and  direction,  and  land- 
sea  air  temperature  differences  in  addition  to  parameters  used  in  deep  water. 

257.  To  test  this  hypothesis,  an  experiment  has  been  conducted  wherein 
these  parameters  and  wind  stress  are  measured.  Wind  stress  is  not  a  simple 
variable  to  measure.  Besides  bulk  coefficient  estimates,  there  are  three  com¬ 
monly  used  methods  to  determine  wind  stress.  They  are  the  profile,  direct, 
and  inertial  subrange  methods.  Of  these,  the  direct  method  is  to  be  preferred 
since  it  requires  the  fewest  assumptions. 

258.  There  are  constraints  on  the  applicability  of  any  model  or  the 
reliability  of  any  experiment  in  cases  of  extreme  unsteadiness  or  horizontal 
inhomogeneity.  The  latter  condition  may  be  of  considerable  importance  in. 
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waters  near  land  with  high  relief  or  strong  temperature  contrast.  In  such 
cases,  considerable  additional  research  is  necessary  to  resolve  an  accurate 
drag-coefficient  closure.  Near  land  of  low  relief,  and  during  high  wind  con¬ 
ditions,  where  gradients  of  air  density  have  less  importance,  the  parametric 
dependence  proposed  here  should  be  adequate. 

259.  An  experiment  was  conducted  in  two  phases  to  provide  insight  into 
the  parametric  dependency  proposed  above.  The  site  chosen  was  the  CERC’s  FRF, 
an  area  of  low  relief,  frequent  high  winds,  and  extensive  logistical  support. 
This  includes  a  pier  which  extends  500  m  seaward  of  the  beach  and  which  can  be 
used  for  over-water  stress  measurements.  The  first  phase  of  the  experiment 
took  place  during  what  is  known  as  the  DUCK85  experiment.  It  was  a  prelimi¬ 
nary  exploration  of  nearshore  uniformity  of  wind,  wind  steadiness,  and  ade¬ 
quacy  of  various  platforms  for  mean  wind  and  wind-stress  measurements.  Re¬ 
sults  indicated  that,  most  of  the  time,  unsteadiness  and  nonunif ormily  are  noc 
major  problems  at  this  site  and  that  nonuniformity  is  principally  due  to  land 
topographic  effects  for  strong  onshore  winds.  On  occasion  when  these  effects 
are  important,  they  can  be  recognized  and  compensated  for  since  variations  in 
the  nearshore  wind  field  appear  to  be  highly  systematic.  Results  indicated 
that  the  permanently  installed  FRF  anemometer  may  be  suspect  in  very  high 
winds.  Also,  air-flow  blockage  effects  may  be  substantial  at  lower  elevations 
on  the  seaward  end  of  the  FRF  pier. 

260.  The  second  phase  of  the  experiment  took  place  in  the  fall  of  1986 
during  what  is  known  as  the  SUPERDUCK  experiment.  In  this  experiment,  bulk 
meteorological  parameters  were  measured  along  with  wind  stress.  This  is  being 
determined  by  the  direct  and  inertial  subrange  methods  from  instruments 
mounted  on  the  pier-end  tower  at  elevations  high  enough  to  avoid  effects  of 
flow  distortion  by  the  pier.  Sufficient  redundancy  in  instrumentation  was 
employed  to  assure  reasonable  continuity  and  accuracy  of  results.  Simultane¬ 
ous  observations  of  sea  state,  water  currents,  and  other  relevant  parameters 
were  obtained  by  other  investigators.  These  parameters,  in  combination  with 
variables  measured  in  the  wind  experiment,  will  provide  a  basis  for  testing 
various  model  hypotheses. 

261.  Data  from  the  SUPERDUCK  experiment  is  believed  to  be  unique  in  its 
intent  and  completeness,  and  so  will  provide  valuable  insight  and  guidance  in 
the  evaluation  and  design  of  wind-stress  closures  for  nearshore  hydrodynamic 
models . 
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APPENDIX  A:  NOTATION 


(  i  Time  average  of  (  ) 

,  a.^  Dimensionless  similiaritv  model  constants 
A  Kolmogorov  constant  for  momentum 

A^  Kolmogorov  constant  for  humidity 

A_  Kolmogorov  constant  for  temperature 

,  Dimensionless  similarity  model  constants 

C  Water-wave  phase  speed 

Cp  Drag  coefficient 

Water-wave  phase  speed  associated  with  peak  spectral  density 

C  Phase  speed  at  which  maximum  air-sea  interaction  occurs 

max 

d  Water  depth 

D  Turbulent  kinetic  energy  gradient  diffusion  term 

2 

e'  Instantaneous  turbulent  kinetic  energy  per  unit  mass 

f  Cyclic  frequency 

f  Coriolis  parameter 

c 

g  Gravitational  acceleration 

G  A  function  in  the  simplified  Wallops  spectrum 

H  Characteristic  wave  height 

s 

i  Direction  index 

j  Direction  index 

k  Radian  wave  number 

k  Characteristic  size  of  physical  roughness  elements 

R 

L  Monin-Obukhov  length  scale 

m  A  function  in  the  simplified  Wallops  spectrum 

p  Total  atmospheric  pressure 

p  Mean  atmospheric  pressure 


AI 


p'  Fluctuating  atmosphe-ic  pressure 

a  Total  specific  humidity 

n  Mean  specific  humidity 

q'  Fluctuating  specific  humidity 

Mean  surface  specific  humidity 
R  Ideal  gas  constant  for  air 

s'  Resultant  of  fluctuating  horizontal  wind  speed  components 

S  Variance  spectral  density 

S  Humidity  variance  spectral  density 

q 

S^,  Temperature  variance  spectral  density 

t  Time 

T^  Averaging  time 

T^  Readjustment  time 

T  Temperature 

T  Mean  temperature 

T'  Fluctuating  temperature 

T  Mean  wave  period 

o 

T  Mean  surface  temperature 

o  r 

Tp  Spectral  peak  period 

T^  Virtual  temperature 

T^  Mean  virtual  temperature 

T^  Fluctuating  virtual  temperature 

u  Wind  speed  in  x-direction 

u'  Fluctuating  wind  speed  in  x-direction 

u'  Measured  fluctuating  wind  speed  in  x-direction 

TT1 

1/2 

u*  (T/o)  friction  velocity 

u.  Friction  velocity  estimated  from  measurements 

*m 


A2 


r 


w 

m 

x 


z 

a 

z 

o 

z 

r 

a 

9 

v 

s 


q 

f 

T 


9 

‘a 

9 

r 

< 


Mean  velocity  in  x-direction 

Particular  measured  mean  wind  speeds 

Reference  mean  wind  speed 

Surface  mean  wind  speed 

Fluctuating  wind  speed  in  v-direction 

Vertical  velocity 

Fluctuating  vertical  velocity 

Measured  fluctuating  vertical  velocity 

Primary  horizontal  coordinate 

Secondary  horizontal  coordinate 

Vertical  coordinate 

Particular  vertical  coordinates 

Vertical  elevation  at  which  U(z  )  -  C 

a' 

Surface  roughness  length  scale 
Reference  elevation 
Coefficient  in  Charnock's  model 
Phillips  constant 
Incomplete  gamma  function 
Boundary  layer  thickness 
Finite  increment  operator 

Viscous  dissipation  rate  of  turbulent  kinetic  energy 

Molecular  dissipation  rate  of  humidity  variance 

Molecular  dissipation  rate  of  temperature  variance 

Azimuth  direction  in  horizontal  plane 

Azimuth  of  mean  flow  at  air-boundary  interface 

Reference  horizontal  wind  direction 

von  Karman  constant 


A3 


> 


Wavelength 

a t e r  wavelength  associated  with  peak  spectral  density 
Latitude 

Dynamic  viscosity 

Kinematic  viscosity 

o  Total  air  density 

c  Mean  air  density 

o’  Fluctuating  air  density 

Standard  deviation  of  sea-surface  elevation 

r  Wind  stress 

i  Dimensionless  mean  shear 

m 

i  Dimensionless  mean  humidity  gradient 

q 

c  Dimensionless  mean  temperature  gradient 

V  Angular  deviation  from  level  of  measuring  system 

Dimensionless  integral  of  $ 
m  m 

Dimensionless  integral  of  <t> 

q  q 

;  Dimensionless  integral  of  4>T 

Radian  frequency 

Frequency  at  which  maximum  air-sea  interaction  occurs 

max 

Sea-surface  elevation  peak  frequency 

P 

:  Planetary  rotation  rate 


A4 


